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ABSTRACT 


An opening in the ionosphere may be created by the release of 
reactive chemicals that promote ion-electron recombination. The 
following sequence of events is initiated by a gas realease. As 
the gas expands, it is rapidly cooled. If the temperature of the 
vapor falls below the saturated vapor temperature, some of the . 
gas condenses. During the day, the condensate is re-vaporized by 
sunlight, but at night the condensed molecules remain as clusters 
that do not effectively interact chemically with the ionospheric 
plasma. 

When the vapor becomes sufficiently tenuous, it is reheated by 
collisions with the ambient atmosphere and its flow is then governed 
by diffusive expansion. As the injected gas becomes well mixed 
with the plasma, a "hole" Is created by chemical processes. In the 
case of diatomic hydrogen release, depression of the electron concen- 
trations is governed by the charge exchange reaction between oxygen 
ions and the hydrogen, a reaction that produces positive hydroxyl 
ions. This latter ion rapidly reacts with the electron gas to produce 
excited oxygen and hydrogen atoms. Enhanced airglow emissions will 
result from the transition of the excited atoms to lower energy 
states . 

Owing to the ion depletion, the remaining electrons are no longer 
effectively cooled by Coulomb collisions and the electron temperature 
in the depleted region rises sharply causing a thermal expansion of 
the plasma and a further reduction in the local plasma concentration. 
Replenishment of the plasma in the depleted region takes place in 
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part by plasma flow from the protonosphe're. Such a flow produces 
a tube of reduced plasma that extends into the magnetosphere and, 
for cases of extreme depletion, may extend into the conjugate 
ionosphere. 

Flights of space vehicles, such as the Space Shuttle, into the 
ionosphere may permit the controlled depletion of the F2-layer by the 
injection of selected gases. Estimation of the optimum modifica- 
tion experiment is based on calculations of temporal and spatial 
characteristics' of the ionospheric hole for various types and amounts 
of injected gas. Criteria for judging candidate gases include the 
amount of ionospheric plasma depleted per unit mass of released gas, 
the chemical reaction rates, and the resistance to photodissociation 
during the day and to condensation during the night. Based on these 
criteria, H 2 is recommended for gas release experiments. Water 
vapor and carbon dioxide are acceptable for daytime releases, but 
lose efficiency at night due to condensation. 

The amount of gas released partially determines the size and 
duration of the depletion. The release of 100 kg of diatomic hydro- 
gen into the ionosphere produces an F-layer depletion lasting over 
2 hours. 

Measurements of the temporal variation in the spatial structure 
of the ionospheric hole may provide information about the dynamic 
processes that maintain the ionosphere. These processes include 
chemical dissociation, diffusion, and photoionization. 

Artificially produced ionospheric holes may significantly affect 
radiowave propagation. During the early period of chemical 


iv 



modification, the hole may act as a lens that focuses or bends VHP 
or HP radio waves. In the later period, a few percent plasma density 
depression may extend over the entire geomagnetic field line. Such 
field-aligned irregularities may guide HP and VLF radio waves for 
i nterhemi spheri c propagati on . 
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Chapter I 


INTRODUCTION 

This research examines the ionospheric depletion produced by 
gases that react chemically in the F2-layer. The gases that promote 
0^ ion and electron recombination in the ionosphere are common in 
rocket exhausts. It is therefore not unexpected that the first 

* 

observations of chemical depression of the F2-layer plasma were made 
while' a rocket was firing at ionospheric altitudes. However, since 
expanding gases produce local depletion by simply displacing the 
ambient atmosphere, it was not always known that the ionospheric 
depletion by rocket exhaust was predominantly due to chemical 
effects. 

On February 17, 1959 around 1100 EST, Vanguard II penetrated • 
the F-region. At the time of the passage of the missile, ionosphere 
sounders located at Cape Canaveral and Grand Bahama Island recorded 
unusually disturbed ionograms. These ionograms were interpreted as 
resulting from a hole created' in the F-region by the missile's 
exhaust gases. This hole lasted for half an hour or more. It was 
theorized that the expansion of the hot exhaust gases had "punched" 
a hole in the ionospheric layer [Booker, 1961; Barnes, 1961]. 

A similar effect was observed by Faraday rotation measurements 
obtained from the Scout rocket (NASA ST-7/P-21) launched from Wallops 
Island, Virginia at 1238 EST on October 19, 1961. During the third 
stage burning period, when the rocket was at an altitude of 200 km. 
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Faraday rotation of signals emitted from the. rocket indicated an 
apparent reduction in ionospheric .electron densities in the wake- of 
the rocket [Jackson et al., 1962]-. Again it was speculated -that the 

4 

deficiency of electrons in the wake of the rocket was due perhaps 
to a displacement of the ambient electrons by the exhaust gases. 

The theory that the rocket produced a depletion of the F-layer 
solely through the displacement of the ionospheric plasma by the 
rocket exhaust gases was inadequate to explain the observed effects. 
The ,size and duration of the opening in the ionosphere indicated 
another process was responsible for the formation of this opening. 

It is now clear that chemical reaction between the rocket exhaust 
gases and the F-region ions was the main mechanism for production 
of the ionospheric hole. 

Chemicals were released into the E and F- regions by the Air 
Force during Project Firefly 1962 [Rosenberg, 1964]. Of the many 
controlled releases during this project, only one dealt with a 
chemical which might accelerate the ion and electron loss process. 
This was the injection of 5.44 kg CO 2 at an altitude of 234 km. 

With this gas, the loss of electrons in the F-region is promoted by 
the following chemical process: 

CO 2 + 0"^ ^ O 2 + CO 
Og + e 0 + 0 
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No measurable distrubance of the ionosphere could be attributed to 
the CO 2 real ease because it coincided with the passage of a large 
meteor that perturbed the ionospheric record and that masked the 
effects of the COg. In this author's opinion, failure of the COg 
release may also be attributed to the condensation of a significant 
portion of the COg* after the release. 

Except for the one CO^ release, chemical experiments in the 
ionosphere, during and after Project Firefly, were designed to 
explore processes other than acceleration of dissociative recombina- 
tion. Discussion of ionospheric modification by the release of 
electron attaching materials such as sulfur hexafluoride [Rosenberg, 
1964] or by the release of easily ionized material such as barium or 
strontium [Haser, 1967] is beyond the scope of the report. 

It was not until the late 1960's that the' reaction rates between 
the ionospheric ion constituents and possible neutral perturbing 
agents were measured accurately in the laboratory. An excellent 
summary of ion-molecule reaction rates is given by Ferguson [1973]. 

The interest in chemical depletion of the ionosphere was 
re-kindled in 1975 with the launch of the Skylab orbiting spacecraft 
by a Saturn V rocket. It was reported by Mendillo et al. [1975a, b] 
that a large and rapid decrease was produced in the ionospheric 
total electron content within 1000 km of the Saturn V burn. The large 
hole was attributed to chemical reactions Initiated fay the H 2 and H 2 O 
in the rocket exhaust. 

The development of the Space Shuttle and large payload rockets 
makes controlled releases of large quantities of reactive gases 
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practical. Thus, with increased knowledge of the F-region chemistry, 
and increased payload capabilities of space vehicles, the use of gas 
releases for scientific exploration of the ionosphere is more feasible 
than it was a decade ago. 

There are three general applications for experiments involving 
release of gases that promote ionospheric depletion. First, communi- 
cations may be disrupted or improved by modification of the F-layer. 
Second, the ionosphere may be used as a laboratory for the measurement 
of the perturbations resulting from the gas release. Observations- 
of the hole created in the ionosphere can provide important information 
about the chemical and diffusive processes that govern the F-region 
dynamics. Third, the depletion of the ionosphere can allow the pene- 
tration of radio signals that normally are reflected by the F-layer. 
These signals include cosmic radio waves of interest to radio 
astronomers and ducted VLF radio waves of interest to magnetospheric 
physicists. 

The organization of this report is illustrated by the flow 
diagram of Figure 1. Following this Introduction, four chapters deal 
with the theory necessary for a complete understanding of the response 
of the ionosphere to chemical releases. Chapter II is a discussion 
of chemicals that may be used to reduce the F-layer 0^ and electron 
concentrations. Chapter III is’ a theoretical treatment of the 
released gas dynamics from the time of injection to the time of 
thermal izat ion with the background atmosphere. Chapter IV presents 
the theory of diffusion necessary for modeling the flow of the 
released gas after thermal izati on has occurred. Also, the diffusive 
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CONCLUSIONS 


Fig. 1. Structure of Report. 










flow of the perturbed ionosphere is described in this chapter. 

Chapter V presents the energy balance equations used in modeling the 
disturbed thermal processes resulting from the release of reactive 
gases . 

The theory presented in Chapters II through V becomes the basis 
of a numerical model used to simulate the disturbances produced in 
the plasmasphere- In Chapter VI, the fluctuations in density and 
temperature caused by hydrogen releases into the F-layer are modeled. 
Chapter VII is a description of the extraordinary airglow emissions 
accompanying the depletion. Chapter VIII illustrates several uses 
of the depleted plasmasphere for guiding or focusing radio waves. 
Conclusions and recommendations for further research are presented 
in Chapter IX. 

This research has contributed to scientific knowledge by 
borrowing from the diverse scientific disciplines of chemistry, gas- 
dynamics, aeronomy, and electromagnetic theory to provide a detailed 
description of the ionospheric depletion experiment. Some specific 
contributions in the area of theoretical analysis of the controlled 
modification process are: 

1) Comparitive analysis of candidate chemicals for depletion 
of the F2-layer. Gases such as H 2 O, H 2 > CO 2 and NO 2 cause signifi- 
cant depletion of the ionosphere because of their high reaction rates 
with 0^ ions. 

2.) Theoretical study of condensation of expanding gas clouds. 
Condensation may limit the usefulness of released gases. After the 
gas is injected into the ionosphere, it expands and cools. Gases, 
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such as H 2 O, COg, and NOg, with high saturated vapor temperatures 
win partially condense during the early period of expansion. At 
night, the condensed droplets will not re-evaporate. The injected 
chemical in droplet form loses effectiveness as an ionospheric 
modifier, 

3) Treatment of the transition between free molecular expansion 
and diffusive expansion- of a released gas cloud. The injected gas 
cloud themalizes by collisions with the background atmosphere. 

Thermal izati on occurs within 10 seconds of the time of initial gas 
release. 

4) Development of a numerical model of the thermosphere for 
simulation of perturbations in density, temperature, and airglow 
produced by the release of reactive gases. Several phenomena 
accompanying the ionospheric depletion include changes in the electron 
temperatures and enhancement in the airglow emissions. During the day, 
the temperature of the electron gas may rise by 1000 K. The tempera- 
ture of the ions drops slightly in the perturbed regions. An expand- 
ing airglow ring is formed around the point of gas release. The 
intensity of the ring may be five times the unperturbed airglow 
Intensity. 

The ionospheric depletion experiment permits: 

1) Measurement of chemical and diffusion coefficients in the 
ionospheric "laboratory'*. Coefficients not available from ground- 
based laboratory measurements may be obtained. 
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2) High resolution observation of extragalactic radio sources 
by using the ionospheric hole as a lens. A refracting radiotelescope 
using the ionospheric lens may have a factor of ten better resolving 
power than existing radiotelescopes. 

3) Ducting of HF and VLF radio signals in the protonosphere 
above the depleted ionosphere. The depletion of the ionosphere 
produces ducts in the overlying protonosphere. Radio waves guided 
by this duct may be used for interhemi spherical communication or for 
controlled wave-particle interaction experiments. 

4) Measurement of ionospheric electric fields and neutral winds 
by observing the perturbed airglow emissions. The expanding airglow 
ring is a new type of chemical tracer that may be used at any time 

of day. 
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Chapter II 

CHEMICAL PROCESSES IN THE IONOSPHERE 

A knowledge of Ionospheric chemistry is necessary to predict 
the F-region plasma depletion caused by the release of certain gases. 
Fortunately, the F-region chemistry is well understood, thus permitting 
accurate modeling of artificially induced depletions. 

This chapter reviews the dominant chemical mechanisms in the 
region of interest and discusses in detail the plasma depletion 
caused by chemically increasing the 1on-electron frecombi nation rate. 

A. F-Reqion Photochemistry 

The plasma that constitutes the ionosphere is created by ioni- 
zation of gases in the neutral atmosphere by the sun. Photoioni- 
zation by shortwave solar extreme ultraviolet (EUV), i.e., radiation 
with wavelengths less than 1027 A, is the main mechanism for main- 
taining the daytime ionosphere. The photoionization of the neutral 
atmospheric constituents (i.e., mono and diatomic oxygen and diatomic 
nitrogen) proceeds according to these reactions 

0 + hv ^ 0^ + e~ 

O 2 + hv ^ Og + e“ (1 ) 

■ Ng + hv - 5 - Ng' + 

The rates at which these reactions occur depend on the frequency (v) 
and intensity of the solar extreme ultraviolet radiation, on the ioni- 
zation cross sections and on the local gas concentration. 
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Variations in soft x-ray and EUV fluxes produce corresponding 
variations in the F-region ionization level. The effects of the 
changes in solar radiation produced by solar flares and eclipses 
have been discussed by Garriott et al , [1967], Evans [1965], and 
others. 

Ionization by energetic charged particles is known as corpus- 
cular ionization. The energetic particles, guided by geomagnetic 
field lines, produce ionization by collisions with neutral particles. 
Corpuscular ionization may be important when there is little or no 
photoionization (i.e., in the night F-layer or in the polar winter 
F-layer). Electron impact ionization is an example of this ionization 
mechanism. 

0 + e"* 0"^ + 2e" - 13.618 eV 

O 2 + e"* ^02 + 2e“ - 12.07 eV (2) 

N 2 + e"* nJ + 2e" - 15.58 eV 

where e represents an energetic electron. 

The rates of these reactions depend on the concentration and ioni- 
zation cross section of the neutral gas, the latter being strongly 
dependent on the energy of the impacting electron. Natural sources of 
energetic electrons are photoelectrons produced in the photoionization 
process (1), high energy electrons produced during magnetic storms, and 
energetic electrons precipitating from the magnetosphere into the 
high latitude ionosphere. 
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Fluxes of plasma into and out of the protonosphere may have a 
significant effect on the ionosphere. The protonosphere is composed 
of ionized hydrogen atoms {i.e. protons). The neutral atmosphere 
in the protonosphere is composed of the light elements H and He. 

ions that diffuse to lower altitudes where monoatomic oxygen is 
abundant may undergo the charge exchange reaction 

+ 0 0**’ + H rate = 2,2 x lO'^^ [0] sec'^ 

[Banks and Kockarts, 1973] (3) 

to produce atomic oxygen ions. The rate of this reaction is pro- 
portional to the square root of the H ion temperature, T^. , and 
proportional to the atomic oxygen concentration [0]. Plasma flow 
from the ionosphere to the protonosphere, subject to the charge 
exchange reaction 

0+ + H + 0 rate = 2.5VlO“^^ T^^ [H] (4) 

[Banks and Kockarts, 1973] 

Here, T^ is the neutral temperature. 

The direction of the plasma flow depends on the time of day. At 

night the flow is downward with being converted to by the 

reaction in (3). The nighttime ionosphere is supported by this 

mechanism. During the daytime, an overabundance of 0^ is produced 

-|- Hh 

in the F-region. The 0 that diffuses upward is converted into H 
according to (4) and fills the protonospheric reservoir. 
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Fluctuations in the protonosphenc content are communicated to 
the ionosphere by flux' variations. Natural events, such as magnetic 
storms may deplete the protono spheric reservoir. A depleted protono- 
sphere will tend to produce a reduction in downward fluxes during the 
night and an increase in upward fluxes during the day. In either 
case, a reduced F-region ionosphere may result. 

Recombination of ions with electrons occurs either through 
radiative recombination 


+ e" -V X + h (5) 

or through dissociative recombination 

XY"^ + e" ^ X + Y {6)' 

Radiative recombination (5) involves the collision of an ion (X ) 

with an electron (e”). In the collision process a neutral atom is 

formed and energy is released as a photon. The radiative recombina- 

-12 3-1 -f- 

tion rate is about 10 cm sec , The 0 -electron recombination 
reaction is 0^ + e” -j- 0 + hv. 

Dissociative recombination (6) occurs when an electron collides 

with polyatomic ion which splits into two neutral particles. Because 

"7 3 -1 

of the large recombination rate for this process (~10 cm sec ), 
it is the principle mechanism for ion-electron recombination in most 
of the ionosphere. The dissociative recombination reactions for 0^ 

and Ng are 
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oj + e" 0 + 0 a(Op = 2.2 X 10“^ (300/1^)*^ ctn^ sec"^ • 

(7) 

nJ + e~ ^ N + N aCNg) = 3 x 10’"^ (300/Tg)*^ cm^ sec"^ 

The recombination rate decreases with electron temperature (T ) 

B 

-372 

because the electron-ion collision frequency varies as • 

-I- 

The loss of 0 in the ionosphere occurs principally through 
charge- exchange reactions with the neutral polyatomic constituents 
of the atmosphere. These ion-neutral reactions are 

+ Og oj + 0 + 1.55 eV yiOg) = 2.0 x 10"^^ cm^ sec"^ 

( 8 ) 

o’*" + N2 ^ NO'*’ + N + 1.1 eV+ y(N 2) = 1^2 X' 1Q-12(^)V2 ^^^3 ^^^-1 
where T is the neutral temperature. 

The charge exchange rates, y( 02^ higher than the 

radiative recombination rates for 0 but are significantly lower 
than the dissociative recombination rates of (7). The ion charge 

* 4 " + * 

exchange products given by (8), NO and 02> are mainly lost 
through dissociative recombination 

NO^ 4 e" -> N + 0 , a(NO^) = 4.1 x 10~^ cm^ sec”^ 

{9) 

O2 + e" 0 + 0 a(02.) 2.2 x 10“^ cm^ sec"^ 
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TABLE 1. Perturbing gas reactions in the ionosphere 
(See Appendix A) 


c> 



Weakly Reactive 

Moderately Reactive 

Strongly Reactive 




NjO 

COg 

HOg 

Hg 

HgO 

Ion-n\oleciiTe 
Reaction with 
o" 

0^ + N, * no'*' + N . 

+1.11 eV. 

o'*’ + 0, -r O,'*’ + 0 

+^1.55 eV. 

0’’’ + N,0 -r H,0'^ + 0 
+“^.72 eV. 

+ NO’^ + NO 
+ 5.9 eV. 

O'^ + CO, -» 0,"^ + CO 
+‘^1.18 eV. 

O'^ + NO, NO,’’’ + 0 
+ 3784 eV. 

0* + H, OK’*’ + H 
+ .37 eV. 

0* + H,0 ->■ H,0* + 0 
^ +“^1.01 eV. 

(rate: cm^/sec) 

1.3[-12] 

2. [-11] 

2.2C-10] 

1.2E-9] 

1.6E-9] 

2. [-9] 

2.3E-9] 

Ion-molecule 
Reaction with 
n" 


N"^ + 0, ■*. 0, + N 

+^2.4B eV. 

•* ttO’^ + N 
+ 6.88 eV. 


n’’ + CO, CO,’’ + N 
+ :78 eV. 

+ NO + NO’^ + N 
+ 5.3 eV. 

«■* + H, -<• HH’* + H 
+ .66 eV. 

n’* + H,0 H,0’* + N 

“ +n.94 eV. 

(rate: cmVsec) 


7. [-12] 


1.3[-9] 

8.E-10] 

7.[-10] 

2.6E-9] 

Other 

Ion-molecule 




CQ,^ + 0 0,’’ + CO 

^ +“^1.33 eV. 

D,’’ + HO, -r HO,’’ + 0, 
^ * 2t29 eVT 

oh’* + H, H + H,o’’ 
^ + 1.82'^eV. 

H,0’’ + H,0 + H,0’* + OH 
^ +•'.43 eV. 

Reactions 
(rate: cm^/sec) 




1- 0* + CO, 
+ .15 evr 

2.6C-10] 

6.6E-10] 

AH’* + 0, OH + 0,’* 

+ 1-12‘^eV. 

2.E-10] 

1.7E-9] 

Primary Electron- 
Ion Recombination 
Reactions 

HO + e N + 0 

+ 2.75 eV. 

oY + e" r. 0 + 0 

+ 6.95 eV. 

NO + tt •* N + 0 

+ 2.7b eV. 

0,’’ + e“ -r 0 + 0 

+6.95 eV. 

NO,’’’ + e" * NO + 0 

+ 4.81 eV. 

OH"* + e" -r 0 + K 

+ 8.77 eV. 

HyO"* + e" -r OH + H 

+ 7.46 eV. 

(rate cm /sec) 

4-U-7] 

2. [-7] 

4.1[-7] 

2.[-7] 

4. [-7] 

l.[-7] 

3.C-7] 

Other Electron- 
Ion Recombination 
Reactions 


no’’ + e‘ H + 0 

+ 2.75 eV. 

+ e" ■> NO + N 

+7.93 eV. 

COjj^ + e **■ CO + 0 

+ 8.28 eV. 

lO’’ + e" ■* H + 0 

+ 2.75 eV. 

W’" + e~ -► N + H 

+ 9.41 eV. 

H,0’’ + e" * H,0 + H 

+V.04 eV. 

(rate: cinvsec) 


4.1[-7] 

4. [-7] 

3.8[-7] 

4.1E-7] 

l.E-7] 

l.E-6] 

Photo- 

dissociation 

Reactions 

■ 

Oj + hv V 0 + 0 

HgO + hv ->■ HO + H 

-♦ Hg + 0 

CO 2 + hu -»■ CO + 0 

n 

Hg + hv -«* H + H 

in 

(rate: sec'^) 


4. [-6] 

2.[-6] 

1.9E-7] 


» 0 



NOTATION: [-x] represents 10"'^ 

























Variations in temperature and composition of the neutral atmo- 
sphere produce changes in the ion composition. The neutral atmosphere 
density increases as altitude decreases. At some altitude below the 
F2 layer, the 0^ and N2 densities become sufficiently high so that 
the charge exchange with Eq. (8), produces polyatomic ion 

densities that are greater than the^ 0^ densities. This region is 
called the FI layer. 

By considering the dependence of ion loss on temperature and 
composition several means of artificial modification of the ionosphere 

may be suggested. Heating of the F-region will promote increased ion 

r 

densities due to the reduced recombination rate but will also tend 

■ 4 * * 

to lower ion densities due to thermal expansion. The loss of 0 
through charge exchange can be promoted by the introduction of 
relatively small quantities of foreign gases with high affinities 
with o'*'. This latter type of chemical modification is discussed in 
the next section. 

B. Chemical Perturbation of the Ionosphere 

In Section A, it is stated that the chemical loss is controlled by 
charge exchange with and Og. As shown in Table 1, numerous gases 
react faster with 0^ than do 0^ or Consider the ion molecule 

reaction with diatomic hydrogen- 

Hg + 0"^ + H yCHg) = 2x10"^ cm^ sec‘J (10) 


15 


LIFETIME (seconds) 



[X]/[N2](cm“^/cm’'^) 


* 1 * 

Fig. 1. Change in 0 lifetime with the introduction of chemically reactive gases into the 

ionosphere. Ambient Np concentration is 3 x 10^ cm“3. Op concentration is 10' cm"^. 
Undisturbed O’*" concentration is 2 x 10^ cm“3. 



Because is three orders of magnitude more reactive than N^, the 
release of a small amount of diatomic hydrogen will rapidly reduce 
the local 0^ population. 

Significant decrease in average lifetime of an 0^ ion occurs 

when 

Km 

DO . 

[N2] 


where 


[X] = concentration of released gas, 
CN 2 ] = concentration of Ng, 

Kj^^ = (0^, 1 ^ 2 ) reaction rate, and 

= (0^, released gas) reaction rate. 


For the example of the release of the inequality becomes 


[H2] 

[iy ' 


1 

1540 


The curves in Figure 1 show the reduction in O'** lifetime for 
different values of the ratio of released gas to nitrogen conentra- 
tions. It appears that substantial reduction in the ionospheric 0^ 
concentrations can be caused by the release of modest amounts of 
certain gases. 

In the unperturbed F2-layer, the electron loss coefficient is 
essentially the 0^ recombination coefficient. This is because the 
ion-molecule charge exchange reactions with N 2 and O 2 (8) are 
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slow compared with the dissociative recombination of the resulting 
molecular ions NO^ and 0^ (9). 

This is not necessarily true in the perturbed ionosphere. The 
reaction rates between 0^ and certain gases are comparable to 
polyatomic ion dissociative recombination reaction rates. As an 
example, consider the dissociative recombination of the OH^ pro- 
duced by (10). 

OH^ + e~ 0 + H a(OH^) 1 x 10^ cm^ sec"^ (11) 

J 

* 1 * 

If is released into the ionosphere, rapid depletion of 0 
occurs and large quantities of OH are produced. This OH then 
reacts with the electron gas to produce neutral atoms. There is a- 
significant time delay between the 0^ depletion and the electron 
depletion when gases with high charge exchange reaction rates are 
injected into the F2 region ionosphere. 

Examination of Table 1 indicates that the perturbed ionosphere 
chemical processes are complicated by ion molecule reactions with 
the foreign polyatomic ions. Again consider OH^, the product of the 
0^-H2 charge exchange reaction. OH'*’ may react with the injected 
gas, 

oh'*' + HgO'*' + H + 1.82 eV y(0H'*',H2) = 1.4 x 10"^ cm^ sec"^ 

( 12 ) 

or react with the atmospheric Og 
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OH'^ + O 2 oJ + oh + 1 .1 eV 


YCOH'^jOg) = 2 X 10"^° cm^ sec"^ 

(13) 


The polyatomic Ions produced in these reactions have dissociative 
recombination reaction rates which are comparable to the recombination 
rates of their parent ion, OH^. The production of secondary poly- 
atomic ions can be neglected in the modeling of chemical ionospheric 
perturbation because of the rapid dissociation of these ions. 

During the daytime, photodissociation of the injected gases 
must be considered. If the gas is photodissociated into non-reactive 
components it is not useful as an ionospheric modifier. Consider 
NOg which has the highest photodissociation coefficient in Table 1, . 
A comparison of its charge exchange time constant, = V([0 ] 

1.2 X 10"®), with its solar photodissociation time constant, Xp^ = 
1/(8, 1 X 10"^), shows that it is not an effective gas if Xpp < 
or if [0*^] < 6,8 X 10® cm"^. This inequality usually is true 
throughout the F-region ionosphere, NO 2 must be discounted as an 
efficient daytime perturbation gas. 

Gases, such as H 2 , with negligible solar photodissociation 
properties are useful as daytime injectants. 

An important consideration concerning chemical modification is 
the depletion efficiency per unit mass of the released gases. 

Diatomic injectants, such as H 2 , remove only one 0 -electron 
pair. Triatomic injectants may remove two 0^-electron pairs per 
released molecule'. Water vapor is an example of this type of gas. 
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H 2 O is dissociated into monoatomic hydrogen and the hydroxyl 
radical upon reaction with the o'*’-el ectron pair. 


H 2 O + o'*" + e" ^ H20'^ + e" + 0 ^ oh + H + 0 


(T4) 


A second 0^-electron pair is removed by reaction with OH. 


I O 2 . + H 

OH + 0^ + e~ j or 

0H'^+ 0 


+G ^0+0+H 


05) 


The rate of this reaction has not been measured. 

Nitrogen dioxide is an inefficient tri -atomic modified because 
it chemically removes only one 0^-electron pair. Charge exchange 
and dissociative recombination reactions involving 0^, electrons, 
and NO 2 proceed rapidly. 


NO 2 + NO 2 + O' 
NO 2 + e' NO + 0 


y(N02) = 1.6 X 10'^ cm^ sec 
a{N0p = 4 X 10”^ cm^ sec"^ 


However, charge exchange with the produced nitric oxide is 
negligible. 

I NO"^ + O 

y(N 0) < 10"^ ^.cm^ sec"^ 

N + O 2 
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The ratio of the gas molecular weight to the number of 0^ 
electron pairs annihilated by several gases is shown in the follow- 
ing table. 

Injected Substance N 2 O 2 N 2 O CO 2 NO 2 

Mol. Wt./(0^ - e”) pair 28 32 44 44 46 2 9 

The highest depletion efficiency per unit mass is obtained by the 
use of diatomic hydrogen and water vapor. 

C. Summary 

The feasibility of artificial modification of the ionosphere has 

been examined with respect to the F-region chemical mechanisms. 

* 

Depletion of the ion population may result from increased chemical 
loss processes. The release of reactive gases is an efficient 
method of increasing the ion-electron recombination in the ionosphere. 
When the ion-molecule reaction rates, photodissociation rates and 
mass depletion efficiency are considered for various gases, H 2 and 
H 2 O are found to be the most useful gases for temporary depletion of 
the ionosphere. 
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Chapter III 

THE SUDDEN RELEASE OF A GAS INTO A TENUOUS ATMOSPHERE 

Chemical modification of the ionosphere can be promoted by the 
release of a dense cloud into the ambient atmosphere. The cloud expands 
rapidly, driven by its thermal energy which is converted into kinetic 
energy, causing a rapid drop in temperature and, possibly, gas 

f 

condensation. When the gas cloud expands to an extremely rarified 
state, condensation ceases because of the lack of collisions between 
the free molecules and the condensation centers. At this stage, any 
particles that have condensed will begin to re-evaporate. 

After being cooled by expansion, the gas is re-heated by collisions 
with the ambient atmosphere. During thermal izati on with the atmosphere, 
the gas cloud temperature approaches the atmospheric temperature and 
expansion is dominated by collisional, diffusive processes. 

This chapter is a discussion of the rapid expansion of the 
released cloud between the time of release and the time of diffusively 
controlled expansion. Chapter IV is an analysis of the later epoch 
of isothermal, diffusive expansion. 

, Gas releases from vehicles passing through the ionosphere can 
be classified into two types. The first type is the spherical flow 
produced by an instantaneous release of gas at a given location. 

Examples of this type of source model are releases involving spherical 
detonations of an explosive gas mixture or pulsed gas jet releases. 

The second type of release is the continuous release of gas over a 
period of time, such as is produced by a rocket exhaust plume. The 


VW 
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The analysis of both types of releases is similar. Only the spheri- 
cal release will be considered in detail. 

The notation is consistent throughout this chapter. Number 
density, velocity, and temperature at time t^ are denoted as n. , 

, and T. respectively. The subscript i Ms constant for all 
variables at a fixed time. 

A. Spherical Expansions 

The problem of spherical gas expansion has been studied for a 
number of years. Point source shape preserving (similarity) solutions 
for expansion into a vacuum have been derived by Mirels and Mullen 
[1963], Grei finger and Cole [1965], and Hubbard [1967]. These 
solutions are based on the assumption of continuum, inviscid flow. 

If the mean free path inside the gas cloud is larger than the cloud ' 
diameter (i.e. Knudson number greater than 1), collisionless solutions 
apply. Gas cloud expansions into a vacuum are considered from the 
free molecular point of view by Mirels and Mullen [1963], and' 
Narasimha [1962]. For models of spherical expansions involving a 
wide range of Knudson numbers, representing continuum to collision- 
less regimes, kinetic theory is employed. Solutions of the Boltzmann 
equation for expansions into a vacuum are produced by Edwards and 
Cheng [1966, 1967], Freeman [1967], Narasimha [1967], Grundy and 
Thomas [1969], Grundy [1969], Thomas [1971], Freeman and Thomas 
[1969], and Cooper and Bienkowski [1965]. A unification of the 
continuum and kinetic theory is presented by Grundy [1969]. 
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Previous papers dealing with spherical expansions into an ambient 
atmosphere are similarly classified according to the method of 
solution. Continuum, inviscid analysis concerning expansions into 
a dense atmosphere are considered by Taylor [1946], Ness and Ambrosiani 
[1968], Cheng and Kirsch [1969], Hoffman and Bert [1974], and 
Baxter and Linson [1976]. Collisionless expansions are discussed by 
Bienkowski [1964]. Kinetic theory solutions for expansions into a 
rarified atmosphere are derived by Erode and Enstrora [1972], Brook 
and Hamel [1972], and Baum [1973, 1974]. A discussion of the 
expansion solutions appropriate for the gas release into ionospheric 
regions constitutes a large portion of this chapter. 

The expansion and simultaneous cooling of a gas cloud is 
easily explained from a molecular point of view. Consider a collec- 
tion of molecules confined in a spherical container with radius r^. 

If these molecules are in translational equilibrium, they will have 
a Maxwellian velocity distribution. Assumirig the gas temperature 
is Tq, the gas particle distribution in space and time is 



where x is the position vector, v the velocity vector, m the 
released gas molecular mass, n^ the number density inside the sphere, 
and k Boltzmann's constant. A velocity-position plot of this 
distribution is shown in Figure 1. 
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Fig. 1. Distribution funct ion for a spherical 
gas cloud, = y2kT^/m 

Macroscopic properties of the gas are foufTd by calculating velocity 
moments of the distribution. The density of the gas is related to the 
"0th” moment of the distribution as follows: 

n(x) = fff f(x,v)d^ V (2) 

00 

where n(x) is the number density of the gas at position x. The 
mean velocity and the temperature are calculated by the first and 
second moments of the particle distribution. 
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00 



/// f(x,v) vd^v 


n(x) 


( 3 ) 


3kT 
0 

m 


f(x,v)v^d\ 

-CO 


- uht) 


(4) 


Examination of Fig. 1 shows that 1) n(x) is constant for |x[ < r^ 
and equal to zero for |xj > r^, 2) u(x) = 0 everywhere, and 3) 
is constant inside the gas sphere. 

If the walls of the spherical container are removed, the gas will 
expand. The change in the distribution function can be quickly 
calculated if one assumes no collisions. In collisionless expansions, 
the flow of the gas is determined by the inertia of each molecule. 
Those molecules with a positive velocity propagate in the positive 
direction at a constant speed. Similarly, the molecules with a 
negative velocity move in the negative direction and those with zero 
velocity remain stationary. Figure 2 shows the collisionless distri- 
bution function at a finite time, t, after the start of the 

I- , 

expansion. The distribution remains constant along particle paths. 
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Fig. 2 . Distribution function for a 

spherical gas cloud undergoing 
collision! ess expansion. 

Examination of the moments of the distribution in Figure 2 shows 
a change in the radial structure of the density, velocity, and temper- 
ature. At a fixed position, x-| , the range of velocities, Av^, 
attained by the molecules is finite. The temperature is a measure of . 
the "width" of the velocity distribution at each position. Since the 
velocity distribution in Figure 2 is "narrower" than the one in 
Figure 1, decrease in temperature is produced by the expansion. 

After a long time, the distribution becomes a narrow strip (Fig. 3). 



Fig. 3. Limiting form of the collision- 
less velocity distribution. 
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Because the width of the velocity range, Av^^, approaches zero, 
the density and the temperature are vanishing quantities. The 
velocity at each point in space asymptotically approaches u = x/t 
where t is the time since beginning of the expansion. 

The inclusion of collisions in the expansion model increases 
the width of the velocity distribution at a point in space. Thus, 
the temperature in collisional expansion remains higher than the 
temperature in collisionless expansions. These qualitative arguments 
about the flow of gas during a spherical expansion can’ be expressed 
in quantitative terms by the solution of the Boltzmann equation. 

B. The Boltzmann Equation 

The evolution of the particle' distribution is found by solving 
the Boltzmann equation - 


M. 

at 




V 


Vf = 



coll 


( 5 ) 


where f is the distribution function, v, ' x, and t the velocity 
space and time coordinates respectively, and represents 

the change in the distribution function due to collisions. The 
collisionless solution to this equation is found by setting the 
collision term to zero. The Boltzmann equation for free molecular 
flow is 
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The solution is a simple function of the initial distribution. 

f^{x,v,t) = f|^{x- vt,v,o) (6) 

The subscript "m" refers to the free molecular solution. This 
result is a mathematical representation of the invariance of the 
distribution along particle paths as mentioned in the previous 
section. 

Using the initial Maxwellian distribution in (1) and the collision- 
less solution to the Boltzmann equation in (6), the moment equations 
(2), (3), and (4) are evaluated. The mathematical details are given 
in Bienkowski (1964). The resulting solutions for density, velocity, 
and temperature are expressed as 
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where is the Initial gas number density, the initial mean 
thermal speed = y2kT^/m, the initial temperature. The radial 
and temporal parameters have been nondimensionalized as follows 


r 




At the center of the spherical expansion, r = 0, these 
equations reduce to 


n/n = 


erf - - 

T 


TTT 


1/2 


exp 


. 1 


(10) 


u(o,t) = 0 


(11) 


T(o,t) ^ T^) exp (-1/t^) 

^0 n/n 


(12) 


The inclusion of the collision term in the Boltzmann equation is 
necessary for realistic modeling of gas expansions. Bhatnagar, Gross, 
and Krook (1954) have proposed the following approximation to the 
collision term 


3t) ■ “("<^0 ■ 


(13) 


where n is the gas number density, w is the gas collision frequency 
and f^ is the Maxwellian distribution expressed in terms of the 
local mean velocity and the local temperature. The BGK model 
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represents the relaxation of a non-Maxwel lian distribution to the 
Maxwellian distribution The gas relaxation time is 1/w. Note 

that the three variables in the BGK model {u, n, and f^) are 
functions of n, u, and T, which are functions of f, the required 
distribution function. 

The Krook equation is the Boltzmann equation with the BGK 
collision model . 


+ V • vf = w(nf^ - f ) (14) 

An expansion of f in terms of f^ is valid for small departures 
from equilibrium. Using this method of expansion, the Chapman- 
Enskog solution of the Krook equation gives the Navier-Stokes 
equations of macroscopic gas dynamics [Vincenti and Kruger, 1965]. 

In the case of inviscid, continuum flow, the Navier-Stokes equations 
reduce to the equations of motion for unsteady flows. 

Equation of Continuity 


+ V • (n7) = 0 


Equation of Motion 


= - - vp 

Dt n " 


First Law of Thermodynamics (Conservation of Energy) 


m_2_ /h J- ^ ^ 

nipt (h + jv ) 


(15) 


(16) 


( 17 ) 
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s^> P is the gas pressure, 

h the enthalpy per unit masstaad the mean square velocity. 

Under the assumptions of a perfect gas that is expanding isen- 
tropically (i.e., with no heat input into the gas) the pressure can 
be related to the density by solution of the first law of thermo- 
dynamics [Morse, 1967]. 


p = const • n^ (18) 

where y is the ratio of the specific heats. 

Making- this substitution into the equations of motion yields the 
following equations, in spherical coordinate form 


3n . 3nv . «nv ^ 
9t 8r 


(19) 


8v 


.8V 


^ = Y const 


,(y-2) ^ 


3r 


( 20 ) 


where v is the radial expansion velocity and r is the radial 
coordinate. 

As pointed out by Stanynkovich [1960], the density gradient on 
the right side of (20) becomes negligible after long times. With 
this assumption, (20) is solved to give v = r/t which was also the 
asymptotic velocity for free molecular flows. The substitution of 
the velocity expression into the continuity equation (19) leads to 
solutions of the form 
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( 21 )- 


nt^ = A(r/t) 


These long time solutions are called "self-similar" because v 

3 

and nt are functions only of r/t. 

Using the perfect gas law, p = nkT, the time variations in 
temperature and density at r = 0 are calculated using (18) and (21) 



( 22 ) 


^ _ 


n 


o 



(23) 


where t^ is the equivalent starting time of the self-similar 

expansion. Zeldovich and Raizer [1967] state that the characteristic 

time scale, t^, is calculated from the long time velocity solution 

t = r /c . If the time reference is shifted so that the gas expan- 
0 0 0 

Sion starts at t = 0, the temporal variation in temperature and 
density for the center of the continuum expansion are approximately 
given by 





3(y-l) 


(24) 


(2B) 

i 
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Fig. 4. Density and temperature variations in a spherical gas expansion. 



where the nondimensionallzed time unity, t, has been used. 

A comparison of the time variation of the density and of the 
temperature for free molecular flow (10) and (12), and for inviscid 
continuum flow, (24). and (25), is illustrated in Figure 4. The 
actual T - n path of the expanding gas cloud lies somewhere between 
these curves. Initially, the large central density of the gas is - 
in translational equilibrium due to collisions so that the T - n 
variation is close to the inviscid, continuum solution. As the central 
gas density decreases, collisions become rare and the free molecular 
solution is approached. 

The analytic expressions derived in this section are representa- 
tive of the temperature and density variations during the gas 
expansion. Because the gas density drops with increased radial dis- 
tance from the center of the expanding cloud, the outer portions of 
the cloud may be in a state of free molecular flow while the inner 
core gas density obeys the ordinary laws of continuum gas dynamics. 

The outer, free molecular shell of the cloud is only subject to forces 
produced by collisions with the ambient atmosphere into which the gas 
is expanding. 

When molecules in the dense inner core of the cloud collide 
with one another, they may adhere to each other. If the temperature 
of the cloud is low, these clusters are thermodynamical ly "stable and 
the expanding cloud condenses into a fog. A discussion of conden- 
sation of expanding vapor clouds is found in the next section. 
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C. Condensation of Expanding Vapor Clouds 


This section presents a discussion of the condensation resulting 

from the rapid cooling of an expanding gas. Each step in the 

condensation and subsequent re-evaporation process is described. An 

analytic example appropriate for a spherical water vapor release with 

25 3 

an initial temperature of 2000 K and an initial density of 5 x 10 /m 
is presented. (The physical data for HgO as well as for are. 
given in Appendix B). 

As a gas cloud expands, the internal pressure and temperature 
drop. At some point in the expansion, at time t-j , the gas pressure 
curve may intersect the saturated vapor pressure curve for the solid- 
liquid or liquid- vapor condensation system. The saturated vapor 
pressure is defined as the vapor pressure of a condensation system in 
equilibrium. The pressure of saturated vapor is governed by the 
Clausius-Clapeyron equation 

n 

TT'TfvVT 1969] (26) 

where is the saturated vapor pressure, T the gas temperature, 

U the latent heat of vaporization, v^ the gas volume and v^ 

the liquid or solid volume. (All volumes are the reciprocal of the 

local densities). Assuming 1) that the perfect gas law, p^v = RT, 

holds, 2) that the condensate density is much greater than the vapor 

density, and 3) that the U is constant, (26) leads to the following 

relation between saturated vapor volume, v , and temperature, T, 

, y 
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Vg = (27) 

where B is a constant. The saturated vapor temperature-vol ume 
(T-v) curve for H 2 O based on measured data [Handbook of Chemistry 
and Physics, 1959] closely matches the form of Eq. (27). The T^v 
curves representing both solid-vapor and liquid- vapor systems for 
H 2 O are shown in Figure 5. 

When the gas cloud expands beyond the state of saturation, it 
becomes super-saturated and condensation rapidly forms around any 
centers of inhomogeneity such as ions or dust particles. This is 
called heterogeneous nucleation. If there are no impurities in the 
expanding system to act as condensation centers, the gas expands in 
a super-saturated state until self-nucleation begins. This process, 
called homogeneous nucleation, involves embryos composed of clusters 
of gas molecules. An excellent introductory discussion of homogeneous 
nucleation is presented by McDonald [1962, 1963], reprinted in the 
text by Abraham [1974]. For a comprehensive treatment of both hetero- 
geneous and homogeneous nucleation see Zettlemoyer [1969]. 

Fundamental to the understanding of homogeneous nucleation is the 
concept of critical cluster size. The saturated vapor pressure of a 
nucleated cluster is dependent on the radius of the cluster. Clusters 
will small radii of curvature have higher surface vapor pressure than 
those with large radii of curvature. Consequently, a small cluster 
in a solid-vapor system may evaporate while a large cluster in the 
same environment may grow. The critical cluster size with respect 
to a given situation is the transition size between evaporating and 
growing clusters. The critical cluster radius is given by 
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Fig. 5. 




Fig. 6. Etnperical critical supersaturation ratio curve. 



O 

«• 



TIME 


Fig. 7. The supersatu- 
ration ratio and nucle- 
ation rate as a function 
of time. 
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(28) 



where a is the. surface tension, Vg the volume per molecule in the 
cluster, T and p the vapor temperature and pressure, respectively, 
and p^ the flat surface vapor pressure. 

Significant homogeneous nucleation occurs at critical supe‘^- 
saturation. The thermodynamic state (i.e., pressure, temperature, 
density) at the point of critical super-saturation is found by calcu- 
lating the point at which the homogeneous nucleation rate becomes one 

3 

embryo per cm per second. The homogeneous nucleation rate of 
droplets from a super-saturated system was derived by Frenkel [1946] 



where a is the surface tension, v the molecular volume in the 
liquid, N Avogadro's number, p the supersaturation pressure, and 
r* the critical nucleus radius. The degree of super-saturation is 
defined as S = p/p^. The point of critical super-saturation is 
insensitive to the definition of "significant" nucleation because of 

O 1 

the exponential nature of (29) (i.e., setting J = 1 cm" sec" or 
4 -3 -1 

10 cm sec leads to approximately 'the same result for degree of 
critical super-saturation, S = p /p ). In the case of water vapor, 
comparison of the theoretical variation of with T with cloud 
chamber measurements shows good agreement only at temperatures above 
260 K [Madonna et al., 1961]. Consequently, measured values for the 
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critical super-saturation pressure ratio for H^O (Figure 6) are 
used in numerical examples. Knowing both the saturated vapor pressure 
and the critical super-saturation ratio as function of temperature, 
the critical T-v curve for H^O is derived (Figure 5). 

After critical super-saturation is achieved, condensation centers 
are formed and the gas rapidly nucleates until a state of equilibrium 

9 

saturation is reached. The transition time from critical super- 
saturation to equilibrium saturation is on the order of a microsecond. 
The total number of condensation centers formed in this transition 
remains fixed as the gas continues to expand. This is schematically 
illustrated in Figure 7 adapted from Feder et al. [1966]. Simul- 
taneously, plotted versus time, are saturation level, S, and the 
homogeneous nucleation rate J. As the degree of super-saturation 
increases with time, the nucleation rate remains low until critical 
super^saturation is reached. At this time, t = t 2 , a spike in the 
nucleation rate is produced. This spike is quenched by the sudden 
lowering of the saturation level. 

The number of condensation centers per unit volume is found by 
integration of the nucleation rate over time. 

to 

N = / J(t)dt == J At, (30) 

Q IHuA 

P 

where J is the maximum nucleation rate and At, the effective 
max 

duration of the nucleation maximum. The cluster density for water 

vapor after critical super-saturation is reached, can vary between 
12 1 R R 

10 to 10 clusters/cm [Feder et al. 1966], depending on the 
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expansion rate, initial temperature and other thermodynamic factors. 

For the water vapor example of this section N = 6.3 x 10^^. 

Defining a growth factor, g(t,t') as the number of molecules 
at time t in a cluster born from a nucleus which appears at time t', 
an expression for the amount of material condensed is written 

-t 

x{t) = / J{f)g{t,f)dtVn (31} 

0 

where n is the fluid density and the degree of condensation, x(t) 
is the ratio of the molecules in the condensed state to the total 
number of released molecules. The relation between the vapor density 
and the fluid density is n^ = (l-x)n. Because the nucleation embryos 
are -almost instantaneously formed at the time of critical super- 
saturation, the nucleation rate is taken to be 

J(t) = .6(t-t2) = N6(t-t2) (32) 

where t 2 is the time of critical super-saturation, 6(*) the 
impulse symbol and At the equivalent width of the nucleation spike. 
Equation (31) becomes 

x(t) = N/n g(t,t 2 ) = V g(t,t 2 ) (33) 

Assuming that the embryos expand with the free gas molecules, the 
ratio of cluster density to fluid density remains unchanged. 

W/n = V = constant 
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Once the gas is in equilibrium saturation, it obeys the follow- 
ing adiabatic equation 

[c^(l-x) + CgX] dt + Rt(l-x) dv/v - [U - (c 2 -c.|)T] dx = 0 (34) 

where c^ is the specific heat of the vapor at constant volume, 
the specific heat of the condensate, R the gas constant, U the 
heat of vaporization and v the specific volume of the vapor-solid 
fluid. The gas volume is related to the fluid volume by v = v (1-x). 
Using this expression, the Claussius-Claperon relation (27) takes the 
form 


v/(l-x) = (35) 

Equations (34) and (35) are solved subject to the boundary conditions 

V = v^, T = when x = 0. 

When the state of the gas returns to saturation after passing through 
critical super-saturation, it is assumed that the T-v curve returns 
to the saturation isentrope crossed before going to the super-saturated 
state. Thus, the boundary conditions, V .and T are the volume and 
temperature when the gas becomes saturated but did not condense at 
time t.| . 

The degree of condensation is most easily written as a function 
of fluid temperature. 
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■ TEMPERATURE (K) 



F1g. 8. Temporal variations in specific volume and temperature of an expanding H 2 O vapor cloud 

with an initial pressure of one atm and an initial temperature of 2Q00i(. Temperature and 
specific volume of the bacground atmosphere are lOOOK and 3x1 m^/kg, respectively. 




where 


x(T) [(T^-T) (a^b^ +tT^} + (ai+b^) log^(T/T^)] 


(b 


(36) 


R + c, c, -• {c, + R) 


Note that as temperature approaches zero, x approaches unity, (A 
derivation of this important expression' appears in Appendix C.) The 
temperature is related to the volume by (35) and the volume is related 
to the expansion time through (24). 

The expansion-condensation process is illustrated in Figure 8 
for water vapor. At time t = 0 steam is released and the gas 
expands and cools adiabatically. After time t-^, the vapor becomes 
super-saturated. At the critical point, t = t 2 > the latent heat 
of vaporization is released during the rapid condensation to ice 
crystals and T jumps to the saturation temperature. The vapor-ice 
system then expands and cools following the saturation isentrope. 

The release of heat and the rise in vapor temperature during the 
condensation process is easily explained from a kinetic point of 
view. The free gas molecules that are attached to the nucleated 
clusters tend to have lower energy than the molecules that are 
reflected or ejected by the clusters. Thus, through nucleation, the 
lower energy gas molecules are sorted out and the temperature of the 
remaining vapor molecules is raised. 

The gas expands along the saturated vapor isentrope given by 
(36) until the gas is far removed from thermodynamic equilibrium. 
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In equilibrium, the temperature of the fluid is close to the 
saturation temperature. As the fluid becomes rarified, the rate 
of adsorption of gas molecules onto the clusters falls below the 
evaporation rate of the clusters and a lower cluster temperature is 
needed to maintain the cluster size. Thus, when the temperature 
necessary to support a condensation center falls significantly 
below the saturated vapor temperature, condensation ceases. At 
this time, t^, the maximum amount of injected material will have 
condensed. 

The flux of molecules towards the surface of a nucleated cluster 
is given by 


P-P 

^ [Hirth and Pound, 1963] 
y2TT mkT^ 

(37) 

where <p is the flux due to adsorption and evaporation, P^ the 
equilibrium saturated vapor pressure, P the equivalent flux vapor 
pressure at surface temperature T^ given by P = (T^/T^)^^^py, 

Ty the vapor temperature, P^ the vapor pressure, and m the 
molecular mass. This formula is used to determine the time of 
maximum nucleation and the evaporation rate of the cluster after 
nucleation ceases. 

The equilibrium growth rate is calculated as follows: Under 

equilibrium conditions, T^ and P are very nearly equal to T^ 
and Py, respectively. The equilibrium vapor pressure is derived 
from (27) 
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(38) 


P 


e 


R 

B 


-U/RT^ 


The effect of the curved cluster surface is neglected. The vapor 
pressure is given by the perfect gas law equi- 

librium the vapor pressure has the form of (38) (i.e,, Py=R/B 
exp[-U/RTy]). The flux of molecules into a cluster becomes 



2kT 


where v = 


V 0 “ 

'' the mean thermal velocity, and , T 


^s-^v 


the 


m • s 

super cooling factor. When 6 is substantially greater than zero, 

the equilibrium condition is no longer valid and condensation ceases. 

The growth rate is the flux multiplied by the surface area of the 

cluster. 


Hn ? r 

= 47rr 4) = — — — n _ exp . (40) 

1 /3 

where r„ is the cluster radius = r g , g’ the number of molecules 
per cluster, r^^ the molecule radius, and T^ = U/R. A similar 
expression is given by Zeldovich and Raizer [1967]... 

Noting that n^ = n 2 (l-x). (t 2 /t)^ and -g = x/v, Eq. (40) is 
written as 
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( 41 ) 


^ vl/3 x2/3n.x) ( yt)2 (l-e''''‘''^''=) ' 

where w , the collision rate for free gas niolecule's at the time 
m 

2 - 

(t 2 ) of critical supersaturation = 47ir^vn2 and n 2 is the gas 

density at the same time. The cluster number, v is found from the 

ratio of the cluster number density, N, to the fluid number density, 

-9 

n 2 - In the water vapor example, v = 1.4 x 10 clusters/gas 
molecule. 

The solution of (41) for 0 results in 


0 



l09e 1 


L 


t^ dx/dt 
x^'^^(l-x)A. 


(42) 


where 


A = 



3 




0 is much less than unity except when 


t^dx/dt 

x^^^{l-x)A 


> 1 


(43) 


This is the criterion for determining the time when condensation 
ceases. The values of x and dx/dt are derived from (25), (35), 
and (36). Figure 9 illustrates the temporal variations of x for 
the H 2 O vapor release. Using the data from these figures and (43), 
it is found that condensation stops at t^ = .11 sec. The temperature 
and degree of condensation at this time are T^ = 195 K and x^ = .24. 

g 

The number of molecules per cluster is given by g^ = x^/v = 1.71 x 10 
indicating a cluster diameter of .2 microns. 
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Three forces affect the motion of the nucleated clusters. 

Gravity exerts a downward force, displacement of the atmosphere 

exerts an upward force according to Archimedes' principle, and 

aerodynamic drag exerts an additional upward force. The force 

equation is ? = - + gm^ + where g is the acceleration- of 

gravity, m the mass of the cluster = Airr^ p /3, m the mass of 
c • c c a 

3 

the displaced air = 47rr P /3, r the cluster radius, p the 

C a O . Q 

cluster mass density, p^ the atmospheric density, and the drag 
force. Since the cluster is much more dense than the atmosphere, 
the buoyancy force is negligible. Assuming a spherical 
particle, the pressure drag is B = STiyr v [Hoerner, 1965] where y is 

I V 

the viscosity, and v the upward particle velocity with respect to the 

atmosphere. The acceleration of the cluster is dv/dt = - B + B /m = ■ 

r c 

- g - 9yv/(2r^p^}. The viscosity of the atmospheric constituents 

0, 0^ and N2 is approximately given by y = 4. x lO"^ T*^^ kgm"^(s‘'^) 

[Banks and Kockarts, 1973]. The equation of motion for the spherical 

particle becomes dv^/dt = - g - 1 .8 x 10”^ T’®^ v/(r^p ). 

c c 

For the example of a .2 micron diameter ice cluster in a 1000 K 

F-region atmosphere, the equation of motion becomes dv/dt = 

7-> 

- 9 - 2 X 10 V. 'The particle approaches steady state velocity with 
time constant x =.5x10^ sec. The cluster velocity approaches a 
terminal velocity of v = - 4.5 x 10'^ m/sec. Because of this low 
velocity, the falling motion of the nucleated clusters may be 
neglected. 

After thermodynamic equilibrium is disturbed, the condensed 
particle clusters freely evaporate. The evaporation flux is given by 
(37) with the external vapor pressure term dropped 
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(44) 


-P. 


= 


v^^irTT 

m 


The cluster surface vapor pressurej including the effects of the 
spherical radius of curvature, is derived from (28) 


P 

e 


= P 

CO 



(45) 


The flat surface vapor pressure, P^, is given by (38) resulting in 


P 

e 





(46) 


The evaporation rate is the evaporation flux multiplied by the cluster 
surface area.. 


^ = - f exp 






(47) 


The variation in cluster temperature is a function of cooling by 
evaporation and of heating by the sun's rays. .Heat is removed from 
the clusters by evaporation because the high energy molecules escape 
leaving the cooler, low energy molecules attached to the cluster. 

The temperature change is given by the heat equation 


dQ 

gm c 


2 dt 


= gmU + Q^kVe 


(48) 
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Fig. 10. Ice absorbtion coefficient (k) in cm~^ versus wavelength in microns. 






n of ice clusters. 




where Q is the heat energy of the cluster, C 2 the specific heat 

of the cluster, U the change in energy of the cluster associated 

with evaporation, the solar constant = 1340 watts/m^, R the 

average solar radiation absorption coefficient, V the cluster 
3 

volume = 4/37rr^g and e the solar heating factor depending on the 
sun's position (e = 0 at night and e =1 at noontime). The 
average absorption coefficient is calculated as follows 

00 

- Jo R(x)kU)dx 

k = — (49) 

V“R{X)dX 

•^0 

where k is the wavelength dependent absorption coefficient and R(x) 
is the solar power radiated per unit wavelength. The spectrum of the 
absorption coefficient for H 2 O is shown in Figure 10 as given by 
Irvine and Pollack [1968], Averaging over the solar spectrum yields 
kH^0=7770.-’ • 

Substitution of (48) into (49) gives 



where p is the mass density of the cluster. The cooling and evapo- 
ration of the nucleated clusters are found by solving the coupled 
first order nonlinear differential equations (48) and (50) subject to 
the initial conditions: 
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(51) 


T = Tg, g = Xj/v at. time tg 

The evaporation of nucleated ice clusters is illustrated in 
Figure 11. The difference between the nighttime and daytime cluster 
persistence is striking. After 1000 seconds, only 11% of the night- 
time ice cluster has evaporated whereas, during the day, 98% of the 
same cluster has evaporated. From these calculations, it is concluded 
that nucleated gas molecules have negligible effect on the ionosphere 
during the night, but they will enter into recombination reactions 
during the daytime after some delay. 


0. Transition Between Translational Equilibrium and Frozen Flow 

After condensation equilibrium is disturbed and condensation 
ceases, the expanding gas temperature is no longer influenced by 
cluster cooling. The vapor is in translational equilibrium and expands 
according to the continuum isentropic expression 


n 

"3 





3{y-1) 


(52) 


where n^ and T^ are the density and temperature of the gas at the 
time of deviation from saturated vapor equilibrium. The rapid cool- 
ing of H2O vapor after this time is shown in Figure 8. 

Under translational quasi -equilibrium conditions, the gas 
collision or relaxation time = l/w, is small in comparison with 
the time scale, x^, of the gas expansion. When Xj^ » x^ the flow 
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becomes frozen in a nonequilibrium state. The transition time, 

t^, between equilibrium and frozen, free molecular flow is calculated 

by equating 'the relaxation time with the characteristic expansion time. 

The relaxation time is the ratio of the internal mean free path 
to the mean thermal velocity. ^ 


" c 


m _] 1_ 

^ 2TTd^ n/r 


where X = 7=^ is the mean free path inside the cloud, d is the 
^ v2Trd^n 

molecular diameter, n the gas number density, k Boltzmann's 
constant, T the gas temperature, and m the molecular mass. 

If the released gas does not go through a state of condensation, 
expressions for T and n from (22) and (23) are used in the 
expression for the relaxation time to yield 


.*0 (‘^ 0 ^ 


R V, 




where t = 


If the gas does condense, the expansion after the terminal nucleation 
time, t^, can be found by using (52) for temperature and density 
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The characteristic expansion time is taken to be the' expansion 


velocity time constant.. 


1 _ dv 

Tg ■ Vdt 


( 56 ) 


Since v = r/t in both continuum and collisionless expansions into 
a vacuum, = t. 

Equating and and solving for the transition time t^ 

gives 

2 

t^ = t^ (v^ird^ (without condensation) 

( 57 ) . 

2 

t^ = t^ (^/^^^d^t 2 v^kT^/m n^) (with condensation) 

( 58 ) 


Another, equally valid definition of the time separating the 
continuum flow from the collisionless flow is obtained by letting the 
Knudson number be unity. 


K 


n 



= 1 


( 59 ) 


where is the mean free path inside the cloud and d^ is the 
cloud diameter. The cloud diameter is taken to be an equivalent 
diameter of a sphere with the cfentral cloud density 




( 60 ) 
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Equating d with X and solving for n yields ■ 

( 61 ) 


(62) 

Except for a factor of 2^^^, this is the same result as (57) with 
Y= 5/3. Thus, an estimate of the transition time between continuum 
and frozen flow for noncondensing gases is given by 

‘4 = ■ (63) 

For water vapor, this has the numerical value of t^ = .44 second.' 

The extiraate transition time including condensation, (58), yields 
a value of t^ = .67 seconds for the water vapor release. Since the 
transition time calculation is only approximate, values from (63) are 
sufficiently accurate in most cases. 

E. Interaction with the Background Atmosphere 

As the gas cloud expands, atmospheric particles will be imbedded 
in the outer shell of the cloud. The heating of the gas cloud by 
picking up atmospheric particles has been examined by Hoffman and 
Best [1974] and Baxter and Linson [1976]. This work only considers 


/t 

By using (23) n = 

t = t (2/STTdVn^)^''^ 
0 0 0 
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the atmosphere-to-cloud energy transfer when the rarefied atmosphere 
is well mixed with the injected cloud. At this time, the gas molecules 
in the cloud are in a state of free molecular flow with themselves. 

In this section, the gas expansion between the times of initial 
free molecular flow, t^, and diffusive flow, tg, is examined. 

The initial gas distribution at t^ is modeled as a uniform sphere 
with density .and temperature T^. The radius of the sphere is 
defined so that the total number of particles is conserved (i.e., 
N^(l-x) = where r^ is the equivalent cloud radius at time 

t^, is the number of molecules released, x is the maximum 
degree of condensation, and n^ is the central number density). 

With this assumption, the free molecular flow expressions given by 
(7), (8), and (9) are used. The time and distance are normalized 
in the usual way: t = tc^/r^ and r = |x}/r^ where c^ = /ZlTT^/nT. 

Since only the evolution of densities and temperatures near the center 
of the cloud will be considered, little error is made by using this 
model . 

Assuming an unperturbed background atmosphere and negligible 
self scattering, the released gas distribution satisfies the follow- 
ing kinetic equation. 

11+ V *^f = - f) + 6(t) Q^(x) (j)^(v) (64) 

where f(t,x) is the gas distribution function, to = c /A , the 

p 

ambient collision frequency, A^ = l/(v^ 7 rd n), the mean free path, d 
the gas molecules diameter, c^ = v/2'kT~7"C the mean thermal speed, n 

a a a 



the released gas density, the spherical source function, 

^a ” |“C(v - atmospheric velocity 

relative to the gas release, and <[>^ the Maxwell source velocity 
distribution of (1). This is the Krook equation with an atmospheric 
collision term and a spherical source term. Equation (64) is solved 
following the procedure used by Baum [1973]. 

The time and space derivatives are removed by operating with the 
Laplace-Fourier transform as defined in (65) and (66). 

00 CO 

g(p,k) = fd^x / dt g(t,x) (65) 

-OO Q 

g(t,J) = g(p,k) (66) 

-«> (2T7)^ c 

The transformed Eq. (64) becomes 

(p + ik-v + w^)f = Wg(!)gn = (j>QQQ(k) ■ (67) 

the overbar indicates transformed quantities. The solution for 
is written as 


f = 


p+ik*v+(o 


( 68 ) 


j3:»- 


Using n - Jf d v, the transformed density function is found to be 


. j3~^ 

- ^ 4>^d v 

n = Q„(l<) J — ° ^ 

p+ik*v+w 


1 - 


(j) d^^ 
“a / ^ 


- 1-1 


p+il<-v+(i) J 
a 


(69) 
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This density function is split into two parts 


n = 


n + n 
a c 


where 


and 


Sa ■ 


* d^v 


p+il<*v+aj. 


(70) 


fi,. = n G(p»k) 
c a 


(71) 


G(p,lc) = w. 




® p+il^-v+cj 




r 

p+ik*v+(jj, 


l-l 


aj 


(72) 


The first term, n , is related to the collisionless solution to 
the Boltzmann equation. The collisionless Boltzmann equation with a 
source term is 


9f. 


3t 


+ v.Vf|^ = 6(t) Q^(x) 4 .q(v) 


(73) 


The Laplace-Fourier transform of this equation is manipulated to yield 




f = 0 0 


m , 

p+ik-v 


(74) 
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Integration over velocity space yields 


r.= /fd3? = 


m 


%it) ! 


5+if*V 


Comparison with (70) shows that 


(75) 


"a<P> ° * “a> 

Finally we invoke the inverse Laplace-Fourier transform, Eq. (26), 
to get 


n 


a 


' "n, 


The first term in the density function is the free molecular flow 
solution multiplied by exp( - w^t) which represents damping by the 
atmosphere. The free molecular solution for spherical expansion, 
has already been derived and is given in (7). Hence, 



L 


(76) 


The second term, n » in the density function represents the 
collision effects of the atmosphere. After many collisions, n 
vanishes and the scattered gas molecules, represented by n^, 
constitute the expanding cloud. The multiplication of the two 
functions, n^ and 6(p,k), in (71) transforms into convolution 
resulting in 
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( 77 ) 


t o 

«€ = njt,x) * G(t,x) = '/ dt^ / d'^x-| njt^x^) G(t-t^, x-x^) 

where G(t,x), the scattering Green's function is the inverse 
Laplace-Fourier transform of G(p,lc). 

The exact form of G(t,x) is mathematically complex and the 
exact convolution in (77) is a formidable task. For large-times, 
Baum [1973] has derived the following asymptotic form for G(t,x). 


G(t*,r*) 


g-(r*)^/2t* 

(2TTt*)^^^ 


(78) 


where r* = |x*-W,t*j, x* = , t* = tco,, and W, = u,/c,. 

a A a a a a 

a 

(In the remaining part of this section, time, distance and velocity 
are nondimensionalized by l/u, , and c, — the characteristic 

atmospheric quantities.) The use of this aymptotic Green's function 
simplifies the calculations and does not introduce appreciable 
qualitative error because n 2 is significant only after long times. 
The long-time expression for the free molecular term, n^, is 


g(r*)2/(bt*)^ 
"m^Xa 3/2(,^*j3 


(79) 


where b = c^/Cg and = N^(l-x) is the number of gas molecules 
uncondensed. This is the free molecular density solution for an 
impulsive source release. The expression of the gas density of all 
finite, free molecular gas releases approaches (79) in the long-time 
limit. 
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Convolve exp(-u^t) with G(t*,r*) to get 



/ e"^‘ G(T,R)dt‘ 

0 


( 80 ) 


where R = |x*4^(t*-t' ) 1 , T = t*-t + 1/2(bt')^, and G(RJ) is 
given by (78) 

The composite solution is found by adding the damped free molecular 
solution, Eq. (76), to the expression for the scattered gas density, 

Eq. (80) 


n(t*,x*) 




) 


f ' e~^'G(T,R)dt‘ 
0 

(81) 


where a = r-/X b = c,/c . . x* = x/X , t* = tw^, r* = |x*i 

o a . ^ s . a ■ _r 2/2T 

R = lM,(t*-f)|, T = t*-f + l/2(bt')^, and G(T,R) = ^ ^ 

The composite density solution consists of two modes, the 
collisionless mode and the scattered mode. This bimodal nature is 
evident if the gas is released at a supersonic velocity with respect 
to the atmosphere (W^>1). The center of the collisionless mode 

•• a 

continues in the inertial reference frame of the released gas. The 

V 

center of the scattered mode will remain fixed with respect to the 
atmospheric point of release. Figure 12 illustrates the conversion 
from the collisionless mode to the scattered mode for a relative 
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Fig. 12. Interaction of a moving gas ball with a stationary 
atmosphere. An expanding hydrogen cloud containing 
No molecules is released at Mach 2 into an atmos- 
phere with meanfree path Xa and mean collision 
frequency Most of the injected gas becomes 
thermal i zed with the atmosphere after 10 collision 
times. 
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release velocity of Mach 2. To construct the Isopycnic contours in 
this figure, the integral in (80) is evaluated by the Gaussian 
quadrature method. For times less than one mean collision time, the 
spherical gas cloud is only slightly perturbed by the atmosphere 
(Figure 12(a)). At about 3 collision times (Figure 12(b)), approxi- 
mately half of the gas is in the free molecular state and the other 
half is moving with the atmosphere. The scattered gas molecules, are 
hotter than the molecules that have not collided with the high 
temperature atmospheric gas. The high thermal velocities of the 
scattered mode result in a higher expansion rate for the scattered 
mode compared to the collision! ess mode. This effect is represented 
by the widening of the density contours in going from the collision- 
less to the scattered regions. After many collision times, the 
injected gas flow approaches a spherically symmetric distribution 
about the point x* = where x* is the position of the point 

of gas release in the inertial frame of reference (Figure 12(c),). 
Note the change of scale in Figure 12. 

When the gas expansion is controlled by atmospheric collisions 
(t* » 1), Eq. 80 becomes 



This expression can also be derived assuming diffusive expansion 
(Chapter IV). 

The temperature of the released cloud is dependent on the loss 
of internal energy to radial expansive motion and energy transfer 
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from the ambient atmosphere through collisions. The temperature is 
found by calculating the second- velocity moment of the distribution 
function in (68) which is repeated here 


f = 


P+lK*V+(jJ 


Every variable in this expression is known. Density n, the inverse 
Laplace-Fourier transform of n, is given by (81). 

Considering the collisionless Boltzmann equation 


3f_ 


m 




and the solution for f using the L-F transform 

m ^ 


(j) Q 

^ _ ^ 0^0 
"" ' p+il<-v 


Comparing this equation with (68) we see 


f = f_+f (p+w_) 
c m a 


where 


(83) 


(84) 


(85) 


f — a a 
C 


p+i'^-V+03 
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Thus, the particle velocity distribution function is the sum of two 
distribution functions as follows 


f = 


f + f 
c m 



t 


(86) 


The free molecular distribution function, f^, has been previously 
derived (6). 

The collisional contribution to the distribution function, f , 
is found from the following convolution 


f 

c 


(L-F)- 


p+ik*v+u. 


] a^a 


(87) 

where (L-F) ^ indicates the inverse Laplace-Fourier transform 
operation. It is easy to show that 


(L-F) 


-1 


p+it-v+ci). 


= H(t) 


6(x-vt) e ^ (88) 


a-J 


where H(t) = 0 t<0 

1 t>0 

The result of the convolution, expressed in nondimensionalized units, 
is 


f 


c 


1 


a 



/ n[T-| ,x*-v*(t*-T,)] e ' dx, 

0 1 I . I 

(89) 


where v* 


v/c^ and n[x^ ,x*-v*(t*-x^ )] 


is obtained from (80). 
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The temperature is found at the center of the cloud (x = 0) 
under conditions of no relative atmospheric motion (W = u^ = 0). 
The temperature variation is composed of a damped free molecular 
contribution and an atmospheric collision contribution. 


T = 


i ,[l ill i /// 

•>03 «>CO —CO 


-w t 


m 


+ T, 


(90) 


T|^ is given by (12). In nondimensional ized form 




erf 


bt* 






JL 

n/, 


(91) 


T is found by calculating the second moment of f -as given by (89), 
^ c 




T_ - 


- 3k 


/// V 


I /// fe 


(92) 


At X* = 0, the term inside the integral in (89) becomes 


n[x^ ,-v*(t*-T^)] exp[-(t*-T^ )] = ^3 


exp{-[v*(t*-T^ )/(bx^] -x^} 


(bt,)3 


x-j exp{-[r*(t*-x^ )] /T-X2> 
+ f dx« 


0 


(ttT) 


3/2 


(93) 


70 



2 

where T = 2 {t-j-T 2) + (bT 2 ) and the long time damped free molecular 
density (79) has been used. Substitution of (93) into (89) and 
integrating over velocity space results in 


^ ^ / dxi « - 


^a nX ^ '0 
a 


IT 


3/2 


""1 (br^)^ 




e-*T 


(94) 


The central density variation, n(x* = 0), takes the long-time form 


n(o,t*) 


A 

>3 


e-t* 






-T 


1 dx 


1 


IT 


3/2[2(t*-x,)+(bt)2]^/^ 


(95) 


Combining (94)' and (95), we get an expression for the temperature 
increase produced by collisions with the ambient atmosphere. 


T_ = 


t* -(t*-T,) 

"4 


-Tn 


(A^+(t*-x^)^/^ •'0 [A2+(t*-T^)]^""J 


-I 




-T2 

d„e A. 


5/2 


-t* 


.t* 




I “^1 


-T1 


(96) 




' where 

A-| = (bx-j)^ b-j = (bx*)^ 

B 2 = 2(x*-t-|) + (bx-j)^ 

The density and temperature expressions derived in this section are 
used to formulate the T-V variations of Figure 8 during the period when 
the released gas is thermal izing with the atmosphere. 
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Fig. 13. Temporal variations in specific volume and temperature of an expanding H 2 cloud 
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Fig. 14. H 2 saturated vapor curve. 



F. Discussion 


The expansion of molecular hydrogen differs from the expansion of 
water vapor for several reasons. Since the H2' liioTecule is diatomic, 
its energy is partitioned into fewer degrees of freedom than the 
tri atomic H^O molecule. Thus, a greater portion of energy goes into 
the translational motion of the H2 molecule as indicated by the 
higher specific heat ratio for H2. Also for a given temperature, 

H2 molecules have higher average velocities than the H2O molecules 
because of the molecular weight difference. The larger values of 
specific heat ratio and mean thermal velocity for H2 results in a 
higher expansion and cooling rate over H2O. This is evident from a 
comparison of the time marks for the H2 {Figure 13 ) and H2O 
(Figure 8 ) expansions. 

Condensation of H2 is' improbable because of the extremely low 
values of saturated vapor temperature (Figure 14 ). Since no latent 
heat of vaporization is released, the expanding Hg gas temperature 
drops to a very low value before heating by atmospheric collisions 
occurs (Figure 13 ). Of the four most reactive gases (H2» H2O, CO2 
and NO2). molecular hydrogen is the only one which does not condense 
appreciably after release. 
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Chapter IV 

DIFFUSION IN A NON-UNIFORM, CHEMICALLY REACTIVE ATMOSPHERE 

Diffusion is the statistical process by which particles inter- 
mingle as a result of their- random motion. Diffusive motion of a gas 
is a function of the thermal velocity, and the collision frequency of 
the gas molecules. The diffusive migration of gas molecules results 
from gradients of concentration or temperature. 

A. Diffusion in a uniform Atmosphere 

Diffusion in a uniform, one-dimensional medium can be modeled as a 
random walk. Consider a particle released at the origin of a one- 
dimensional coordinate axis at time zero (Figure 1). Imposed upon the 

PARTICLE AT x«0 WHEN t=0 

1/2 1/2 

° * I 

■~H \ ! ! 1 1 1 \ 1 ► 

-4Ax-3Ax-2Ax -Ax 0 Ax*2Ax 3Ax 4Ax x 

Fig. 1. One-dimensional random walk model of diffusion. 

particle is probability one-half of decreasing its position by Ax in 
time. At and probability one-half of increasing its position by Ax in 
time At. This is a model of a medium with mean free path Ax and mean 
collision time t. If U(x,t) is the probability that the particle is 
at position x at time t, the following difference equation results 
from the law governing the motion of the particle: 

U(x,t+ t) = 1/2 U(x- x,t) + 1/2 U(x+ x,t) (1) 

The initial condition is U(x=0, t=0) = 1 and U{x/0, t=0) = 0. The 
temporal change in probability at a fixed location is written with the 
aid of (1). 
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(2) 


U(x,t+At) - U(x,t) ^ ^ 
At 2 At 



Taking the limit of Ax*H) and' At-^ while keeping' the quantity 
constant gives 




(3) 


where the following definitions of partial derivatives are used 


lim 

At^f) 


U(x,t+At) - U(x,t) 


3t 


2U(x,t) + U(x-Ax,t) _ a^U 

2 ■ 2 

Ax ax 


If many particles are released, the single-particle probability 
function U(x,t) is equivalent to the macroscopic multi -particle 
concentration. Thus, (3) is the simple diffusion equation for flow 
in a one-dimensional uniform medium. D, the diffusion coefficient 
with dimensions of length squared per time, is the square of the 
distance between collisions divided by-twice the time between collisions'. 
In the atmosphere, the distance between collisions is the mean free path 
and the time between collisions is the reciprocal of the mean 
collision frequency w_. The atmospheric diffusion coefficient is 

a 

written as 
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The mean free path is related to the collision frequency by 


= c = /2kT /m 

ad a di B. 


( 5 ) 


where c^ is the mean thermal velocity, the atmospheric tempera- 
ture and m the molecular mass. The diffusion coefficient becomes 

O 


D = 



"^a“a 


( 6 ) 


In three dimensions, the diffusion equation (3) is 


9n _ 
3t 





(7) 


where n(x,y,z,t) is the partial number density. The release of a 
large number of particles at a point implies an impulsive initial 
condition. 


n(x,y,z,o) = 6(x,y,z) 


where is the total number of particles released. With this 
boundary condition, the solution to the 3-dimensional diffusion 
equation is 


(4irDt) 


3/2 


(x^+y^+z^ ) 
e 4Dt 


n(x,y,z,t) = 


( 8 ) 



Substituting the mean free path and collision frequency factors for D 
(4) and using the nondimensionalized time and distance of the previous 
chapter (i.e., t* = tw and x* = results in 

a, a 

N 5 

n{x*,t*)= ^ ° 0/2 exp[~x*V(2t*)] (9) 

a 

This is identical to the asymptotic expression for a uniform sphere 
expanding into an ambient atmosphere after the released gas has 
thermal ized with the atmosphere (Chapter Illj Equation (82)). 

B. General Diffusion Equation 

The flow of the injected gas has been calculated neglecting the 
effects of gravity, of a non-uniform atmosphere and of chemical 
reactions. These effects are included by- solving a general diffusion 
equation derived from the equations of continuity and momentum for 
the diffusing gas cloud. 

The equation of continuity for a neutral gas is 


V[n?] = p - L (10) 

where n is the relased gas density, v the velocity of the gas, 

P the production or source term arid L the loss term. The continuity 
equation is derived from the Boltzmann equation under the assumption 
of continuum flow. 

The velocity of the gas is given by the momentum equation 
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( 11 ) 



k 


V(nT) 

nT 


+ 0 ) 


a 



where g = Is the acceleration due to gravity, a^ the 

vertically upward unit vector, T the gas temperature, tL the 

. a 

neutral drift velocity, and w the collision frequency. For 

a 

velocity fluctuations with time constants greater than a mean collision 
time, accelerations can be neglected, dv/dt = 0, and velocity is 
found to be 



nT Hg 



( 12 ) 


where D is the diffusion coefficient, and H = kT/mg is the gas 
scale height. The velocity of the gas is composed of two parts: the 
diffusion velocity produced by concentration and thermal gradients 
and the drift velocity due to the motion of the bacground atmosphere. 

The general diffusion equation is found by substitution of (12) 
into (10) 





This reduces to the simple diffusion equation (7) if D and T are 
constant, and gravity, atmospheric motion, production and loss are 
neglected. Equation (13) will be solved analytically but first the 
properties of the, diffusion gas cloud are discussed. 

The expansion velocity of the diffusing cloud is one-half that 
of the cloud expanding into a vacuum. In a uniform stationary 
atmosphere, the expansion velocity given by (12) becomes 
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(14) 


v = -0f 

Substitution of the spherical diffusive expansion solution, Eq. ,{,8) ' 
yields 


V 




r 

2t 


(15) 


to be compared with the vacuum-like expansion velocity v = r/t, 
derived in the previous chapter. 

The motion of the center of mass of the gas cloud due to 
gravity and atmospheric density gradients is calculated as follows. 
The center of mass is defined as 


a(t) = jp ft n(x,t) d^^ (16) 

where a(t) is the time varying position -vector of the center of 
mass, is the total number of gas molecules, and x is the 
position vector to a. point in space. The velocity of the center of 
mass is 


dt N ^ ^ 


0 


at 


(17) 


Using (13) 


da^ ^ 
dt 




0 


dXi + ^ 

n 


d^x 


(18) 
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u 


Assuming vertical atmospheric density gradients of the form = n^^ 
exp{-z/Hg) and thus exponential variation is collision frequency with 
altitude (w = exp(-z/H,) where n^ is the atmospheric density, 
is the density at z = 0, = kT^/mg is the atmospheric scale 

height and m is the atmospheric collision frequency at z - 0), 
the diffusion coefficient is inversely proportional to the collision 
frequency leading to 


0 = exp[z/Hg] 


( 20 ) 


With this substitution, the expression for the velocity of the center 
of the mass becomes 


da _ 
dt 



g f~o 


n(x,t) ^3 

No 




/ D(z) ^W^M'x = 



( 21 ) 


where D{t) is the average value of the diffusion coefficient inside 
the cloud. The buoyancy of the cloud depends on the relative scale 
heights of the atmosphere and the released gas. If the average 
atmospheric gas molecule is heavier than the released gas molecule, 
will be less than and the cloud will rise. Similarly, if the 
cloud is composed of a gas with a larger molecular weight than the 
background it will fall., Thus, the familiar bulk motion of gases in 
the atmosphere are inherent in the solution of the diffusion equation. 
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C. Solution of the General Diffusion Equation, for Point Gas Releases 


After the injected cloud has thermal i zed with the ambient atmos- 
phere, its expansion is indistinguishable from that resulting from 
the solution of the diffusion equation assuming a point release. In 
the early stages of the diffusive expansion, gravity and the atmos- 
pheric density gradients have little effect on the expansion and the 
solution has the spherically symmetric form of (9). When the gas 
expands over many atmospheric scale heights, the effects of gravity 
and atmospheric density gradients become noticeable. 

The diffusion equation for the injected gas is written as 








( 22 ) 


4- + 

where y is the (0 , injected gas) reaction rate, and [p ] the 

0^ density. A Cartesian coordinate system, centered at the point of 

gas release, with the z-axis pointing vertically upward is used. 

Assuming an exponential atmosphere yields D(z) = DQexp(z/H^) 

where D is the diffusion coefficient at the point of gas release 
0 

{z = 0), and H the atmospheric scale height. The last term in 
(22) represents the chemical depletion of the injected gas by reactions 
with the ionospheric 0^. This term is negligible compared with the 
diffusion terms which dominate the gas behavior until the cloud is 
dissipated to an insignificant density level. The chemical loss term 
has appreciable effects only during the last .stage of the gas expansion 
when the ionosphere is well on its way to recovery. Therefore, this 
small term is replaced with an approximation that greatly simplifies 
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the solution of (22). A constant a = y[ 0^] is .substituted for 

4* 4* “h 

y[0 ] where <[0 ]> is the average 0 density in the unperturbed 

ionosphere. 

After substitution for the diffusion coefficient and the chemi- 
cal loss term, the diffusion equation becomes 

z/Ha 

3n _ ^ 0 ® 3 f.dn\ . 3 

3t “ r 3r 3z 

2 2 2 

Cyclindrical coordinates are used where r = x + y . The initial 
condition is n(x,y,z,0) = . The analytic solution of this 

equation without the chemical loss term is given by Yu and Klein 
[1964]. The solution is outlined as follows: Eq. (23) is Fourier 

transformed with respect to the variables x and y yielding 




of 

2 

3r 

a 




3N^ 

_ 


+ 4D. 


1 


a g 





(24) 


where C = exp(-z/2 H_), f and f are the transform variables 

a X y 

corresponding to x and y, respectively, and is the Fourier 
transform of n. This equation is Laplace transformed with respect 
to the variable t. 



+ 4 


^ - H?(f? + ff,)(2n)^ - 


D_ 


(p+a) 


-2r 6(C-1 ) 

a 0 


(25) 
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where p is the transform variable corresponding to t, and 
the Laplace transform of N^. The homogeneous part of (25) is of 
the form 


^ + 39 II + (b+cg^) f = 0 
9g ^ 


which has the Bessel function solution 


• f(g) = 9 ^ 


[Kamke, 1959] 


where v = 1/2 and is any linear combination of 

Bessel functions'^ of the first kind [J(*)] and the second kind 
[Y(0]. The solution of (25) with the boundary conditions as 

or oo is 




/H N \ 

1+— 

1 r 

1 a O' 
\ 

|5 ^ iri 1 


J (2i»paC) H^n( 2 i,/p+ 5 ) e(l-5) 


where v = 
function of 


J^(2ii/^ H^^^(2iv^ 5) 6(5-1) 

(26). 


3 _ i) +4 H^(f^ + f^) (2 tt)^ 
a A y 


is the Hankel 


the first kind and e(x) is zero for negative x and 


unity for positive x. The inverse Laplace transform of (26) is 




tD^ ^ 


-(?^+l)Hg 

tT 


- at 




(27) 
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where a = modified Bessel function of the 

first kind. Applying the inverse Fourier transform one obtains 

J^(2Trrp)dp 
(28) 

where = (a-1)^ + (4irH p)^ 

a 

The integral in (28) can be evaluated numerically, but during 
the time interval of interest, the Bessel function inside the integral 

q 2/ 

can be closely approximated by I (q) == -- — e~^ ' Justification 

^ / 2 ^ 

for this approximation is based on examination of the following series 
for the modified Bessel function 
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This expression gives values of n within a factor of two of those 
given by (28) for times less than half an hour under typical F2-region 
conditions. In the limit as time approaches zero, (31) -reduces to 


n(x,y, z,t) 


( xW+z^) 

■ 


(4TTD„t)^/^ 


which is the spherically symmetric expression for expansion into a 
uniform atmosphere. 

The effects of gravity and the non-uniform atmosphere are seen 

in the time sequence of gas expansion in Figure 2. The time unit is 
2 

defined as H,/D^ at the point of gas injection. The isopycmc 

a 0 

contours containing the center of the cloud are plotted at .01, .1, 
and 1 unit times. The contours are nearly circular for cloud diameters 
less than one scale height. When the gas expands over many scale 
heights, the effects of atmospheric density variations becomes apparent. 
The increase of the diffusion coefficient with height produces rapid 
transport of gas molecules into high altitudes. The flow at low 
altitudes is impeded by the dense atmosphere and a hard bottom forms 
on the cloud. The location of the peak gas density falls with time 
because of the rapid diffusion out of the top of the cloud. 

The atmospheric scale height and injected gas diffusion coef- 
ficient are determined from atmospheric and injected gas parameters. 

The F-region atmosphere contains, in order of abundance, monatomic 
oxygen, diatomic nitrogen, and diatomic oxygen. Because of the 
variations in the relative concentrations of these neutral species. 
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ALTITUDE/H 



-3- -2-10123 


Fig. 2. £xpansion of a gas cloud In a nonuniform atmosphere. 

The cloud becomes highly nonspherica.l after expanding 
over many scale heights. An expanding Hp cloud con- 
taining molecules is released into an atmosphere 
with scale height H . The diffusion coefficient at 

the release point is'^D . 

0 
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the scale height varies with altitude. In the injected gas expansion 
model, the scale height at the point of injection is used. The 
scale height is given by 

H = .848(1 + h /R T/M- (32) 

a 0 e 

where h is the altitude to the point of injection, the earth's 
0 - •= - 

radius, T the neutral temperature, and M the average molecular 
weight. At 300 km altitude, the scale height is approximately 50 km. 

The diffusion coefficient at the gas release point is dependent 
on the local temperature, and the composition of the injected and 
background gases. A reasonable expression for binary gas mixture 
diffusion coefficients is of the form 

D^2 = AT^n (33) 

where is the mutual, diffusion coefficient for the binary 

mixture of gases 1 and 2, A a function only of the diffusing species, 
T the neutral temperature, s a species dependent exonent, and n 
the sum of the densities of the two diffusing gases. When the 
injected gas. is in a state of diffusion, its density is less than 
the ambient atmospheric density. The density used in the diffusion 
coefficient is taken to be that of the background gas. 

The constants A and s, determined from experimental data (when 
available) or from kinetic theory, are given in Table 1 for selected 
reactive injectants. 
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TABLE 1 


Diffusion Coefficients D = AlVn 


System 

A 

s 

A ^ 

s 

0 

A- 

s 


— 

— 

8.29[16] 

.724 

9.69[16] 

.774 


8.29[16] 

.724 

— 

— 

9.69[16] 

.774 

N 2 O 

2.34[17] 

.5* 

2.37D7] 

.5* 

5.87[17] 

.5* 

CO2 

6.58[16] 

.752 

5.77[16] 

.749 

5.87[17] 

.5* 

NO 2 

2.31 [17] 

.5* 

2.35[17] 

.5* 

8.84[17] • 

.5* 


2.8 [17] 

.740 

3.06[17] 

.732 

2.97[18] 

.5* 

H 2 O 

2.04[17] 

.632 

2.02[17] 

.632 

8.46[17] 

.5* 

*Theoretical Values 

Note: [x] means 10^ 




The measured values of the diffusion constants are from Mason and 
Marrero (1970). 

The theoretical values of the diffusion coefficient are based 
on the model of rigid, elastic spherical molecules. This model leads 
to the parameters 

A = .75/{d^+d2)^ kN (m^+m2)/(7nn^m2) and s = .5 

[Chapman and Cowling, 1970]. 

where d-j and d 2 are the molecular diameters of the binary dif- 
fusion system, m^ and m 2 the mass of the- molecules, and k 
Boltzmann's constant. The theoretical temperature dependence exponent 
s = .5 differs from the typical measured exponent s = .7 because of 
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the deviation from the assumed elastic-sphere molecular encounters. 
The diffusion coefficient for a gas in a multi -species mixture is 
given as 




(34) 


where Djj^ is the binary diffusion coefficient for the j 


,th 


and k 


th 


species. Typical values of diffusion coefficients for and H 2 O 

11 2 

in the 300 km atmosphere are D(H 2 ) = 1,1 x 10 • cm /sec and D(H20) 
= 3.4 X 10^^ cm^/sec and D(H20) = 3.4 x 10^^ cm^/sec. 


D. Solution of the Diffusion Equation for a General Gas Source Term 

The solution of the linear diffusion equation (23) for any initial 
gas release may be expressed in terms of the solution for a point 
source. The linear diffusion equation with source term 0(x,y,z,t) 
is written as 


9t 




- an + 3(x,y,z,t) 


(35) 


the source term for an instantaneous point release, 3(x,y,z,t) = 
<5(x,y ,z)6(t) , represents the gas injection at point (x,y,z) = 
(0,0,0) at time t = 0. The solution with this impulsive term is 
called the Green’s function for (35). The solution for a general 
source term is written as 
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t «> 

^ ill! 3{a,3,A,Y) G(x-a, y-3, z-X, t-y) dadedAdy (36) 

0 - 0 ) ■ • 

1 

where G(x,y,z,t) is the Green's function. 

As an example, consider a moving source, such as a rocket, 
continuously releasing gas at a constant rate over m a period of 
time lasting t seconds. Such a source may be modeled as an 
infinite sequence of uniform spherical gas releases. After the 
injected gas has thermal ized with the atmosphere, each release has 
the form of the impulse response or Green's function solution to the 
diffusion equation. This model for the source term is expressed as 

N 

B(x,y,z,t) = G(x-v^t, y-v t, z-v t, t ) for 0<t<t (37) 

‘'p A y z 0 p 

where v^^,v^,v^ are the x,y, and z components for the time dependent 

velocity of the source, t . the time duration of the source, N /t 

y o' y 

the number of particles released per unit time and t^ the time 

delay in the Green's function to account for thermal izati on with the 

atmosphere and the finite width of the emitting source. If t is 

0 

set to zero, the source function is a moving impulse. 

N 

6(x,y,z,-t) = 6(v-v t, y-v t, z-v,t) for 0<t<t,, (38) 

A J Z 1-1 

If the released gas distribution is identical to the Green's function 
response to an impulsive release some t^ seconds before, the 
release will expand as if it was created by a point release at time t^. 
Mathematically this is expressed in terms of the following theorem: 
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If 3{x,y,z,t) = G{x,y,z,t^) 6{t), then 

t “ 

n(x,y,z,t) = / /// 6(a,3,A,t ) 6{y) G(x-a, y-3, z-X, t-y) dadpdAdy 

0 -oo ^ 

00 

- /// G(a,g,X,t ) G(x-a, y-B> z-A, t) dad3dAdy 

-00 

= G(x,y,z, t+t^) 

This theorem is used in the derivation of the gas expansion for a 
moving source. 

The moving source model, (37), is substituted into (36) yielding 
the expression 

N min(t,t ) «> 

n(x,y, 2 ,t) = # / /// G{a-v y, B-v y, z-v y,t ) 

V -^0 X y z 0 

(40) 

G(x-a, y-B, z-A, t-y) dadBdXdy 


where m1n(t,t ) is the minimum value of time t and t . By 

Eq. (39), 

N min(t,t ) 

n(x,y,z,t) = X- / ^ G(x-v y, y-v y, z-v y, t +t-y) dy (41) 

-V o X y z 0 

The expression for the diffusive expansion of a point release in 
a non-uniform, reactive atmosphere, (31), may be used as Green's 
function but for simplification of the mathematics consider diffusion 
in a uniform atmosphere. Green's function for diffusion in a uniform 
non-reactive atmosphere is 
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Fig. 3, Density contours for a source of gas moving at 10 km/sec 
continuously releasing gas at a rate of 10 kg/sec. 
Release starts at zero seconds and shuts off at ten 
seconds. 





(42) 


■G(x,y,z,t) = exp[-(x^V+z^)/(4Dt)]/(4TTDt)^/^ 


The solution for diffusive expansion of a moving source into a uni- 
form atmosphere becomes 


N min(t,t ) -^^‘2 3/2 

n(x,y, 2 ,t)= ir-y ^ exp[-(x-VY)^/4D(t -t-Y)]/]4TTD(t +t-Y)]'^'^ ay 

Vo 0 0 


(43) 


The diffusive expansion of a rocket plume in a uniform atmosphere is 
illustrated in Figure 3. The rocket travels horizontally, from left 

4 

to right with a constant speed of 10 m/sec releasing H 2 at a rate 

of 10/kg/sec. The gas release lasts 10 seconds. The isopycnic 

/ 

contours along the rocket trajectory are shown at 5, 10 and 15 
seconds after initial gas release. While the rocket is releasing 
gas, a high concentration resides near the rocket nozzle (Figure 3(a) 
and (b)). After the release is terminated, the concentration is more 
uniformly distributed along the trajectory (Figure 3(c)). In a 
uniform atmosphere, the exhaust gas distribution approaches a uniform 
sphere in the long-time limit. 


E. Diffusion of Charged Particles 

The flow of ions and electrons is affected by magnetic fields, 
electric fields, density gradients, and collisions. Charged particles 
are constrained to move along magnetic field lines if magnetic forces 
are stronger than collisional forces. In the F-region ionosphere, 
diffusive flow of plasma is so constrained. 
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The plasma velocity along geomagnetic field lines is 


V 


-D 


<tn(npTp) 

3s 


sin I 


+ V. 


(44) 


where Op is the ion or electron density, Tp is the plasma 
temperature = (T +T. )/2, H is the plasma scale height = 2T k/(mg), 

e 1 p p 

I is the magnetic field dip angle, s is the distance along the 
magnetic field line, D is the effective ambipolar diffusion coef- 
ficient = 2kTp/(mu) and Vp is an externally -imposed drift velocity. 

The continuity equation for ionization trapped on magnetic 
field lines is 


^ 3 ^ 

8t 3s 


P - L 


(45) 


The field line path is related to the vertical axis coordinate by 
■s = z/sin I. The diffusion equation“i^found to be 



3t 3Z 


sin I D n_ 


^a'JogeCHpTp) 
3z 



np sin I Vp 


+ P - L 


(46) 

where L represents plasma loss due to recombination or charge 
exchange (Chapter II). It should be emphasized that the coordinate 
2 is the vertical distance to a point on a given magnetic field line 
and the derivative with respect to z is taken along this line. 

The interaction between the injected gas and the ionospheric 
plasma is represented by the coupled differential equations (22) and 
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(46). The injected gas density is depleted by expansion and chemical 
react on with the F-region 0 . If chemical depletion is small, (31) 
is used. Equation (46) can either represent depletion of 0 by the 
injected gas or variation in the density of the minor ion produced 
by reactions with 0 and the injected gas. Before these equations 
can be solved, a description of the ion, electron, and neutral tempera- 
tures must be obtained. 

F. Discussion 

’ ■ 1^2 has a higher diffusive expansion rate than H 2 O because of the 
relative low mass and collisional cross section of the hydrogen 
molecule. The H 2 diffusion coefficient is 23 times the H^O 
diffusion coefficient for expansion of either gas in the same atmos- 
phere. 

The comparatively rapid expansion of hydrogen may lead to diffi- 
culties in producing a desired reduction in the F-region ionization. 
Because the lighter molecules have higher expansion velocities, the 
concentration of light gases rapidly decreases and less localized 
chemical reaction occurs. Consequently more molecules must be released 
for a light, rapidly expanding gas than for a heavy, slowly expanding 
gas to produce the same depletion at a fixed point in space. 

If a strong localized depletion of the ionosphere is desired, 

H 2 O injection may be preferable to H 2 injection^-provided the 
condensation effects discussed in Chapter III are tolerable. 
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Chapter V 

THERMAL BEHAVIOR OF THE DEPLETED IONOSPHERE 

The F-region thermal structure is characterized by the neutral, 

ion and electron temperatures. The relationship between these 

temperatures is determined by the heat source into and the collisional 

/ 

interaction between the neutral, ion and electron gas species. Each 
of these gases acts as a heat source or sink for the other gases. 

A. Thermal Structure of the Neutral Atmosphere 

Because the neutral gas in the F-region is many orders of magni- 
tude more dense than the plasma in this region, its thermal structure 
is independent of the ion or electron thermal structure. The main 
heat source for the neutral gas is solar ultraviolet radiation. The 
photoabsorption by neutrals may produce excited gas molecules which 
thermal ize by collisions with other neutrals. The absorption of 
ultraviolet may photoionize instead of excite the neutrals as dis- 
cussed in Chapter 1. Alternate sources of neutral heating are Joule 
heating by electric currents and tidal dissipation of energy. 

Cooling of the neutral atmosphere is either by conductive or 
radiative processes. Heat flows by conduction from the high tempera- 
ture rarified upper atmosphere to the cooler lower atmosphere. The 
local heat conductivity is approximately proportional to the temperature 
to the .7 power [Banks and Kockarts, 1973]. Emission of photons by 
excited neutral atoms results in radiative cooling of the atmosphere. 
Much of the radiation is emitted in the infrared spectrum. The 
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effectiveness of cooling by the release of optical energy is reduced 
by re-absorption of infrared by the atmosphere. 

Because of the complexities and uncertainties of the physical 
processes determining the atmospheric temperature profile, an empiri- 
cal model is used in this work. 


T = - [T^ - T(z^)] exp[-a(j)] [Antoniadis, 1976] (1) 


where T is the exospheric temperature, T(z ) is the temperature 
at a lower boundary at altitude z^^, a is a constant determining the 
shape of the curve and <J> = (R+z^){z-z^)/R+z) is the geopotential 
height given in terms of the Earth's radius R and altitude z. 
Solving the steady state diffusion equation with this temperature 
profile results in 


«i = n.(z^)exp(ac{)) 


T(z ) 
^ o' 



T^[exp(a4)) - IJ + T{z^) 


where y- = 1-43 m./aT , n. is the concentration and m. is the 
1 1 1 1 

molecular weight of the i^^ neutral constituent (i.e., ^ 2 ^ ^ 2 * 

The neutral temperature and concentration profiles remain essentially 
unchanged throughout the period of ionospheric chemical modification. 


B. Thermal Structure of the Electron Gas 

The electron temperature in the depleted region is significantly 
perturbed as a result of the release of reactive gases. The electron 
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gas heat source and cooling processes are significantly affected in 
the region of reduced plasma density. 

The heat balance equation for electrons is expressed as 




3 , 

2 n^kv 



+ n kT 
e 


e 9s 



( 2 ) 


where s is the coordinate -parallel to the magnetic field line, n^ 
the electron' concentration, the electron temperature, X the 
electron thermal conductivity, Q the heat energy input, and [ the 
energy loss due to heat transfer to other species. One-dtmensional 
heat flow along magnetic field lines may be assumed in regions where 
the ratio of collision frequency to gyro frequency is small. This 
condition is satisfied in the ionosphere above 150 km altitude. The 
energy transport terms due to convection ^|- nkv and compression 

9s^ often neglected in modeling of the unperturbed ionosphere 
because they are usually small in comparison with the heat production 
and loss terms. However,’ these terms must be considered in modeling 
of the depleted ionosphere because of the large fluxes of plasma 
into the ionospheric hole; 

The electron thermal conductivity is determined by particle 
collisions. In a fully ionized plasma, the conductivity is dependent 
on electron-ion collisions and is given by 

X = 7.7 X 10^ Tg^/^ ev cm”'' sec'"' K"'' (3) 

This expression provides an adequate description of the conductivity 
in the unperturbed ionosphere. In the weakly ionized plasma in the 
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depleted region, electron-neutral collisions become important and a 
more general expression must be considered. 


X = 


7.7 X 10^ 

1 + 3.22 X 10^/ — ) 2 
\"e/ k 



(4) 


where nj^ and Qj^ are the density and collision cross section of the 
neutral species. The electron thermal conductivity is substantially 
reduced in the ionospheric hole. 

During the day, energy is transferred to the electron gas by- 
collisions with photoelectrons.' Photoelectrons are produced by • 
photoionization of the neutral gas. Expressions for photoelectron 
heating are given by- 'Swartz and Nisbet [1972]. The amount of heat 
energy transferred from photoelectrons to the electron gas is less 
in the depleted region than in the unperturbed ionosphere. This 
results from a larger fraction of the photelectron energy being ' 
transferred to the ions and neutrals. Other sources of energy input 
to the electron gas, such as Joule heating,, are neglected in this 
work. 

The main loss of heat by the electron gas is energy transfer to 
ions. The electron-ion loss rate is given by 


_f: n n. [T -T.] , 

L . = 7.7 x 10”° = — %rr— eV cm"^ sec” [Banks and Kockarts, 1973] 

Cl A T '3/ L. 




(5). 
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where is the ion mass in atomic units and n. the density 

i U 

of the i^ ion species. The total loss rate to the. ions is the 
sum of the loss rates to each individual ion species. 

The ionospheric depletion produces a simultaneous reduction of 
the heat input and the ion cooling rate. The cooling rate drops 
faster than the heat input. The depression of both the ion and electron 
population has a larger effect on the electron to ion cooling than on 
the photoelectron heating. An overall increase in electron tempera- 
ture results. 

Direct energy transfer from the electron gas to the neutral gas 
is small in the unperturbed ionosphere. Collisions with neutrals 
become important in the ionosphere when the ion population is 
significantly reduced. Cooling may result from energy transfer to 
atomic oxygen as given by 

Leo = ^-2 X (V^n) 

Thermal electron collisions with the neutral gas may produce excited 
atoms which will release optical radiation (see Chapter 6). Monatomic 
oxygen also removes heat from the electron gas through fine structure 
transitions represented by the reactions 

e + O(^Pj) e 

The fine structure loss is written as 
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L = (5.92 - 4.68 x 10^^ T ) (9.06 + 6.57 x lO’^ T ) 

X 10"^^ ngn(o)' (Tg-T^)/T^ ergs cm"^ sec“ (8) 

[Herman and Chandra, 1969-] 

C, Thermal Structure of the Ion Gas 

The ion temperature value is bounded by the electron and neutral 
temperature values i The electrons provide a heat source and the 
neutrdls provide a heat sink for the ions. The ion heat source is 
identical to the electron loss given by (5). Other sources of heat 
input to the ion gas, such as electric fields and chemical reactions, 
can be neglected in the F-region ionosphere. 

The cooling of the ions is dominated by collisions with neutrals— 
mainly atomic oxygen. The ion-neutral heat transfer is given by 

L,n = 6 X 10-’'^ (T,-T„) 

Collisions between similar species of ions and atoms may result in 
cooling by resonance charge exchange. This reaction for 0 and 
atomic oxygen is: 


o’*' + 0 ^ 0 + 0 ^ 


(9) 


*r « 

The charge exchange cooling rate for 0 is given by 


Lee = "1% 2-1 » If’" 


(10) 
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The Ion temperature is found by equating the heat source term 
to the heat loss term. This is justified because the convection, 
compression and conduction terms in the ion heat balance equation are 
negligible. Equating the source and sink terms gives 



+ L. 


ce 


( 11 ) 


Taking 'to constant and numerically equal to 45, T^ 

is found to be 



T + (an/n„) T 
n ^ e n* e 

1 + {an^/n^) 


where a = 6 x 10 


6 


( 12 ) 


The ion collisional cooling rate decreases in the depleted region 
because of the reduction in ion density. However, the ion heating is 
reduced- to a greater degree because of the reduction in both electron 
and ion populations. The overall effect is a slight reduction in ion 
temperature in the depleted region. 


D. Summary 

The release of reactive chemicals, such as hydrogen or water 
vapor, into the F-region ionosphere enhances recombination and thus 
may cause a significant depression in the local ion and electron 
densities accompanied by a rise in the electron temperature. This 
increase in electron temperature is a result of the reduction in ion 
cooling of the electron gas in the depleted region. During the day, 
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the energy input into the electron gas decreases slightly, but thermal 
electrons are not effectively cooled by ions in the region of reduced 
plasma density because the heat transfer from electrons to ions is 
dependent on the -product of the ion and electron densities. The rise 
in electron temperature during the night is less dramatic because of 
the absence of photoeleptron heating. Also, the electron thermal 
conductivity decreases with reduced electron densities partially 
insulating the region of high temperature electrons from its surround- 
ings. 


104 



Chapter VI 

SIMULATION OF THE MODIFICATION EXPERIMENT 

The descriptions of the chemical, diffusive and thermal processes 
given in Chapters II, IV and V are combined to produce a model of the 
chemical ionospheric depletion. This model is used to simulate the 
temporal and spatial variations in plasma density and temperature result- 
ing from the release of reactive gases. Only, modification by the release 
of molecular hydrogen is discussed in this paper but perturbations pro- 
duced by other gases may be simulated by changing the physical parameters 
in the perturbation part of the model. 

A. Ionospheric Model 

The ionospheric model incorporates the most recent F-region theory. 
This model has evolved over the last decade with contributions from 
several, scientists. As more knowledge about ionospheric processes has 
been gained, the model has been updated. The present version of the 
ionospheric model is an extremely useful and sophisticated tool. 

The geometry used in the model is illustrated in Figure 1. The 
X and y axes lie in the local surface plane of the earth. The 2 axis 
points vertically upward. The magnetic field line through the ^2 
release point lies in the x-z plane and passes through the origin. The 
main plasma constituents along the field line are indicated in Figure 1. 

The interaction of the various functional blocks in the model is 
shown in Figure 2. The blocks are categorized according to the processes 
they describe — neutral, plasma, and thermal processes— and according to 
the region of space they occupy — low (100-500 km) region, middle (500- 
3000) region, and high (> 3000 km) region. Two blocks are required to 
describe the extraordinary processes in the perturbed ionosphere. 
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z (VERTICAL) 



Fig. 1. Three dimensional coordinate system 
used in the plasmaspheric model. 
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Fig. 2. Block diagram of the model of the mid-latitude ionosphere and protonosphere used 
in the simulation of the chemical modification experiment. 













starting from the high altitude region, the function of each block is 
discussed. 

Following a procedure introduced by Park and Banks [1974, 1975], the 

region above 3000 km is treated as a reservoir of plasma and heat. The 

temperature of the protonosphere above 3000 km is assumed to remain con- 

/ 

stant. This is the upper boundary temperature. The plasma content above 
3000 km varies as proton and electron fluxes flow in or out of this region 
The upper boundary .concentration is derived from this protonospheric 
content. 

The plasma density between 500 and 3000 km is determined from the 
0 charge exchange reaction 

0"^ + H J 0 + (1) 

*f + 

and from the upper H boundary condition and the lower 0 boundary 
condition. The electron temperature is calculated assuming constant 
heat flux from the upper and lower temperature boundary conditions. The 
model of the neutral atmosphere consists of monoatomic hydrogen, mono- 
atomic oxygen and helium. Details of the model between 500 and 3000 
km are given by Park and Banks [1974, 1975]. 

The 100 to 500 km region is an updated and expanded version of the 
model described by da Rosa [1966] and Waldman [1971]. The 0^ ion con- 
tinuity equation is solved by using the 500 km flux from the middle 
region and chemical equilibrium at 100 km as boundary conditions. Day- 
time photoionization of the neutral atmosphere and charge exchange with 
the 0^ ion forms the basis of the molecular ion solution. During 
the day, photoelectrons produced by solar EUV act as the main electron 
gas heat source. At night, there is no photoelectron production and 
the electron temperature is maintained by conduction from the protono- 
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sphere heat reservoir. Thermal electron cooling Is determined from 
electron to ion and electron to neutral heat transfer. The upper thermal 
boundary condition is heat flux. At 100 km it is assumed that the neutral 
and electron temperatures coincide. Plasma motion by an imposed neutral 
wind is included in all ion and neutral continuity equation solutions. 

Two foreign constituents will be present in the modified F-region— 

1) the released gas and 2) the ions produced by reactions with this gas. 
For the molecular hydrogen release, the foreign constituents are Hg 
and OH"^. 

The Hg continuity has been discussed in detail in Chapter III 
and, in the paper by Bernhardt, Park and Banks [1975]. The j’nteraction 
between the Hg gas model and the unperturbed ionosphere model Involves 
chemical conversion of the 0^ ion into the OH molecular ion. The 
low region neutral atmospheric density determines the diffusive expansion 
of the gas. 

The solution of the OH^ continuity equation is described in detail. 
This solution is typical of the procedure used in solving the other par^ 
tial differential equations in the ionospheric model. The equation 
describing the motion of the OH'*' ion including diffusive and chemical 
processes is obtained from Equation (46) of ChapterlV. 

3[^ = + k^CO'^'lEHg] - kg[0H'*']ng (2) 

where n = [OH'*'] + [O'^] +[N0‘*‘] + [o/] + [n/I. the electron density, 
w (the ion velocity) ^ 

the Imposed drift, D the diffusion coefficient for OH^ in the atmos- 
phere, Hp the oh'*' scale height, Tp the plasma temperature 
(Tg + T.), k^ and kg are the (0'*',H2) and (0H'*',e) reaction 
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rates respectively. 


The OH diffusion coefficient is 


D = D.^Cl + Tg/T.)sin^I 


( 3 ) 


where I- is the magnetic dip angle and is the ambi polar diffusion 

coefficient given by 


(oh") = 


[O2] 


[Ng] 

^N2Ti 


[ 0 ] 

Vi 


(4) 


T, 


fO!^ Yn > 
^2 - 


and 

Y ' 
0 

are: 

= 7.62 X 

lo'S 

cm~^ 

-1 

sec 

K-l 

= 7.42 X 


cm"^ 

-1 

sec 

K-1 

= 1.25 X 

loi® 

cm ^ 

-1 

sec 

K-1 


With the above substitutions, the OH continuity equation is 
written in the form 


M + A + Bn + V + P (5) 

where n is the OH^ concentration, and A, B, C, and P are time 
varying functions of the concentrations and temperatures at each point 
along a fixed field line. Finite differences are used to approximate the 
derivatives. The notation for values of the quantities in time and space 
will be a superscript k denoting the value at time kAt sec and a 
subscript m denoting position z = mAz + 100- km. 
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Consider the Taylor expansion about time t = (k + l/2)At sec: 


m m2 3t m 


( 6 ) 


also 


r,k+l/2 ^ „k+l M z^xk^-T 

m m “ ’ 


2 '3t^m 

Equating these two expressions for gives 


m 


.(7) 


m m 
At 


_ 1 r/9n.k+l , /9n-ikl 
■ 2 


( 8 ) 


The substitution of Equation (5) results in 


r„k+l k-i 

1 

At 2 


r ^2 


-|k+l 


2 

L 9z 

+ A— + Bn + Cn^ + f 
dZ 

■]. 

D-^+a|J+ Bn.+ cn^+ P 
• 3z"^ J 


m 


(9) 

The finite difference approximations to the partial spatial 'derivatives are 


1 3n^ k ^ m+1 ^m-1 

^3z^m 2Az ’ 


/3^n\k 
V oJ, 


9z 


2'm 


«k o„k , „k 
^m+1 ~ %-l 

(Az)^ 


(10) 


These expressions are used to remove the derivatives from Equation (9) 
giving 


nk +1 „k 
n - n ^ 
m m _ 1 


At 


- k+1 _^k+l . „k+l ok+1 

.k+1 Vl ■ Vl . m ,„k+l 

(Az) 2 2AZ ^Vl 

+ -t- ^ 

mm m ' m ' m m 


A '' 

4. m / k k ^ , ok „k , pk, k^2 , pk 

2Az ' m+1 m-1 ' mm mm m 


m-1 ' 


2n^ + rt , 
m m-1 


(Az)' 


(11) 
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This formulation of the numerical solution of partial differential equa- 
tions is due to Crank and Nicolson [1947]. 

k+1 

Equation (11) is quadratic in the future density variable n^^ . 

The equation is linearized by expanding the quadratic term in a Taylor 
series 


, k+1 2 

<"m > 


/ k,2 

(n ) + At — rz — 

' m at 


k 3" 


K St 


kv2 


(n )" + 2At n„ 
' m' m 


, k+1 k . 

b (n - 
k ' in m' 


At 


„ k k+1 , k 2 

"in - ‘"m> 


( 12 ) 


With this substitution. Equation (11) is written in the form 


k+1 k+"[ k+1 ^ 

"m-l “m ^ "m "m+l \ 


{1 < m < M - 1} 


(13) 


where a^, 6^, and 6^^ are functions of D^, A^, P^, 

k k k 

At and Az. is a function of n^^, %-!’ '^m+1 ’ 


C , P , At and Az. 
m m 




m m 


Equation (13) represents a system of linear equations which have 
unique solutions if two boundary conditions are specified. The boundary 
conditions for the OH^ continuity equation are zero ion concentration 
at 100 and 500 km. (i.e, nQ^^ and nj^^^ = 0) 

The system of linear equations has three unknown values at positions 
(m-l)Az, mAz, and (m+l)Az. Such a system is called tri diagonal as 
is evident when the linear equations are expressed in matrix form 
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^1 

^1 

0 

0 

0 


"l 


«1 

“2 

^2 

^2 

0 



k+1 

"2 



0 

“3 

^3 

^3 



* 


• 

• 



“m-2 

Pm-2 ^M-2 


• 


• 

0 

• 

• 

• . 

Vl ^M-1 


„k+l 

"m-1 


«H-1 


(14) 


The following efficient method is used for solving the tri diagonal 
system. Two sets of variables, s^^ and q^, are defined recursively 



^1 


Pi 


^m R 


^m 


ot s ^ 
m m m-1 


"m 


^m ~ %\-l 

Pm ■ Vm-1 


{2 < m £ M-2} 


(15) 


n T, n o • ’ • n^jn, are then calculated from 
m- 1 m-2 . 2 I 


n 


1^+1 = ~ ” M -1 ^ M -2 

M-1 3m_-[ - a^_-| 


„k+l _ „ ^ „k+l 

n ” Q "" s ^ 1 1 

m m m nrf] 


{M-2 < m < 1} 


(16) 


Further detail on the solution of tri diagonal systems can be found in 
Todd [1962]. 

The ionospheric model initially creates an unperturbed plasma and 
temperature profile along a magnetic field line. The field line is 
separated into steady state and time varying regions. The region is con- 
sidered to be in steady state where the recombination or the diffusion 
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time constants are greater than the time steps (At) used in the solution 
[Waldman, 1971]. ' This procedure insures numerical stability. 

The hydrogen gas model is used to simulate the release of a spheri- 
cal cloud of perturbing gas. The program time-s.teps are made small 
(At ~ 1 sec) and the initial transients in the chemical depletion are 
simulated. At later stages in the modification, -the time steps are 
increased to one minute. The simulation continues until the hydrogen 
is exhausted and the ionosphere has recovered. 

B. Results of the Simulation . 

The density and temperature variations produced by the release of 
into the daytime and nighttime F-region are now examined in detail. 

The first part of this section describes the effects of the release of 100 
kg of at an altitude of 300 km and at a geographic latitude of 40® N. 
The field line through the release point is at L = 3 with a dip angle 
of 70.5®. During the day, neutral winds are assumed to be zero. At 
night, a 100 m/sec meridional wind is used in the model. This night- 
time equatorward wind both lifts the nighttime F-layer and moves the 
released gas cloud with respect to the magnetic field lines. The second 
part of this section examines the effects of varying the altitude and 
amount of injection. 

Figure 3 shows the evolution of the night-time F-region profile 

along the field line through the hydrogen release point. The upward 

plasma flux at 500 km is indicated at the top of each profile. A deep 

hole is created within 1 minute after the injection, thus splitting 

+ 

the F-layer into two distinct layers. Increased loss of 0 induces 
an increase in the downward plasma flow from the topside ionosphere 
which helps to maintain the uppermost of these two layers. The lower 
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NIGHT-TIME ELECTRON CONCENTRATION 



Fig. 3. Time sequence of nighttime F-region profiles illustrating the effects of 
hydrogen injection. 



layer is cut off from the protonospheric flux by the diatomic hydrogen 
gas cloud blocking the magnetic flux tube. This lower layer, isolated 
from any plasma source, disappears after 660 seconds. 

The temporal variation in the altitude and maximum density of 
the lower layer may be used in determining the mechanisms 
maintaining the night-time ionosphere. If protonospheric flux 
is the only maintenance mechanism, the lower peak density will decay 
exponentially through chemical recombination. If there is appreciable 
night-time production, i.e., by corpuscular ionization, the lower peak 
density levels off at a certain point after the initial period of 
exponential decay. ObservatioVis of such behavior with an -incoherent 
scatter radar or an HF sounder should provide quantitative information 
on production and loss rates in the night-time ionosphere. 

One half hour after gas release, the ionospheric profile has almost 
fully recovered, but the peak density is somewhat lower than the initial 
peak density for two reasons: (1) the recombination rate is still 

slightly higher because of residual molecules, and (2) the plasma- 
sphere which has been drained by increased downward flux, cannot provide 
the necessary plasma to maintain the original layer. 

The formation of the persistant nighttime depletion is shown in 
Figure 4. Contours of constant electron concentration are plotted in 
the meriodonal. plane. In the northern hemisphere, north is to the right 
of the plots. As the ionospheric hole forms, the effects of the magnetic 
inclination and the neutral wind become evident. Tilting of the isopycni 
contours is produced by the depletion of plasma flowing along inclined 
magnetic flux tubes. In the later stages of depletion (times greater 

V 

than 600 seconds), the F-layer to the left (south) of the release point 
is more depleted than the F-layer to the right (north) of the release 
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Fig. 4. Meridional plane variations in the nighttime F-region electron concen- 



point. This is a conseqijence of the imposed N-S neutral wind of 100 m/s 

moving the gas cloud. After 4800 seconds the F-la^er is recovering 

from the effects of the release of 100 kg H 2 . 

Figure 5 shows a sequence of electron profiles for the same condi- 
tions as Figure 3 except for the daytime ionosphere with an overhead 
sun. The unperturbed profile at 0 seconds shows a large F2 layer above 

an FI layer. One minute after the release of H^, the F2 ppak is reduced 

by an order of magnitude while the FI layer is unaffected. As at night, 
two layers are formed above and below the ionospheric hole. The daytime 
lower layer does not disappear as it did during, the night because of the 
photoionization plasma source. The daytime increase in plasma tempera- 
ture immediately after gas release produces a slight increase in the upward 
plasma flux. Eleven minutes after injection, the level of this flux 
has decreased below the unperturbed flux level because the F-region 
depletion has been communicated up the field line. After one half hour, 
the shape of the F-layer has recovered but the peak density has not. 

The daytime meridional plane contours of electron concentration 
in Figure 6 show asymmetry and tilting resulting from the inclination 
of the geomagnetic field and nonuniformity of the neutral atmosphere. 
Becctuse the H 2 cloud diffuses more rapidly upward into regions of lesser 
atmospheric density, its shape becomes nonspherical . The nonsphericity 
of the gas cloud coupled with the inclination of the magnetic flux tubes 
produces greater depletion to the north as compared with the depletion 
to the south of the release point. As .with the nighttime depletion, the 
daytime F-region is still recovering from the H 2 release after 4800 
seconds. 

The difference between nighttime and daytime recovery is demon- 
strated by a comparison of Figures 4 and 6. The contours in the region 
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Fig. 5. Time sequence of daytime F-region profiles illustrating the effects of hydrogen 
injection. 
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Fig. 6. Meridional plane variations in the daytime F-region electron conentration 



of the ionospheric hole open downward at night (time 180 sec) and upward 
during the day (time 600 sec). The nighttime plasma flux from above 
pushes the ionospheric hole down whereas the daytime photoionization 
source tends to replenish the depletion region from below. 

Small changes occur in the nighttime electron temperature. (Figure 
7). The slight increase in T^ is a consequence of lower electron to 
ion cooling resulting from the plasma density reduction. Magnetic 
Inclination effects are apparent in the two dimension Isothermal plots 
in Figure 8. 

The rise in the electron temperature is more dramatic during the day 
than during the night. Examination of the daytime electron temperature 
profiles along the field line through the point of H 2 injection (Figure 9) 
illustrates the effect of the lowered ion cooling. The local electron 
temperature minimum at the altitude of the F-layer concentration peak 
(300 km) rapidly disappears as high temperature electrons produced by 
thermalized photoelectrons enter the depletion region. The rise 1n‘ the 
electron temperature produces a decrease in the downward thermal heat 
flux from the overlying protonosphere. The maximum electron temperature 
occurs about 10 minutes after H 2 release. At this time the temperature 
at 300 km has increased fay 1500 K, 

The spatial extent of the thermospheric modification is shown in 
Figure 10. Before the gas release, time 0 seconds, a slight meridional 
variation in the daytime electron temperature results from the latitudinal 
variation in the solar zenithal angle. One minute after H 2 release, a 
100 km wide region of temperature enhancement (Figure 10) coincides 
with a similarly sized reduction in electron concentration (Figure 6). 
After ten minutes, the width of the supra-normal temperature region has 
increased to 250 km. Later, the electron contours relax toward the 
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Fig. 7. The rise in nighttime electron temperature. 
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Fig, 9. Time sequence of daytime electron temper 
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unperturbed state but even after 2 hours, significant thermospheric 
perturbation exists (Figure 11). 

Coincident with a rise in the electron temperature is a decrease 
in the ion temperature. The heat source Tor the ions are the electrons. 

In the depleted ionosphere, the electron-to-ion heat transfer is reduced. 
The nighttime unperturbed ion temperature profile is shown in Figure 12. 
After the reactive gas release, the ions are cooled by heat transfer to 
the neutral atmosphere (Figure 13). The magnitude of the ion cooling is 
similar during the day (Figure 14 and 15). 

The ion temperature reduction has a significant effect on the topside 
ionospheric, H concentration. The reduction in the F-region 0 concentra- 
tion produces a reduction in the concentration under conditions of 
Isothermal chemical equilibrium because reaction of O'*’ and H is the 
main mechanism for production in the F-region. However, comparison 
of the profiles before and after nighttime H 2 injection reveals a 
decrease in H only in a small region at the bottom of the H layer 
(Figures 16 and 17). Throughout most of the layer, the concentration 

increases. This is due to changes in loss and diffusion rates pro- 
duced by ion temperature variations. The (H^,0) reaction rate is pro- 
portional to the square root of the ion temperature. Thus, the reduc- 
tion in ion temperature tends to produce an increase in H^. The 
diffusion coefficient is proportional to the ion temperature to the 5/2 
power [Schunk and Walker, 1970] resulting in a reduction in the diffusive 
transport of H through the depletion region. These two effects produce 
a bulge in the H profile in regions where production of H by 0 is 
not important. After 2000 seconds, the protonosphere has been suffi- 
ciently drained so that the high altitude concentration is noticeably 
reduced but the lower altitude bulge persists (Figure 17). 
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Fig. .12.. Nighttime unperturbed 
temperature profile. 
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Fig. 14. Daytime unperturbed tempera- 
ture profile. 


Fig. 15. Daytime ion cooling. 
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Fig. 16. Nighttime -and 0"^ profiles along the 

field line through the point of H2 real ease. 
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Fig. 17- Same as Figure 16 but at a later time. 
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The bulge occurs at higher altitudes in the daytime plasma- 

sphere. Since 0 is a more important source of low altitude H in 

the day than at night, the depression in the 0^ concentration pro- 

+ 

duces a depression in the H concentration up to an altitude of 

1200 km (Figure 18). The daytime enhancement occurs at 1800 km 

where the 0^ and concentrations are comparable. 2000 seconds 

after the gas release, the upward daytime plasma flux has increased 

the protonospheric content, resulting in a larger bulge in the H'*’ 

+ 

concentration at 1800 km (Figure 19). Note that the perturbed H 
profiles in Figures 16 through 19 are referenced to the unperturbed 
H"'' profile just before H 2 injection. (The unperturbed concentrations 
during the times indicated in the figures would be slightly less during 
the night due to the continual drainage of the protonosphere and would 
be slightly more during the day due to the upward plasma flow into the 
protonosphere.) In summary, the release of H 2 into the F-region 
produces 1) a reduction of at low altitudes as a result of the 
o'*’ depletion, and 2) an increase in at high altitudes due to the 
reduction in ion temperature. 

The depletion of the F-layer produces a relative reduction in the 
protonospheric content above 3000 km. A temporal sequence in the 
formation of the nighttime protonospheric depleted tube or trough is 
shown in Figure 20. The protonospheric content reduction is plotted 
against the meridional cross-field distance at the base of the tube 
relative to the field line passing through the point of release. The 
asymmetry in the trough is due to the magnetic field inclination and 
to the motion of the neutral atmosphere. Since the nighttime neutral 
wind blows from north to sourth in the northern hemisphere, the tube 
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Fig. 18. Daytime H and 0^ profiles along the field 
line through the point of H2 release. 
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Fig. 19. Same as Figure 18 but at a later time. 
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Fig. 20. Nighttime depletion of the protonosphere. The reduction in protono- 
spheric content from 3000 km to the equator is plotted against the 
distance between adjacent field lines in the meridional plane on the 
surface of the earth. 
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Fig. 21. Daytime depletion of the protonosphere 


is more depleted on its southern side. The central content reduction 

t 

of the trough approaches a value of 4.2 x 10^^ cm“^ {a 2.4% reduction 
in the unperturbed content) with a time constant of ~1 hour. The half 
width of the tube is 75 km at the surface of the earth. This corresponds 
to an equatorial tube half width of 390 km. 

The daytime protonospheric trough formation is illustrated in 
Figure 21. The tube depletion is skewed toward the northern -side 
because of the non-sphericity of the released gas cloud and the incli- 
nation of the field lines as mentioned earlier. The trough formation 
in the daytime protonosphere is less rapid than in the nighttime 
protonosphere because of the thermal plasma expansion produced by the 
rise in electron temperature. Since the daytime O’*" concentration is 
larger tha the ni httime O'*' concentration, more reacts in the 
daytime region near the release point. The diameter of the cloud, 
when a fixed amount of the injected gas has reacted with the ionosphere, 
is larger at night than during the day. Thus, the daytime trough is 
narrower and has a larger central content reduction than the nighttime 
trough. The central content reduction of the daytime trough illustrated 
in Figure 21 approaches a value of 8.8 x 10^^ cm“^ (4,9% reduction) with 
a time constant of -3 hours. The half width of the tube at the surface 
of the earth is 60 km corresponding to an equatorial half width of 
300 km. 

The protonospheric depletion may be controlled by varying the 
amount of gas released. Both the depletion level and the equatorial 
half width of the protonospheric trough increases with the release of 
larger quantities of H 2 , but the increase is not linear. Figure 22 
illustrates the tube content reduction and the tube radius as a function 
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Fig. 22. Nighttime protonospheric depletion 
produced by hydrogen released at 
300 km altitude. 



Fig. 23. Daytime protonospheric depletion 
produced by hydrogen released at 
300 km altitude. 
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of the amount of nydrogen released at an altitude of 300 km into the 
nighttime F-region. A tenfold increase in the amount of H 2 released 
approximately doubles the content reduct on and similarly doubles the 
tube radius. The variations in the size of the L = 3 daytime trough 
with the amount of injected at 300 km is shown in Figure 23. 
Figures 22 and 23 demonstrate that the daytime trough is narrower and 
more depleted than the nighttime trough. 

If multiple releases of at a fixed point in the depleted 
magnetospheric tube are possible, any desired protonospheric tube 
depletion may be achieved. If after each H 2 injection, the F-layer 
recovers, the protonospheric- content reduction will be the depletion 
produced by one release multiplied by the number of releases but the 
width of the depleted tube will be the width produced by one release 
alone. For example, consider the release of 500 kg into the daytime 
ionosphere. If all 500 kilograms are released at once, a 500 km tube 
with a depletion of B% will be fomed. However, if the gas is released 
as five 100 kg injections separated by 24 hours (to allow F-layer 
recovery), a 22% depletion tube with a half width of 300 km will be 
produced. In Chapter VIII, it will be demonstrated that a narrower 
and more depleted protonospheric trough traps a wider range of radio- 
waves. Thus, by multiple H 2 releases, radio wave guides through the 
protonosphere. may be created. Unfortunately the random movement of the 
depleted flux tube by ionospheric electric fields makes multiple 
releases of reactive gases into the center of the tube highly unlikely. 

The protonospheric depletion may also be affected by the position 
of the gas release. Figures 24 and 25 are graphs of the nighttime and 
daytime tube content reductions for releases of 100 kg H 2 at altitudes 
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Fig. 24. Nighttime protonospheric depletion produced 
by 100 kg H 2 released between 100 and 500 km. 



Fig. 25. Daytime protonospheric depletion produced by 
100 kg H 2 released between 100 and 500 km. 
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Fig. 26. Minimum nighttime F-layer peak density pro- 
duced by the injection of 100 kg H2. 



Fig. 27. Minimum daytime F-layer peak density pro- 
duced by the injection of 100 kg 
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ranging from 100 to 500 km along the L 3 field line. Even though 
the peak of the ionosphere is at 300 km, the maximum tube depletion 
occurs for 330 km releases at night and for 240 km release during the 
day. 

The height of the release also has strong effect on the F-layer 
peak density. Figures 26 (night) and 27 (day) show the minimum value 
of the ionospheric peak density .resulting from the release of 100 kg 
at varying altitudes. Again, the maximum effect does not occur for 
releases right at the F-layer peak. 

C. Summary - Stability of the Trough 

The coupled plasma spheric-thermospheric model has been used to 
simulate the atmospheric perturbations produced by the release of 
diatomic hydrogen into the F-region ionosphere. The gas injection 
' produces a reduction in the* 0^ and electron concentrations, a slight 
increase in the H concentration in the /O - H transition region 
and a relative reduction in the and electron concentration in the 
magnetosphere. The increase in the nighttime electron temperature is 
small in comparision with the increase in the daytime electron tempera- 
ture. The magnitudes of the slight reductions in ion temperatures 
during both night and day are about the same. 

Further investigation needs to be done concerning the stability 
of the protondspheric trough. If strong plasma gradients are created 
at the edge of the trough, electric fields may build up tending to 
cause plasma transport across the magnetic field lines and produce a 
kink in the trough. The problem of confining a neutral region inside 
a plasma is the dual of the well studied problem of confining a plasma 
inside a magnetic bottle for use in thermonuclear fusion. Thus, well 
established theories on plasma stability may be applied to this problem. 
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Chapter VII 


AIRGLOW EMISSIONS RESULTING FROM CHEMICAL DEPLETION OF THE IONOSPHERE 

In the previous chapter, it was demonstrated that the release of 
reactive chemicals into the F-region ionosphere produces 1) supra-normal 
electron temperatures, and 2) increased concentrations of foreign poly- 
atomic ions. These two disturbed conditions of the atmosphere can cause 
observable increases in airglow emissions. At the time of this writing, 
no measurements of enhanced airglow coincident with the release of deplet- 

t 

mg gases have been made but the theory behind airglow processes is well 
understood so that the results of this chapter can probably be confirmed 
by experimental measurements. 

A. Airglow in the unperturbed ionosphere 

Naturally occurring airglow is defined as the emission of radiation 
by gases that have been excited by solar electromagnetic radiation. Energy 
that is stored and re-radiated after a short delay produces dayglow. 
Nightglow results from a long delay in the release of stored energy. 

Airglow may be produced by the absorption of radiation into a 
molecule or atom producing a transition from its ground state to an 
excited state followed by re-emission and transition back to the ground 
state. If the re-emission occurs of a wavelength that is the same as 
the excitation wavelength, the excitation process is called resonance 
scattering. This is contrasted to fluorescence which involves re- 
emission at a wavelength that differs from the excitation wavelength. 

Excited particles formed by photodissociation of molecules may emit 
radiation upon transition to a lower state. Photodissociation of dia- 
tomic oxygen by the Schumann-Runge continuum (1350 - 1750 A) produces 
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the first excited state of monatomic oxygen 0(^D) which emits 6300 A radia- 
tion upon transition to the ground state 0(^P) [Rob1e et al., 1975]. 

Molecules formed by photodissociation may release energy by re-asso- 
datlon. The oxygen atoms formed by dissociation of 0^ may undergo 
3-body recombination during the night to produce airglow emissions. in 
the 2600 to 3000 A range. [Hunten, 1971]. 

Chemical energy is stored in a molecule by photoionization. This 
energy is released by the formation of excited atoms in dissociative 
recombination reactions. The decay of these excited atoms produces air-' 
glow. Consider the following dissociative recombination reactions; 

+ e" ^ O(^p') + 0(^P) + 6.85 eV 

or 0(^p) + (^0) + 4.9 eV 

or ^ 0(^P) + 0(^S) + 2.76 eV 

or ^ 0(^D) + OC^ D) + 3.01 eV 

or (O^D) + 0(^S) + .79 eV 

Laboratory measurements indicate that. the branching ratios for the produc- 
3 1 1 

tion of P, D and S oxygen atoms are 1 : .9 : .1 respectively [Zipf, 

1970]. Thus, for every O2’*’ ion annihilated, one ground state 0(^P) atom 
is formed, .9 excited 0(^D) atom is formed, and .1 excited 0(^S) 
atom is formed. States 0(^D) and O'(^S) emit radiation at 6300 A and 
5577 A respectively (Figure 1). 

The production of an excited state may be forbidden by the law of 
spin conservation. The dissociative recombination reaction 

■ N0‘^(^i:‘^) + e' ->.N(^S) + 0(^D) + .79 eV 

does not conserve spin and therefore is not a mechanism for the production 
of 0(^D). However the reaction 
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NO'^(^S^) + g“ N(^D) + 0(^P) + .38 eV 

satisfies the spin conservation law. ' Therefore, dissociative recombina- 
tion of N0^{^E^) is a source of N(^D) which emits 5300 A radiation upon 
transition to the ground state. 

The kinetic energy stored in particles may be transferred to other 
particles by collisions. The photoelectrons produced in photoioniza- 
tion reactions may excite atoms and molecules through inelastic colli- 
sions. Similarly, in regions of high temperatures, the tail of the 
Maxwellian thermal velocity distribution may contain sufficient electrons 
to produce significant excitation of neutrals. 

Once a particle is in an excited state, it may either radiate by 
spontaneous emission or be deactivated by transferring its energy to 
another particle. Thus, deactivation, or “quenching," of excited 
particles limits the amount of airglow produced. For example, 0(^D) 
atoms can be quenched by nitrogen molecules and by electrons. 


0(^0} + -»• 0(^P) + N 2 rate: k 2 = 7 x lO"^^ cm^sec”^ [Moxon, 1970] 

0(^D) + e” ^ 0(^P) + e“ rate: k^ = 2 x 10”^ cm^sec"^ [Seaton, 195 ] 

( 2 ) 

Many spectral lines are found in the F region airglow. In the rest 
of this section, only the 6300 A line attributed to 0(^0) will be con- 
sidered. This line is produced by photodissociation of O 2 S photoelec- 
tron impact with 0, dissociative recombination, and thermal electron 
impact with 0. The four components of the daytime 6300 A airglow are 
compared in Figure 2. The two components produced by Schuman-Runge photo- 
dissociation of O 2 and by photoelectron impact are from Roble [private 
communication, 1975]. The last two components in Figure 2 are results 
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Fig. 2. Daytime 6300 airglow components in the unperturbed F-region. 


of a simulation of the unperturbed daytime ionosphere using the model 
of the concentrations and temperatures described in Chapter VI. During 
the night, only the dissociative recombination airglow component is 
significant. Calculations of other airglow lines, such as 0(^S) 5577 A, 
N(^D) 5200 A, and 02(^S) 7619 A, are similar to the 0{^D) 6300 A line cal- 
culations described below. 

1 

The dissociative recombination component of the 0(D) 6300 A airglow 
is found by considering 1) the rate of polyatomic ion production by 
reaction of 0^ and neutrals, 2) the rate of dissociative recombination 
of these ions to produce 0(^D), 3) the quenching rate (Equation 5.2), 
and 4) the spontaneous emission rate given by 

0(^D) -> 0(^P) + hv rate = A 2 -j= .01 sec~^ 

The continuity equation describing the 6300 A emission resulting from 

1 + 

the production of 0(D) according to the O 2 dissociative recombination 
reactions given in Equation 5.1 is 

rflCO’"] - IC2[0('d)J[N2] - l<3[0(^D)]ng - 
where y is the number of 0(^D) atoms produced by the dissociative 

4 * 

recombination of one O 2 (y 1), 3-| = k-][ 02 ] is the recombination 

coefficient due to O 2 , k-j = 1.5 x 10 cm /sec, k 2 and k^ are 

specified in Equation (5.2) and Cpj^ is the photon volume emission 
-3 -1 

rate (photons cm sec ) resulting from dissociative recombination. 

Under steady state conditions, the volume emission rate becomes 

^DR ^ (k2[M2^ *^4^e^'^^21^ photons/cm"^sec"^ (4) 
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The volume emission rate for 0(^D) activated by thermal electron 
collisions is calculated in a similar manner. The excitation process 
is represented by 

0(3p) + Q* ^ 0(^D) + e rate = $21 (5) 

where e* is a high energy electron. 

Collisional excitation of 0(^D) is strongly dependent on temperature 
because of the exponential nature of the Maxwell thermal velocity dis- 
tribution. Electrons must have energies above 2 eV to induce transitions 
from the ground state, 0(^P), to the next highest level, 0(^D) (Fig. 1). 
The excitation coefficient is calculated from 


s 


12 


f 

''12 


V f(v) dv 


( 6 ) 


where Q ^2 velocity dependent collisional cross section for 0 

V is the velocity 
f(v) is the^ velocity distribution 
v ^2 minimum velocity needed to excite the atom. 

For a Maxwellian distribution this equation becomes [Seaton, 1955] 


where 



= 2 


x 10 ^ 


-(E2-E^)/kTg 


cm^/sec 


= minimum excitation energy = 2 eV 
= electron temperature 
= Boltzmann constant 


(7) 


The volume emission rate is calculated considering excitation, quenching, 
and spontaneous emission. 
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( 8 ) 


Gje s^2'^g[0]/[l . + (k2[N2l + k3n^)/A2^] 


Profiles of the daytime 6300 A volume emission rate due to dissocia- 
tion recombination, and due to thermal electron collision, 

are shown in Fig. 2. Under typical conditions, is greater than 

Ej^, but when becomes greater than 3000 K, Ej^ may be significant. 

Integration of the emission rate along the line of sight gives the 
observed brightness. The brightness is measured in terms of a columnar 
emission rate as given by 


I = y e(z)dz .-lO"® (9) 

0 

where I is intensity in Releighs. 

3 

e(z) is volume emission rate in photons per m per sec. 
z = distance along line of sight in meters. 

B. Airglow in the perturbed ionosphere 

Chemical energy stored in an injected gas may be released in the 
form of heat and light. Observations of airglow are a measure of the 
changes in chemical energy resulting from aeronomic reactions with the 
perturbing gas. 

Any of the reactive gases mentioned in Chapter II will produce 
polyatomic ions which, upon dissociative recombination, will produce 
enhanced airglow. In this section only OH^, the charge exchange product 

■f* 

of 0 and molecular hydrogen, is considered. Energy is released in the 
dissociative recombination process 
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H 2 + O"^ + H ■ RATE: 2.0 x 10“^ cm^ 

OH'^(V) + e“ H + 0(^P) + 8.77 eV 

-> H + 0(^D) + 6.81 eV RATE; kg 1 x 10“^ cn? s"^ 

H + 0(^S) + 4.6 eV (^0) 

Spin is conserved in all of these reactions. The branching ratios for the 
production of the three lowest states of oxygen from 0H{ r ) are unknown at 
the present time and are taken to be ^ ^ for 0(^P), 0(^D) and 

(O^S) respectively. 

An interesting form of deactivation results from the chemical reac- 
tion between 0 and The reaction with the oxygen in its ground state, 

0( P), does not proceed rapidly because the reaction is slightly 
endothermic 

0(^P) + H 2 OH + H - .08 eV rate: kg = 1.7 x 10~^^ cm^/sec 

( 11 ) 

The energy possessed by the 0(^D) state is enough to make the reaction 
exothermic [Schofield, 1967] 

0(^D) + H 2 OH + H + 1.88 eV rate: kg = 2 x 10"^^ cm^/sec. 

(12) 

While this is not quenching In the strict sense, it does remove excited 
oxygen atoms. Thus, the addition of molecular hydrogen to the ionosphere * 
layer produces a modified volume emission rate (compare with (4)). 


[8-|[0“'] + kg[0H^]ng/3] 


+ 

3 e 


Simulated nighttime profiles of along the field line through 

the point of H 2 injection are shown in Figure 3 for various times rela- 
tive to the time of Injection. After one minute, the peak nighttime 
volume emission rate has increased by a factor of 10. At this time, a 
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large bite-out in the emission profile results from the severe depletion 
of 0 around the point of release. As the 0 and H 2 concentrations 
diminish, the peak airglow emission level drops below the unperturbed 
emission level. After one half hour, the nighttime 6300 A emission rate 
atthe point of injection is one third the emission rate just prior to 
injection. 

The volume emission rate is enhanced in regions where H 2 and 0^ 
exist simultaneously. This is illustrated in the sequence of meridional 
plane emission maps shown in Figures 4 through 7. Figure 4 shows the 
6300 A airglow layer just before gas injection. One minute after injec- 
tion, a brighttorus around the field line through the point of injection 
is formed (Figure 5). The intensity inside the ring is low because of 
the chemical depletion of the 0^ along the central magnetic field tube. 
The torus expands with the H 2 gas cloud producing a "smoke ring" effect 
(Figure 6). The emission on the left (southern) side of the ring is 
stronger than the emission on the right (northern) side of the ring 

because, at night, an equatorward neutral wind is blowing the left 
side of the’ cloud into regions of undepleted 0^. 

The nighttime variation in the 6300 A intensity found by integration 
along the field line through the point of H 2 injection is shown in Figure 
8. The unperturbed nighttime intensity level is 18 Rayleighs in this 
example. 60 seconds after the H 2 release the intensity rises to a 
maximum of 120 Rayleighs. Because of the 0^ depletion, the intensity, 
level drops below the unperturbed intensity level reaching a minimum 
value of 11 Rayleighs after 1800 minutes. Finally, as the ionosphere 

recovers, the intensity level approaches the unperturbed level. 

The large regions of increased electron temperatures resulting from 

the release of H 2 during the day produce increased airglow. Using 
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Fig. 4. Unperturbed nighttime airglow layer. 
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Fig. 5. Bright ring formed one minute after injection. 
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Fig. 6. Expanding A6300 airglow ring. 
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Fig. 7. Asymmetry in nighttime airglow ring due to equatorward neutral 
wind. 
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Fig. 8. Temporal variations in intensity through the center of 
the expanding ring. 
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Equations 5.7 and 5.8 (modified to include deactivation), the thermal 
electron impact volume emission rate profiles are calculated. These 
profiles are compared with the daytime dissociative recombination air- 
glow profiles in Figure 9. One minute after release, the dissociative 
recombination 6300 Aairglow component dominates the thermal electron 
component. However, ten minutes later, after the and O"'’ have been 
exhausted and remaining electrons have been heated, becomes larger ' 

than 

The expanding "smoke ring" effect is seen in the daytime emission 
maps in Figures 10-13. The daytime ring is more symmetric than the night- 
time ring because no neutral wind was used in the daytime model. High 
altitude, field aligned increases in the emission rate are produced by 
high temperature electrons 6 minutes after gas release (Figure 14). 

As at night, the daytime airglow intensity increases rapidly after 
chemical release (Figure 15). The peak intensity after 60 seconds is 
due to dissociative recombination of OH*^. After this time, the dayglow 
decays at a slower rate than the nightglpw because of the contribution 

from colli sional excitation by the electron gas. Without thermal elec- 

« 

tron impact, the dayglow intensity level would decay below the unper- 
turbed level (Figure 15, dashed curve). 

C. Summary - Use of Airglow Enhancement Measurements 

The effects of ionospheric depletion are measurable from airglow 
fluctuations. The temporal and spatial variations in airglow emission 
are directly related to the physical processes producing and resulting 
from F-region plasma reduction. The 6300 A line intensity arising 
from the first excited state of monatomic oxygen is sensitive to vari- 
ations in the F-region temperature and polyatomic ion concentration. 
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Fig. 9. Perturbed daytime X6300 airglow components. 
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Fig. 10. Unperturbed daytime airglow layer. 
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Fig. n. Bright ring formed one minute after injection. 
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Fig. 12. Expanding X6300 airglow ring. 
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Fig. 13. Perturbed daytime airglow layer containing the effects of 
dissociative recombination and thermal electron impact. 
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Fig. 14. Daytime airglow modification containing 
high altitude perturbation by thermal 
electron impact and remnants of the 
expanding ring at the F-layer peak. 
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Fig. 15. Temporal variations in intensity containing dissociative 
recombination and thermal electron impact effects. 
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Fig. 16. Structure of the expanding airglow ring in the presence 
of external neutral winds or electric fields. 
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other lines In the airglow spectrum such as the 5577 A line 
0(^S), the 5200A line of N(^D), and the 7619A line of have 

different sensitivities to fluctuations in the state of the F- region. 
Simultaneous observations of variations in several spectral lines may 
permit isolation of the thermal and chemical airglow excitation mech- 
anisms. 

The morphology of the expanding airglow ring is an indication of 
the electric fields and neutral winds in space surrounding the point of 

gas release. The changes in the ring as observed from the ground are shown 

schematically in Figure 16. Without any externally imposed E-winds or 
fields the ring will be circularly symmetric around the gas release 
point. A neutral wind will move the injected gas cloud toward one side 
of the depleted flux tube resulting in an elongated airglow ring with 
maximum brightness in the direction of the neutral wind. An electric 
field will move the partially depleted plasma tube toward the edge of 
the injected gas cloud. Thus, a neutral wind will move the optical 

portion of the ring away from the gas release point in the direction of 

the wind and in contrast, an electric field will move the hole in the 
ring away from the gas release point in the direction of the ExB tube 
drift. By measuring the variations in the structure of the airglow 
ring with respect to the point of gas injection, the magnitude and 
direction of neutral winds and electric fields defonning the airglow 
ring may be estimated. 
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Chapter VIII 

INFLUENCE OF THE CHEMICALLY DEPLETED IONOSPHERE ON THE 
PROPAGATION OF ELECTROMAGNETIC RADIATION 

The artificially depleted ionosphere and magnetosphere will pro- 
duce variations in refractive index \?i?hich affect the propagation of E-M 
radiation. Changes in radio waves propagating through the modified 
plasmasphere are indicative of physical changes occurring in this 
region. 

A wide spectrum of radiowave signals is affected by the region 
of depressed plasma density. Section A describes the reflection of HF 
radio waves from the bottomside ionosphere in the presence of a trough. 
Section B describes the trapping of HF and VLF waves in a depleted 
plasma tube formed by the draining of the protonosphere into the hole 
created in the ionosphere. Trapped waves may interact with the plasma 
producing amplification' and/or particle precipitation. The transit 
time of a signal between two points in a depletion duct is related to 
the plasma density .along the duct. VHF and HF wave propagation through 
the F-region trough is examined in Section C. These waves suffer 
changes in amplitude, phase and polarization. Measurements of these 
changes may be used to estimate the time evolution of the depletion 
region. The ionospheric hole may be used as a lens for focusing 
extragalactic radio signals. 

A. Ground-based investigation of the depleted ionosphere by electro- 
magnetic waves. 

HF ray trajectories may be significantly affected in the vicinity 
of a trough. Generally, rays reflected from the sides of the trough 
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are defocused. Focusing occurs for rays. propagating through the 
depleted region (see section C). Radio waves with frequencies greater 
than the critical frequency at the trough center may escape through 
the ionospheric ho.le. Doppler frequency shifts in the returned HF 
waves indicate trough movement. Therefore, if an ionosonde is used 
to measure the structure of the chemically depleted ionosphere, angle- 
of -arrival and Doppler information may be necessary for proper 
interpretation of the ionograms. A theoretical study of reflection 
of high frequency ionospheric radio waves in the vicinity of a trough 
has been presented by Helms and Thompson [1973]. 

Incoherent scatter radars are valuable tools for studying electron 
densities, electron and ion temperatures, and plasma drifts in the 200 
to 800 km region. The power of radar signals scattered from the 
ionosphere is proportional to the electron density at the scattering 
point. The width of the received power spectrum is directly related 
to the electron-to-ion temperature ratio. The Doppler shift of the 
power spectrum is a measure of the plasma motion along the radar line- 
of-sight [Evans, 1969; Carpenter and Bowhill, 1971]. Thus, Thomson 
(incoherent) scatter radars may be used to determine density and tempera- 
ture perturbations resulting from the injection of reactive gases. 

B. HF and VLF Radiowave Trapping in Depletion Ducts 

This section provides a simple theory of ionospheric and magneto- 
spheric guided propagation involving depletion ducts. -This mode of 
propagation differs from HF propagation in earth-detached radio ducts 
[Carrara et al . , 1970]. Raypaths following field aligned troughs (as 
opposed to raypaths guided by the ionospheric layer parallel to the 
earth's surface) are considered. The self-focusing property of ducted 
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waves has been examined [Sodha et al., 1972]. Only the existence’of 
trapping based on comparison of the ray curvature and the duct curva- 
ture is considered here. 

A rarified plasma tube extending through the magnetosphere 
becomes a duct when it guides electromagnetic waves. In the guiding 
process, the raypaths are bent by the gradients of refractive fndex 
resulting from the variations of plasma density in the tube. A ray 
is said to be trapped if it tends to be bent toward the center of the 
tube when it Is propagating parallel to the tube. 

The trapping criteria are found by equating the radius of curva- 
ture of the magnetic tube with the radius of curvature of the ray at 
each point along the tube. If the ray radius of curvature Is smaller 
than the radius of curvature at the upper side of the tube, the ray 
is bent toward the center of the tube and is trapped. If the ray 
radius of curvature exceeds the upper side radius of curvature, the 
ray is not bent sharply enough to remain in the duct {Figure 1). 
Similar arguments are appropriate for^^bending of rays at the lower 
side of the tube. The transition between trapped and untrapped rays 
occurs for the marginally trapped ray which has a curvature equivalent 
to the duct curvature. 
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Since the trough is aligned with the geomagnetic field, the radius 
of curvature changes along the raypath. For a centered dipole approxi- 
mation of the magnetic field, the duct radius of curvature is given by 

R„ (4-3 

"o ■ 3 (2-R/Ro' 3 (l+sin^x) ‘ ’ 

where is the duct radius of curvature, R the geocentric distance 
to a point on the field line, the geocentric distance to the 
equatorial point on the field line,, and x the invariant latitude 
(Figure 2). The radius of curvature is maximum at the equator (x=0) 
with a value of R^/3 



Fig. 2. Centered dipole .magnetic 
field geometry. 


The ray radius of curvature is determined from the wave normal 
radius of curvature and from the relationship between the ray angle 
and the wave normal angle. The wave normal radius of curvature, 
determined from the differential phase velocity across the propagation 
path, is found to be 
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( 2 ) 


p = ^ = --P = -U 

/wn de dVp/dr dy/dr 

where is the wave noma! radius of curvature, dx the incre- 
mental distance along the path, d0 the change in wave normal angle 
along the path, Vp the phase velocity, dVp/dr the change in phase 
velocity across the propagation path and y the refractive index 
(Figure 3). 



The ray radius of curvature is defined as 


- ^ 
ray " d3 


(3) 


where dg is the change in ray angle. This expression is written in 
terms of the wave normal radius of curvature. 


^ray de dB ^wn dB 

The term d9/dB is the differential change of wave normal angle with 
ray direction and has the value of unity in isotropic media. 
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HF radio waves propagating, nearly parallel to the .magnetic field 
direction (i.e., quasi-longitudinal or ",QL" propagation) have identical 
ray and wave normal directions [Sudden, 1966], It. is clear that this 
condition of quasi -longitudinal (.QL) propagation applip.s to trapped HP 
radiation in field aligned troughs. For this condition, de/d 3 = 1 
and the ray radius of curvature is given by 

*^ray dy/ dr 

The QL approximation to the refractive index is 


2 

y 


= 1 



( 6 ) 


where X = f ts the wave frequency; f^ = ^-Vn^e^/Ce^m), 

the electron plasma frequency; n^ is the local electron density; 
e the electron charge; the dielectric constant of .free space; 
m the electron mass; = f^/f; fj^ = pos -the longitudinal 
component of the electron gyrpfr.equency; B is the imposed magnetic 
flux densi.ty; and the angle between the magnetic field and the 
wave normal. Jhe plus and minus signs denote the extraordinary and 
ordinary modes of propagation, respectively; Practical expressions 
for X and are 


X = 80.6 n^/f^ 

-{7) 

Y, = ’2.8x1 cos S/f 
L m 
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where the units for and B are electrons per meter^ and Webers 
2 

per meter , respectively. 

The refractive index is. further simplified by noting that is 
much less than unity for HF propagation in the magnetosphere. With this 
approximation. 


2 

P 


- 1 


- X 


(8) 


(The difference between extraordinary and ordinary mode trapping is not 
considered here.) Using (5) and (8), the HF ray curvature becomes 


P 


~ 2(1-X) 
ray - dX/dr 


(9) 


A Gaussian variation in the crossfield electron density is assumed 
and is given by 


2 2 

n^{r,R) = n^(l-d e“'" ) (.10) 

where n^(R) is the unperturbed electron density, d is the relative 

density reduction at the center of the duct, a{R) is the half width 

« 

of the duct, and r is the cross duct distance (r=0 at the center of 
the trough). A 500 kg release of H 2 at L=3 into the daytime ionos- 
phere produces an 8% depletion duct with an equatorial half width of 
500 km (Chapter VI). The numerical examples in this section are for 
a 500 km duct with a variable depletion factor. The n^.(R) function 
is based on a diffusive equilibrium model of the magnetosphere 
[Angerami and Thomas, 1964; Park, 1972]. Figure 4 shows the n^(R) 
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R (earth radii) 


Fig. 4. Electron concentration 

variation along field line 
at L=3.. 


variation along the L=3 field line. The half width of the duct, 
a(R), varies inversely with the square root of the magnetic flux 
density along the tube. The magnetic flux dnesity in the centered 
dipole model is 

B=Bo 


X/ («R/R„) 


where is the equatorial geocentric distance, R the geocentric 
distance to a point on the field line, and the magnetic flux 
density at R = R^. The variation in the trough width is then 



a (-R/R 
0 ^ 0 ^ 

{4-3 


( 12 .) 
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where is the equatorial half width. This variation is shown in 
Figure 5 for a^ = 500 km and = 3 earth radii. 



Fig. 5. Variation of duct radius with geo- 
centric distance. 

Equating the ray curvature and the tube curvature leads to an 
expression for the marginal trapping frequency as a function of 
depletion and position inside the duct. 

p 3ir 

( 13 ) 

where is given in (1) and (7) and {9) have been used. 

The ray curvature is maximized near the point of maximum electron 
density gradient inside the tube {f=a). Here, 


n^ = n^O-de"^/^) 
3r " "o a 


(14) 
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given these substitutions, the marginal trapping frequency becomes 


- \J 40.3 n^[2 + - 2)] (15) 

the HF trapping conditions as a function of tube depletion are shown 
in Figure 6 at four points along the duet. The trapping frequency 
regions are bounded on the high frequency side by the marginal trapping 
frequency and on the low side frequency by the plasma frequency below 
which there is no HF propagation i Note that at a constant level of 
depletion, the range of trapped frequencies increases as the geocentric 
distance to the duct, R, decreases. Also the minimum trapped 
frequency, the local plasma frequency, increases with decreasing R. 

The utility of the marginal trapping criterion is demonstrated 
by consideration of two HF rays with frequencies just below and just' 
above the marginal trapping frequency in a 60% depletion duct. The 
marginal trapping frequency for this duct is found to be 910 kHz from 
Figure 6. Figure 7 illustrates the path of-the ordinary, ray at a 
’ frequency of 900 kHz. The ray is launched at point A inside the duct. 
After an initially smooth trajectory, the ray oscillates inside the 
duct but always remains trapped. After entering the top side ionos- 
phere (point B) the ray reflects and is ducted back to the transmission 
point. The travel time between points A and B is .:14 second. For a 
ray at frequency 925 kHz, the equatorial portion of the duct is not in 
a trapping region (Figure 8) and the ray is not guided to the conjugate 
topside ionosphere. Note that the ray direction angles and the wave 
normal angles are identical along the path of propagation. 
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DEPLETION (7o) 






VLF propagation differs significantly from HF propagation. The 
propagation of the VLF "whistler"’ mode occurs for frequencies less 
than the electron gyrofrequency given by where e is the 

electron charge, B is the magnetic field flux density, and m is 
the electron mass. Only the ordinary mode propagates below the gyro- 
frequency. The VLF approximation to the ordinary mode refractive 
index may be derived from (6) with the result 


p 


A(cos - A) 


(16) 


where A is the ratio of plasma frequency to gyrofrequency, 9^^^ is 
the angle between the magnetic field and the wavenormal, and A = 

The refractive index is strongly dependent on the direction of propa- 
gation. The anisotropy greatly complicates the analysis of ducting of 
VLF signals. 

The relationship between the wavenormal angle and the ray direction 
angle for VLF propagation is most easily described by referencing these 
angles to the magnetic field lines. 


^m " 1 + cos (cos - 2A) 


[Hell 1 well, 1965] 


(17) 


where 8 and 9 are the ray and wavenormal angles with respect to 
m m 

the local magnetic field direction. 

For longitudinal propagation along a field aligned duct, 8^,, - 0- 
This occurs for two wavenormal angle conditions from (17). 


1 ) cos “ 2A 

2) sin 9„ = 0 

' m 


(18) 
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Trapping of waves at 0^^ = cos“^ 2A is called Type 1 trapping and 

at 0 =0 is called Type 2 trapping, 
m 

Smith, et al . [1960] and Helliwell [1965] have analyzed the duct- 
ing of whistlers with the use of refractive index surfaces. The 
refractive index surface is the curve traced out by a vector drawn 
from the origin whose length is p and whose angle with the ordinate 

axis is 0 . The angle between a normal to the surface and a vertical 
m- 

line is the ray direction angle 

Figure 9 illustrates the shape of the wave normal surface below 

and above half the gyrofrequency. If A < 1/2, two minitnums exist 

in the surface for 0 = ±cos~^ 2A. d0„/dg„ is less than zero at 

m m m 

these points. By examining the changes in the refractive index surface 

as a ray propagates to regions of differing plasma density, it can ‘ 

be shown that Type 1 trough trapping occurs for waves propagating with 

wavenormal angles near ±cos~^ 2A. 

A maximum in the refractive index surface (for A < 1/2) occurs 

at 0 =0. Here, d0 /dg is greater than zero. Examination of the 
m m m ^ 

effects of plasma variations on the refractive index surface suggests 
that Type 2 trapping of VLF waves in enhancements occurs for frequencies 
less than half the gyrofrequency and at wavenormal s aligned with the 
magnetic field. 

There is only one minimum at 0^ = 0 in the surface for frequencies 
above half the gyrofrequency. de^/dg^ is less than zero at this point. 
Type 2 trapping in troughs occurs because of this minimum. 

In order to analyze the trapping of VLF waves in ducts aligned 
with the earth's magnetic field, the ray curvature versus the duct 
curvature approach is used. The VLF ray curvature is a function of the 
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Fig. 9. VLF refractive index surfaces. Below the half gyro frequency, 
trough trapping occurs at large wavenormal angles around 
e=0G* Above the half gyro frequency, trough trapping occurs 
around 0=0. 
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plasma gradients, the magnetic flux gradients, the magnetic field 
line curvature, the wave frequency and the direction of propagation. 
Thus, the solution to the VLF trapping problem is more involved than 
the solution to the HF trapping problem. 

It is necessary to express the ray and wavenormal angles with 
respect to the magnetic field in terms of the same angles with respect 
to the horizontal coordinate axis. The differential relationship 
between these angles for propagation a distance dx along a field 
line is 


de = d0„ - 
m 

dx 

Po 

(19) 

dB = dB„ - 
m 

dx 

(20) 


where 6 and B are the wavenormal and ray angles referred to the 
horizontal axis, and is the local radius of curvature of the 
field line. 

The ray radius of curvature is found from (4) 


^ray 


^wn 


^ =n ' ten. 

dB “^wnydB 


dx/dB ) 
PQ / 


From (3), I s 


giving 


*^ray 





Using (3) and (20) 


(21) 


( 22 ) 
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(23) 


d3 

da = m — 

0 

1 - 

^0 

Substitution into (22) yields 


‘^ray ^wn 


de_ 


m 


da 


m 



Solving for the ray radius of curvature one obtains 


(24) 


Tp + 0 ) 


1 + 


\r\ 33 


m 


m 


(25) 


For trapping, IPpayl-lpQl upper edge of the duct and 

IP^ayUlpQ the lower edge of the duct. 

The relationship between de^^^ and d3^ is found by taking the 
derivative of (17) 


de. 


m 


1 + 2 cos 0^(cos 0^ - 2A) + {cos - 2A)^ 


~sin^ 0m + (cos^ 0^ + sin^ e^^) (cos 0^ - 2A) + (cos - 2A)^ 


m 


(26) 


The two conditions for longitudinal propagation (3^ = 0) require 


d0 

for Type 1 trapping cos 6^ = 2A then J 

m d3 

m 


-1 


-1 


12 2 
sin 0JJ, 1 - U 


0 5 A ^1/2 


(27) 
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de 

for Type 2 trapping cos e = 1 then 

Tn 2 


2 - 2A 
1 - 2A 


0 <: A C 1 


(28) 


The wavenortnal curvature is calculated from the refractive index 
using (2). 


-y 


-2A^ A (cos 6^ - A) 
' m / 


wn dp/dr 


dA 

dr 


A (cos -A) - ^ (cos - 2A) 


(29) 


where is the derivative perpendicular to the field aligned 

direction of propagation. Using the expression for A (the ratio of 
plasma to gyrofrequency) in terms of magnetic flux density and plasma 
density and the Type 1 and 2 propagation conditions, (29) simplifies to 


.and 


Type 1: 


Type 2; 


^wnl 


^wn2 


1 

dB 

1 

dne 

Bdr 

■ 2 

n^dr 



1 

dB 

1 



Bdr 2(1-N) " 2 n^dr 


(30) 


(31) 


The crossfield gradient of magnetic flux density is related to the 
field line radius of curvature. A magnetic flux vector can be decomposed 
into a scaler magnitude and a unit vector. 


^ = B t 


(32) 


where B = |^| and T is a unit vector tangent to the magnetic field. 
The radius of curvature of a field line is defined by the equation 
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(33) 


(V X f) X f = — 

Pq 

-J- 

where N is unit normal to the field line (i.e., N*T = 0). In a 
static, current-free region of space. Maxwell’s second equation 
reduces to 

V X B = 0 

Equation (32) is substituted in (34) to give 


V X 


T = -VB X f 

' B 


(35) 


This expression for the curl of f is used in (33). The resulting 
triple vector product is expanded yielding 


VB - (VB-T) T 
B 



The dot product of N with both sides-of (36) yields 


(36) 


->■ 


VB-N -1 





(37) 


The directional derivative VB-N is simply the cross field gradient of 
magnetic flux density dB/dr. Therefore, 


Po ” dB/dr 

This expression is used to remove the magnetic field dependence from the 
expressions for wavenormal radii of curvature. 


187 



-1 


(39) 



Note that for Type 1, longitudinal propagation is independent 

of frequency. 

The expressions for doyd^j^j, (27) and (28) and for (-39) 

and (40), are substituted into (25) to give 


Type 1 p^^^ 

0 S A £ 1/2 



(41) 



The conditions for VLF trapping can now be determined. First, 
consider marginal Type 1 trapping. Marginal trapping occurs when the 
ray just follows the field line. In (41) p^^^ = p^ if and only if 

0 or (l-AA-^) = 0 ‘ (43) 

e 

Therefore, marginal Type 1 trapping occurs if no plasma density gradient 
exists or if A = 1/2. 


iftft 



This is a very significant result. Type 1 trough trapping occurs 

dn 

for any trough > 0 . In the absence of a trough the ray 

naturally follows the field line when propagating with a wavenormal 
angle = cos ^ 2 A, The bending of the ray by magnetic field gradients 
just compensates for the curvature of the field line (see (37)). 

A trough, no matter how weak, reinforces the gradients by bending the 
Type 1 ray toward the center of the depletion duct. 

Type 2 trapping is considered. If the density gradients are zero, 
(•42) becomes 


^ray rTT ‘ 

Thus, with constant plasma concentration, the ray radius of curvature 
is less than the field line radius of curvature. Consequently, 
deflection from the upper edge of the duct is automatic. The problem 
in Type 2 trapping is deflection of the ray on the lower edge of the 
duct. 

The conditions for bending of a Type 2 ray at the lower edge of 
a duct are frequency dependent. By examination of (42), 


if A < 1/2 then 
or if A > 1/2 then 


e 

1-2A 

( i . € • , 

n^di- 

2(l-A)Pjj 


1-2A 


n dr 

P 

2(l-n)p„ 

^ 1 . e . , 


enhancement trapping) 


trough trapping) 


Note that- if -> 0 or A -+• 1 , infinite plasma gradients are required 
for trapping. The effects of density gradients on Type 2 trapping are 
exactly as predicted by the use of refractive index surfaces. 
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The conditions for trapping of VLF waves in ducts are summarized 
in Table 1 . 


Table 1 

Criteria for VLF trapping in field aligned ducts 

FREQUENCY 


0 < A < 1/2 1/2 < A 1 1 



Jm. 

1 

dn^ 

0 


<u 

T3 

O 

s: 

cos 

m 

= 2A 

n dr ^ 
e 


c 

o 

•r- 



trough 



fO 






rG 

Q_ 

o 

Type 

_2 

“"e ^ 

a-1 

. . 2A-1 

CL. 

6m = 

m 

0 

n.dr - 
e 


Hgdr - Pp(T-A.y 





trough 


In both Type 1 and Type 2 trappings, field aligned propagation 
occurs for A = 1/2 independent of plasma concentration. This is 
because the refractive index surface is maximally flat at half the 
gyrofrequency (picture the transition between the upper and lower 
diagrams in Figure 9). 

Using the Gaussian duct model evaluated at the point of maximum 
density gradient ((14)), the depletion required for trough trapping is 
found to be 

d > 0 for 0 i A < 1/2; Type 1 
1 /2 

« ^ (L-l)a 1/2 ^A^l; Type 2 

Above half the gyrofrequency, a Tower duct radius of curvature requires 
a higher depletion level for trapping. ’ 
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DEPLETION (%) 



EQUATORIAL DUCT HALF WIDTH: 500 km 



Fig. 10. VLF Trapping Frequencies. 



R (eorth radii) 



R (earth radii) 


TYPE 1 
TRAPPING 



Fig. n. Magnetospherically reflected whistler Fig. 12. Retrapping of a VLF ray in a depletion 

in a depletion duct: f = 7 kHz, duct: f = 7 kHz, depletion = 20%. 

depletion = 40%. 



R (earth radii) 


TYPE I 
TRAPPING 



TYPE I 
TRAPPING 


Fig. 13. VLF trapping for frequencies greater than half 
the equatorial gyro frequency: f = 20 kHz, 

depletion = 50%. 
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Strong depletion required for trapping near the equa- 
torial gyro frequency: fn = 34 kHz, f = 30 kHz, ' 

depletion = 70%. " 
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The lower and upper limits to the VLF propagation frequency are 
the lower hybrid resonance (LHR) frequency and the electron gyro- 
frequency, respectively. The LHR is given by 


1 

^LHR 


1 


'Hi 


+ r 


m 


'Hi 'He 


[Stix, 1962] 


(45) 


where f^^ is the ion gyrofrequency, f^^ is the electron gyro- 
frequency and f^. is the ion plasma frequency. A practical expression 
for the LHR frequency is 


Tlhr = 6.53 X 10° B 


5.26 X 10^^ + n 

J 

9.68 X 10^^ + n. 


11/2 


(46) 


where B and n^ are in MKS units. A more general expression for 
\hR the effects of heavy ions is given by Smith and Brice 

[1964]. 

The VLF trapping frequencies as a function of tube depletion are 
shown in Figure 10. Between the LHR and half the gyrofrequency, 
trapping is always possible. As frequencies increase above t|^/2, 
larger depletion is required for trapping. A trapped signal propagating 
below the local half gyrofrequency at low altitude will have the 
greatest probability of becoming untrapped in the high altitude 
equatorial regions of the trough where its frequency may be greater 
for f^/2. 

Examples of Type 1 ducted propagation in a trough are shown in 
Figures 11 and 12. Figure 11 illustrates a magnetospherically 
reflected whistler undergoing LHR reflection at points A and B. 
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The travel time from point A to point B is 1.35 seconds. Note that 
for Type 1 propagation the wavenormal makes a large angle with respect 
to the ray direction. 

The natural guiding properties of the geomagnetic field are 
illustrated in Figure 12. The ray leaves the Type 1 trapping region 
and then returns to the duct to be retrapped. Figure 12 illustrates the 
fact that the curves of trapping frequency versus depletion level are 
only useful as approximate criteria. The trapping of a specific ray 
at a given frequency and initial ray angle is most reliably determined 
by ray tracing techniques. 

Type 2 propagation near the equator is illustrated in Figures 13 
and 14. The wave is initiated below at low altitudes, enters 

a Type 2 region (f > f|^/2) and is guided along the lower side of the 
tr ugh into a Type 1 region. As predicted by the marginal trapping 
theory, higher levels of depletion are required for Type 2 equatorial 
trapping than for Type 1 equatorial trapping. 

As the wave frequency approaches the~ gyrofrequency, larger duct 
depletion is required for trapping. In Figure 14, the wave frequency 
is 90% of the gyrofrequency. A 70% depletion is required for trapping. 
The arrows in Figure 14 point to the lower side reflections in the 
trough. For cases of Type 2 trough trapping, the ray stays between the 
center and the lower edge of the trough. 

The guiding of HF radio waves in naturally occurring depletion 
ducts has been discussed by Muldrew [1963, 1967a,b], Muldrew and Hagg 
[1969], Nishizaki and Matura [1972], Sharma and Muldrew [1973, 1975] 
and Cerisier [1974]. Conjugate echo studies involving satellites 
Alouette 1 and 2, ISIS 1 and 2, and FR-1 suggest that these ducts are 
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small in comparison with the artificially created ducts previously 
described in this report. The natural troughs have diameters less 
than one ki.lometer and are on the order of several thousand km in 
length. The natural depletion level is about 1%. Below L = 1.7, 
propagation in natural ducts has not been detected [Cerisier, 1974]. 

The formation of natural ducts is attributed to electrostatic fields 
produced in the E-region and extending up into the magnetosphere 
[Sharma and Muldrew, 1975] or to thunder cloud electricity [Park and 
Helliwell, 1971]. 

Artificial ducts produced by chemical gas release have the 
advantage of precise control of size, location and depletion level. 
Artificial troughs with better wave trapping properties than found in 
naturally occurring ducts may be easily produced. 

As mentioned previously, signals propagating in depletion ducts 
do not usually reach the ground. VLF signals have large wavenormal 
angles at low altitudes that do not penetrate the bottomside ionosphere. 
HF trapped signals are usually reflected because the low trapping 
frequency does not penetrate the topside F-layer. However, in arti- 
ficially created troughs with steep density gradients, radiowave 
frequencies above the ionospheric critical frequency may be able to 
penetrate the F-layer and be trapped in the duct to be guided to a 
conjugate ground receiver. Thus, there exists the possibility of inter- 
hemispheric HF propagation by ducting through the magnetosphere. 

Valuable information may be obtained from HF and VLF ducting of 
signals that do not penetrate the F-layer. HF topside ionograms of 
ducted signals can be used to estimate the electron concentration in 
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the trough. Similarly, variable frequency whistler propagation in 
depletion ducts may be used to determine the electron concentration 
along the raypath as has been done for whistler propagatibn in enhance- 
ment ducts (see, for example. Park [1972]). 

The ability to produce waveguides in the magnetosphere makes 
possible a number of experiments involving wave- particle interactions. 
For certain relationships between a charged particle velocity and an 
electromagnetic wave frequency and phase velocity, strong interactions 
occur [$tix, 1962]. These interactions represent coupling between 
wave and particle energies. The energy exchange between streaming 
particles and waves has been suggested for the amplification and 
absorption of whistler signals (see, for example, Bernard .[1973]). 

The collision-free magnetospheric plasma is not in a state of thermal 
equilibrium and, consequently, contains high energy trapped particles 
with a wide spectrum of pitch angles. Interaction between electro- 
magnetic radiation and trapped particles may produce pitch angle 
scattering. Some of the scattered particles are no longer trapped and 
will precipitate along magnetic field lines to low altitudes where 
they will interact with the relatively dense D, E, and F-region 
atmosphere and may produce observable ionospheric phenomena such as 
enhanced airglow, heating or ionization. Observations at the base of 
the trough may indicate precipitation effects triggered by natural or 
manmade VLF propagation [Rosenberg et al., 1971]. 

Ground-based wave-particle interaction studies involving VLF 
transmission through the F-region into the magnetosphere have required 
signals with small wavenormal angles with respect to the magnetic field 
direction for penetration of the ionosphere and for guiding along 
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Fig. 15. Focusing of a vertically incident plane wave by the 
nighttime ionospheric depletion one minute after the 
release of 100 kg H 2 . The downward propagating wave 
has a frequency of 10 MHz. 
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Fig. 16. Same as Figure 15 but one half-hour after the 

release. ^ 
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natural, randomly placed enhancement ducts [Helliwell and Katsufrakis, 
1974]. The excitation of VLF emissions by transmissions from space- 
craft into the depletion ducts at known locations permits the study of 
wave-particle interactions with signals at all wave normal angles 
under controlled conditions. 

C- ‘ Propagation Through the Depletion Region - The Ionospheric Lens 

The variation in phase velocity in the ionospheric depletion region 
produces bending and phase perturbations in the rays passing this region. 
The ionospheric hole may be modeled as a lens of unknown diameter and 
thickness. The size .of the lens may be determined from the variations 
in phase and amplitude of signals propagating through the lens. 

The focusing property of the chemically modified ionosphere is 
illustrated in Figures 15 and 16. A plane wave generated above the 
ionosphere becomes distorted by propagating through the ionospheric 
lens. The refractive index increase in the depleted ionosphere pro- 
duces a decrease in phase velocity. Rays immediately below the iono- 
spheric hole are bent toward the central axis of the lens resulting 
in an increase in signal amplitude. Surfaces of constant phase are 
indicated by tick marks along the rays in Figure 15. A ring of infra- 
normal amplitude surrounds the focal point of the lens as is indicated 
by the reduction in ray density on either side of the central axis 
of the lens. Rays passing through the ionosphere far from the depletion 

N 

region are unaffected. As the diameter of the hole increases, the 
rays are not as sharply bent, and the focusing is not as pronounced 
(Figure 16). 
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The horizontal extent of the hole and the magnitude of the vertical 
depletion through the hole may be determined from the variations in 
phase and amplitude of a plane wave propagating through the depletion 
region. A Fresnel, diffraction approach is used to relate the size of 
the depletion region to the variations in phase and amplitude. The 
plane wave is assumed to propagate through the ionospheric lens with 
only a delay in the phase front. (This is a reasonable assumption ' 
for VHF propagation.) Just past the ionospheric lens, the phase per- 
turbed signal is decomposed into an angular spectrum of plane waves. 

Each of these .plane waves propagates to the point of reception where 
they are recombined into the received signal. This is schematically 
illustrated in Figure 17. 

A Gaussian phase delay will be used to simulate the effects of the 
ionospheric lens. A unit amplitude plane wave is represented as 

U(x,y,z,t) = exp[j(2TTZ/X - ut)] (47) 


where X is the wavelength and w is the frequency in radians/sec. 

(The complex time factor will be suppressed in future equations.) 

After propagation through a Gaussian phase lens, the complex amplitude 
becomes 

2 2 2 

U(x,y,z) = Gxp[j(2TTZ/X + ee'^^ h (48) 

where a is the horizontal radius of the lens and e is the phase 
perturbation (in radians) produced at the center of the lens. The 
origin of the (x,y) coordinate system is centered on the lens. Taking 
z =» 0 just past the lens. 
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2 2 

U(x,y,0) = exp[jee“’^ ] (49) 

9 9 2 

where r - x + y . The Gaussian phase model is both a reasonable 
approximation to the ionospheric hole and an easily manipulated mathe- 
matical expression. 

For wave frequencies well above the ionospheric electron plasma 
frequencies, e is related to the change in electron content produced 
by the injected gas according to 

, = , y« d, = 4^ ilk = (50) 

where AX is the change in fjj/f produced by the chemical depletion, 
f^ the local plasma frequency, f the wave frequency, k the wave 
number given by 2 tt/X, and AI the reduction in total electron content 
along the wavepath. 

The decomposition of a signal into its angular spectrum of plane 
waves is given by 

A , |-) _//u(x,y,0)exp[-j2H(^ X + f y)] dx dy (51) 

tmOO 

where a and 3 are the direction cosines of each component plane 
wave. Each plane wave is represented as 

B(x,y,z) = exp[j ^ (ax + By + yz)] (52) 

where y = y1-o^-^. The received signal at a distance z is found 
by summing the angular spectrum components taking into account the 
different distances each plane wave travels. This is given as 
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U(x,y, 2 ) = ■/ 


X X 


[Goodman, 1968] 


C53) 


If the angular spectrum of plane waves has a negligible contribution 
at high wavenormal angles, considerable simplification is possible, 
will be assumed that , |) is nonzero only if j-^l , |f|«l. 

This means that most of the energy in the signal is concentrated with 
wavenormal directions close to the z axis direction. Under these 
conditions, the square root in (53) is approximately given by 


^ 2 2 
^1-a -3 1 



(54) 


This is referred to as the Fresnel approximation and is discussed in 
further detail by Goodman [1968]. With this approximation, (51) is 
substituted into (53) and the integration over (^ , |-) is carried 
out to give 


U(x,y,z) = ff U(x,y,0)exp[jiT((x^-x)^ + (y^-y)^)/Xz] dx^ dy^ 


exp[j2Trz/A]/(jXz) (55) 

This is the Fresnel diffraction formula. 

For a Gaussian phase signal given by (49), the half width of the 
angular spectrum is approximately, ci = ^ > and must be much less 
than unity for the Fresnel diffraction formula to be valid. 
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Substitution of the Gaussian phase signal into the Fresnel dif- 
fraction formula and evaluation at x,y = 0 produces the following . 
infinite series [see Appendix D] 


CO St 

U(0.0,z) = exp[j(2Tiz/X +e}] E ^ 

s=0 l2TTaV ® 


(56) 

where P (x) = E and are Stirling numbers of the 

second kind as given by 



(57) 


The first few polynomials are 


Pq(x) = 1 

Pl(x) = 1 
P2(x) = 1 + X 


, P 3 (x) =1 + 3x + x^ 

P 4 (x) = 1 + 6x + 7x^ t 

Pg(x) = 1 + lOx + 25x^ + 15x^ + x^ 


2 

The series converges if (Aze/2Tra ) is less than unity. Assuming that 
the distance to the lens (z) is greater than the lens radius (a), this 
convergence criterion is compatible with the validity criterion for use 
of the Fresnel diffraction formula. 

The variation in the received signal on the lens axis is represented 
by 


U(0,0,z) = (58) 

Without the ionospheric hole, A = 1 and Acf = 0. Equating this with 
(56) gives 
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Fig. 18. Phase thickness [e(2ir radians)] of the iono- 
spheric lens as a function of measured relative 
amplitude and phase, (e = A<j> for 20 log-io(A)=0. ) 
Darkened curve represents phase and amplitude 
variations produced .in a 30 MHz wave propagating 
through the disturbed nighttime ionosphere. 
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Fig. 19. Horizontal size of the Ionospheric lens as a 
function of measured relative amplitude and 
phase. 
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(59) 



Because this equation is complex, it represents two equations with 
two unknowns. The unknowns are c, the central phase perturbation 

2ira^ 

produced by the hole, and the ratio of the lens circular area 

to the square of the classical Fresnel zone radius, 

2 

Contours of constant e and are plotted on a A4) versus 

A graph in Figures 18 and 19. The units of observed phase change are 
2 tt radians and of observed amplitude change are decibels. Except for 
low values of phase, e and A<}) are nearly identical. The deduced 
horizontal width of the ionospheric lens is strongly dependent on the 
measured signal amplitude. 

Also indicated in Figures 18 and 19 is the calculated variation 

in A and A(J) of a 30 MHz signal propagating vertically downward 

through the nighttime depletion region produced by the release of 100 

kg H 2 gas. These variations are also plotted versus time in Figure 20. 

Ten seconds after the H 2 has been released, the signal amplitude has 

increased by 1.5 dB and the phase has increased by 48 radians point 

A). As the ionospheric hole evolves, the amplitude decreases and the 

phase increases. A peak of 80 tt radians in the phase is reached at 900 

seconds after gas release (point B), After this time, the ionosphere 

recovers producing a drop in both phase and amplitude. 

Using Figures 18 and 19, or an equivalent numerical procedure, 

the variations in A(j) and A may be converted into variations in 
27ra^ 

£ and which then may be used to produce estimates of electron 

content reduction and ionospheric hole size (Figure 21). The usefulness 
of this method is limited by the accuracy of phase and amplitude 
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Fig. 21. Electron content reduction and ionospheric hole size deduced from the variations in Figure 20. 
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Fig. 22. The Refracting Radiotelescope. 
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measurements. Since the relative amplitude is larger during the 
initial stages of ionospheric modification, more accurate estimates 
of the horizontal extent of the depletion region are made at this 
time. 

In practice, the measurement of absolute phase is accomplished 
by comparing the received signal phase with a fixed phase reference. 

The reference signal phase may be supplied by the transmission point 
in the form of a microwave signal that is not .significantly affected 
by the ionospheric variations. Alternately, two VHP signals with 
different frequencies, and consequently different phase and amplitude 
sensitivities to ionospheric plasma variations, may be used. Since the 
extraordinary and ordinary modes of propagation in a magneto-ionic 
medium travel with different phase velocities, the change in relative 
phase between the two modes, and equivalently the change in polarization, 
has an information content equivalent to the absolute phase change of 
one mode. The theory outlined in this section may be easily extended 
to convert differential phase or polarization measurements, along with 
amplitude measurements, into estimates of the size of the ionospheric 
hole. 

Low frequency radio astronomy is limited by the reflection of cosmic 
radio waves from the topside ionosphere. Controlled ionospheric 
depletion will permit the penetration of low frequency extraterrestrial 
radio waves for reception in a radiotelescope. 

Low frequency radio astronomy is also limited by the radiotelescope 
aperture size. The half power beam width of a radio telescope is on 
the order of Va where X is the radio wavelength and a is the 
radius of the reception aperture. The construction of extremely large 
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Fig. 23. Frequency of the wave 
focused in the ground- 
based image plane. 



Fig. 24. Approximate beamwidth of 
the lens (A/a). 
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EXTRATERRESTRIAL 
RADIO WAVES 



F-ig. 25. Refracting-refl acting radiotelescope. 



dishes comprising reflecting radiotelescopes for reception of low 

frequency (long wavelength) radio signals is impractical. 

A refracting telescope using the ionospheric lens is proposed. An 

array of antennas is located in the image plane below the hole in the 

ionosphere (Figure 22). The output of each antenna is connected to a 

separate tunable receiver. The frequency of each receiver is varied 

keeping the image plane at the focus of the lens. 

Consider the ionospheric lens with the time varying parameters 

given in Figure 21. The focal length of the lens is found from 

» 

h = 27ra^/eA = aV/(40.3 AI) 

where h is the focal length, a the lens radius, f the radiowave 
frequency in Hz, and AI the reduction in ionospheric content through 
the center of the lens. Letting h = 300 km, the frequency of the extra- 
terrestrial plane wave focused in the ground-based image plane is given 
by 

f = 3480(AI)^/^ 

0 a 

The temporal variations in this frequency are shown in Figure 23. The 
beamwidth of the telescope is less than .2 of a degree (Figure 24). 

A greater degree of flexibility is achieved with the refracting- 
reflecting radio-telescope illustrated in Figure 25. The radio ray 
trajectories are initially bent by the ionospheric lens and then focused 
by a parabolic reflector. This telescope permits a wider range of 
focused reception frequencies than the simple refracting radiotelescope. 
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Some degree of directional variability is achieved by movement of the 
dish. The dish may also serve to cut down interference by blocking 
manmade radio signals. Existing radiotelescope dishes may be used in 
this imaging system. . 

The resolving power of the refracting radio telescope is limited 
by aberrations. in the ionospheric lens. Detailed studies of the lens 
imaging properties will hopefully indicate that the distortion produced 
by the lens is not serious. 
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Chapter IX 


CONCLUSIONS - RECOMMENDATIONS FOR AN IONOSPHERIC PERTURBATION EXPERIMENT 

With the advent of the Space Shuttle System, routine launches into 
the lower thermosphere will be possible. The Shuttle Spacelab can be 
used to release chemically reactive gases directly into the F-region. 

The large payload capability (30,000 kg in 380 cubic meters) of the 
Shuttle makes it an ideal vehicle for experiments involving controlled 
chemical modificat on of the ionosphere. Not only can the Shuttle 
initiate the chemical depletion of the plasmasphere but, it can be 
used to observe the effects of the depletion. Chemical releases may also 
be performed by rockets or satellites, but the Space Shuttle is the pre- 
ferred vehicle for multiple releases of large quantities of injected gas. 

This research has presented an analysis of the fundamental aspects 
of the chemical modification experiment. It has demonstrated that sub- 
stantial knowledge will be gained from observations of the effects of the 
reactive gas release. This research has also suggested conditions for 
the optimal experiment. 

Of all the reactive gases investigated, hydrogen is the most 
promising. H 2 has the highest depletion efficiency per unit mass and 
the highest resistance to photodissociation. Loss of H 2 by conden- 
sation is highly unlikely. Chemical dissociation of is simpler 
to describe than the multi-stage dissociation of triatomic injectants 
such as H 2 O. 

Two disadvantages in using hydrogen instead of heavier gases are 
1) the large volume of the gas in storage, and 2) the rapid expansion 
of the gas upon release. For a pressure of 1 atm and a temperature of 
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3 3 

15*C, TOO kg of H 2 occupies about TO m . This volume reduces to 

3 3 

n m at 100 atm pressure and to 1.4 m for liquid hydrogen. In com- 
parison, the volume of 100 kg of H 2 O vapor at one atm and lOO^C is 

3 3 

170 m and in liquid form, 100 kg of H 2 O occupies 0.1 m . The rela- 
tively low mass of the H 2 molecule results in rapid diffusive expansion 
and less localized depletion than would be produced by the heavier H 2 O 
gas. 

Initial modification experiments should be conducted with hydrogen 
because of the simplicity in describing depletion by diatomic molecules 
and because of the efficiency of H 2 when considering weight and 
condensation limitations. Sophisticated controlled depletion experi- 
ments can later be conducted by H 2 O injection. 

The optimum altitude of gas injection is dependent on the experi- 
mental goals. Maximum depletion of the protonosphere is attained by 
release approximately 50 km above the F-layer peak. The maximum 
reduction in the daytime F-layer peak density occurs for releases 50 km 
below the unperturbed peak. Other effects -such as electron temperature 
rise and airglow enhancement are maximized for chemical releases in 
the vicinity of the F-layer peak.. 

The release altitude may be chosen from orbital considerations. 

The ground track for a space vehicle at 276.6 km altitude with an 
orbit inclination of 50° is shown in Figure 1. Orbits at this altitude 
overlap themselves every 24 hours (17 rotations for a vehicle traveling 
from west to east) and are therefore ideal for multiple gas releases 
over a given geographic location spaced by intervals of time of one 
day or greater. Orbits with large inclination with respect to the 
ecliptic, cover a wide geographic latitude band. 


218 



ro 



Fig. 1. Overlapping ground track for an inclined orbit at F-region altitudes. The spacecraft 
is moving eastward with a velocity of 7.743 km/sec at an altitude of 276.6 km making 17 
orbits about the earth per day. Selected conjugate points are connected by dashed lines. 
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The minimum 'time between orbit crossings decreases with increasing 
latitude. Table I lists the latitude of and the time between the orbit 


TABLE I 


Crossing latitude 

Time 

between crossing 

+ 

6° 

11 

hours 19 minutes 

+ 

17.6° 

9 

59 

+ 

27.4° 

8 

33 


35.1° 

7 

00 

± 

41.1° 

5 

45 

± 

45.1° 

4 

15 

db 

47.9° 

2 

8 

+ 

49.7° 

1 

29 


crossings shown in Figure 1. The short time between the higher latitude 
crossings permits subsequent observation and/or remodification of the 
depletion region by releases from the orbiting space vehicle. 

The gas should be released either over a ground station capable of 
measuring the effects of the deplet on or over the geomagnetic conjugate 
of the station. Local effects produced by the gas release are F-region 
plasma depletion, temperature' change and airglow excitation. Conjugate 
effects include protonospheric and conjugate F-layer depletion. 

Several conjugate locations are indicated in Figure 1. It is 
difficult to find stations that exist over land' and intersect the 
orbital ground track as at points A and A'. Releases at an orbit 
crossing, points B or point C’, may produce noticeable changes in 
ground-based observations of the conjugate ionosphere (points B' or 
point C). If the conjugate point is over the sea (point D') only 
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local effects may be easily. measured (point D). Equatorial releases 
may be investigated at low latitude conjugate observatories such as at 
points E and E'. 

Modification of regions with high-geomagnetic-latitudes may be 
conducted from orbits covering the mid-geographic-latitudes. Consider 
the geomagnetic conjugate stations Roberval , Canada and Siple Station, 
Antarctica (Figure 1). Because of the tilt of the geomagnetic poles, 
an orbit with an inclination of 50° crosses only Roberval. The release 
of reactive gases over Roberval may permit modification of the plasma- 
pause and creation of high-geomagnetic-latitude whistler ducts. Obser- 
vation of the effects of gas releases may be conducted at Siple 
Station, Antarctica, far from the regions covered by the orbital ground 
track. 

The ducted wave portion of the modification experiment must be 
conducted from space-born vehicles because VLF or HF signals trapped in 
depletion ducts probably don't reach the ground. An orbiting vehicle 
could release the gas at an orbit crossing such as point D in Figure 1 
and on returning to the release point emit HF or VLF waves. The rapid 
round trip travel time .3 seconds) of the HF signal may permit 
reception of the ducted echo by the spacecraft. The long round trip 
time (- 3 seconds) of the VLF signal makes it likely that the spacecraft 
will move out of the duct before an echo is received. The effects of 
VLF trapped signals may be detected from high altitude satellites or 
from ground-based observations of precipitation produced by the 
scattering of trapped electrons. 

The nature of the modification, experiment is dependent on the time 
of day. Investigation of the recovery of the nighttime F-layer from 
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chemical depletion will provide information on mechanisms such as 
plasma flow from the protonosphere or corpuscular ionization. The 
nighttime enhancement of 6300A and other airglow lines is associated 
with the dissociative recombination chemistry involved in the modi- 
fication process. 

The analysis of the daytime experiment is complicated by photo-- 
electron and ion production by solar rays. The recovery of the daytime 
F-layer is governed by the mechanisms' that govern the nighttime layer 
plus the photo processes that are present during the day. The daytime 
airglow lines are not only produced by the dissociative recombination 
process but by' particle impact and photodissociation excitation 
processes-. 

Interesting effects may be observed by conducting the ionospheric 

experiment at dawn during equinox. Consider the conjugate points B 

and B’ in Figure 1. At dawn; point B' will be illuminated by the 

sun before point B. Photoionization at B' may produce a plasma 

flow along the magnetic field lines to point B. The importance of 

the dawn production at B' on the ionosphere at point B may be 

* 

determined by chemically depleting the ionosphere at B' just as it 
is illuminated by the sun. The seasonal dependence on plasma inter- 
change between conjugate ionospheres at dawn may be similarly determined. 

The ionospheric hole created during the modification experiment 
may be used as a lens in a refracting radio telescope. This telescope 
may be used to investigate previously unmeasured galactic radio 
sources in the 1 - 10 MHz frequency range. The nighttime release of 
100 kg H2 at the peak of the F-layer produces a lens with an aperture 
radius ranging from 20 to 120 km capable of focusing radio waves with 
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decameter wavelengths. The angular resolution of a radio telescope 
using this lens is better than’ 0.2 degrees. 

Sources observable by the refracting radio telescope are the 
"brightest" in the galaxy. The galactic brightness temperature 
reaches a maximum for sources with frequencies around 10 MHz. The 
nonthermal (synchrotron) component of decameter radiation from the 
galaxy reaches brightness temperature of 250,000 K. (For meter 
wavelengths, the brightness temperature on the nonthermal radiation are 
of the order of 10,000 K. ) 

The "aiming" of the refracting radio telescope in initial experi- 
ments need not be precise. The galactic disk has an angular thickness 
of about 5°. The position of the gas release should be between the 
galactic disk and the receiving antenna. The extraterrestrial signals 
can be recorded using a receiver sweeping the range of focused 
frequencies. The precise location of the lens may be determined a 
posteriori by optical measurements of the airglow ring surrounding the 
ionospheric hole. Further information on the location and the structure 
of the lens may be obtained from ionsonde or incoherent scatter radar 
measurements or from perturbations in satellite radiowave transmissions 
propagating through the lens. 

Narrow-beam observation of the sun at decameter radio wavelengths 
has been previously impossible using reflecting radio telescopes. With 
accurate positioning of a gas release between the sun (angular diameter 
= .5®) and a receiving aperture, high resolution pictures of the radio 
sun may be obtained. 

The expenditure of time and money involved in performing the iono- 
spheric modification experiment can be justified by the enormous amount 
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of information derived from the experiment. By observing the pertur- 
bations in the F-region, one can use the ionosphere as a laboratory 

j 

for measuring selected physical and chemical properties. Diffusion 
parameters may be determined from the evolution of the released gas 
cloud, the ionospheric trough, and the airglow ring outlining the 
depletion region. Chemical reaction rates may be determined from the 
magnitude of the ionospheric depletion and the intensity of the en- 
hanced airglow. Knowledge of thermospheric heating and cooling rates 
may be obtained from measurements of the ion and electron temperature 
fluctuations. 

Imposed winds and electric fields may be determined from obser- 
vations of the effects of chemical releases. Transport of the injected 
gas cloud by neutral atmospheric drifts produces motion of the 
airglow ring. Movement of the depleted plasma tube results in motion 
of the hole inside the airglow ring. Therefore, neutral winds and 
electric fields may be determined from observations of the airglow 
emissions. Measurements of these optical emissions may be taken at 
any time of day. Electric field estimation by tracking the motion of 
the depleted plasma tube may also be obtained by radiowave measurements. 

The depleted ionospheric and magnetospheric regions can be used 
in determining physical quantities not directly connected with the 
formation of these regions. The magnetospheric trough may guide radio 
signals for use in wave-particle interaction experiments. Some phenomena 
that can be studied using depletion ducts are amplification of VLF 
radio waves and scattering of trapped particles by these waves. By 
using the perturbed ionosphere as a lens, observations of extraterres- 
rial radio sources may be made. 
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Several areas of investigation are indicated as an extension of 
this work. Further study of the early effects of the injected gas 
sweeping out the background atmosphere is needed. The effects of 
depletion of the dusk and dawn ionosphere should be examined. The 
study of releases at low latitudes (with nearly horizontal field lines) 
or high latitudes (with field lines opening into the magnetotail) 
are of interest. Extending the mid-latitude ionospheric model to the 
conjugate ionosphere would be useful in investigating nonlocal effects 
of the depletion. An extension of this work is to investigate the 
injection of gases other than H 2 and to investiate airglow emissions 
besides the 0(^D) 6300A emission. Further investigation of the' 
depleted magnetosphere should include a stability analysis of the 
depleted plasma tube and the interaction between the plasma and waves 
trapped in the tube. Finally, the imaging properties of the ionospheric 
lens for use in a refracting radio telescope should be examined in 
greater detail. 
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Appendix A 

CALCULATION OF CHEMICAL PROPERTIES 

Most of the data in Table 1 were obtained from experimental 
measurements of chemical properties. In a few cases the reaction rate 
of a chemical process was inferred from rates of similar processes. 

The reaction rates in all the ion-molecule reactions in Table 1 
are based on measurements [Ferguson, 1973], but the heat released in 
these reactions is not. The heat of reaction is calculated from the 
ionization potential and heat of formation of each substance in the 
reacti on . 

As an example, consider the reaction 0^ + H 2 -»• OH'*’ + H. This 
reaction is subdivided into a series of steps each requiring the 
addition or subtraction of a known amount of energy. 

The ionization potential of 0 is 13.618 eV 

0"^ + H 2 + e ^ 0 + + 13.618 eV 

The heat of formation of H is -52.096 kcal/mole or -2.261 eV. 
Since H 2 is the normal state of hydrogen at 300 K this means 

1/2 H 2 H - 2.261 eV. 

Therefore, 

0'^ + H 2 + e^0+H + H + 13.618 - 2(2. 26T) eV 
= 0 + H + H + 9.096 eV. 


fRscmme page 
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The heat of formation of OH is -9.31 kcal/mole = -.404 eV. The 
heat of formation of 0 is -59.553 kcal/mole = -2.585 eV. 

The formation of OH from O2 and , H2 is characterized as 

1/2 O2 + 1/2 H2 ^ 0 + H - 2.585 - 2,261 

OH - .404 

or 0 + H ^ OH + 4.442 eV. 

Therefore, the exchange reaction between e", 0^ and H2 is’ 

0“^ + H2 + e” OH + H + 9.096 + 4.442 
= OH + H + 13.538 

*i" 

The ionization potential of OH is 13.17 eV. The overall energy 
released is 13.538 - 13,17 = .368 eV giving 

0"^ + H2 ^ OH’*' + H + 3.68 eV. 

This somewhat lengthy procedure has been used to calculate the heat 

released for all the reactions in Table 1. All heat of formation and 

ionization potential data are taken from Franklin and Dillard [1969]. 

The dissociative recombination reaction rates for NO2, OH , 

^ " 4 * ^ 

H2O , N2O and NH are based on rate coefficients for similar ions. 
Biondi [1973] has shown that electron-ion rates are dependent on ion 
complexity: the rates decrease monotonically for ions of increasing 

atomic number and for ions with a higher number of ato|ns. For 
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example, the diatomic ion, OH'^, is assumed to have a rate of 
-7 

1 X 10 because its atomic number lies between two diatomic ions 
with known reaction times 

(i^ate = 1 X 10“®) and nJ (rate = 2.2 x 10“^). 

Photodissociation rates are calculated by Integration of the 
product of absorption cross section times radiation density over the 
solar input spectrum [Banks and Kockarts, 1973 and Bates and Hayes, 
1967]. The photodissociation coefficients in Table 1 are for an 
overhead sun. Data at other times of the day must be scaled accord- 
ing to solar declination. 
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Appendix B 


H 2 O AND H 2 PHYSICAL DATA 


Parameter 

Symbol 

Units 

H 2 O 


Molecular radius 

m 

Angstroms 

2 

1.25 

Molecular mass 

m 

kg 

3.0x10"^® 

3.34x10“^^ 

Molecular volume in liquid 

"b 

m 

3.0x10"^^ 

1.11x10“^^ 

Specific heat ratio 

Y 

— 

1.41 

1.3 

Latent heat of vaporization U 

J/kg 

2.47x10'^° 

6.8x10^ 

Saturated vapor constant 

B 

m^/CK-kg) 

1.31xl0~^° 

1.48x10’^ 

Gas constant 

R 

J/(K-kg) 

462 

4155 

Specific heat: vapor 

c^ 

0/(K-kg) 

1607 

10200 

solid 

^2 

J/{K-kg)‘ 

4186 

6120 

Surface tension 

cr 

Nt/m 

.09 

.00231 

Critical point temperance 

To 

K 

. 374.2 

33.18 

Triple point temperature 


K 

273.16 

13.947 


H 2 O data from American Institute of Physics Handbook [1963] 
H 2 data from Scott [1 964] 


^HECEDUSfO 
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Atipendix C 


EQUILIBRIUM CONDENSATION DURING EXPANSION 

The time varying degree, of condensation is calculated for an 
expanding gas cloud in thermodynamic equilibrium. The specific gas 
volume V is a function of temperature T and degree of condensation 
X as given by 


V = (1-x) BT exp(U/RT) (1) 

where B is the gas dependent saturated vapor constant, and U the 
latent heat of vaporization, and R the universal gas constant. The 
first law of thermodynamics for a perfect gas is 

dQ = dE + PdV (2) 


For an isentropic expansion, dQ = 0. The internal energy differ- 
ential dE is written in terms of partial derivatives with respect 
to T, V, X. 


'I" "(ifh " (f)'’' '"> 


The change in internal energy with temperature at constant volume 
and degree of condensation is a weighted sum of the gaseous and solid 
specific heats. 


— 'l = 

9T/v,x \9T/v,x, solid \3T/v,x, vapor 


X Cj + (l-x)c, (4) 


^RECEDING PAGE BLANK NOT FILMED 
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where is the specific heat of the vapor at constant volume and 
Cg the specific heat of the condensate. (3E/9v)j ^ ~ ^ ^ 

perfect gas. The energy change with condensation is heat gained by 
the change in the amount of each phase less the heat lost by conden- 
sation. 


lv,T= 


(5) 


Finally, the perfect gas law gives 


P = RT/v = RT(l-x)/v 
y 


( 6 ) 


where the gas volume is related to the fluid volume v by the 

* 

expression v = (l-x)Vg. 

Making these substitutions, the adiabatic heat equation becomes 


[c-|(l-x) + C 2 x]dT + RT(l-x)dv/v - [U - (c 2 -c^)T]dx = 0 (7) 

Differentiate (1) to get 

dv = [(l-x)BT exp(-U/(RT^)] + (l-x)B exp(U/RT]]dt +' [-BT exp(U/RT)]dx 

( 8 ) 

Divide by v to obtain 


dv/v = [-U/(RT) + l/T]dT - dx/{l-x) 


(9) 
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Substitute for dv/v in (7) to get 


[c^O-x) + C 2 X + d-x)(RT-U)/T] - [RT + U - {c 2 ~c^) T] dx = 0 


00 ) 


Solve for dx/dT to obtain 


dx x(l+b^T) 

dT " (l-b^T)T ir-b^T)T 


01 ) 


where = (R + c^)/U and b^ = (c 2 -c^ + R)/U 

The solution for x subject to x = x when T = T is 

0 0 




" ° ® '’‘0 ^ t{ T,(b,T,-l) 


dT, (12) 


where F(T) = / dT, = log,[(b,T-l y (b,T^-l )2 t] 

0 


Integrate to get 


x(T) = 


(b^T-l)‘ 


— ^ T (T„-T) Ca,b,+1/(TT^)] 


+ (a^+b^) • log^d/T ) 


(13) 


If Xq is zero and the temperature when condensation begins, 
(13) becomes 
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x(T) - [(T„-T)[a,b,+1/(TT„)] + (a,rt,) logJT/T„)]j 


(- 14 ). 


Equations (1) and (14) specify the temperature and condensation degree 
for each value of specific gas volume. 
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Appendix D 

DERIVATION OF THE DIFFRACTION SERIES FROM THE FRESNEL DIFFRACTION 
EQUATION - SPECIFIC APPLICATION TO THE GAUSSIAN PHASE SIGNAL 

The Fresnel diffraction of a wave 1s given by 

03 CO 9 ? 

U(x,y,z) =/ /u{x-|,ypO) exp[jir((x^-x) + (y^-y) )/xz] dx^ dy^ 

—03 -03 

“exp[j27Tz/X]/(jAz) (1) 

where X is the wavelength, (x,y) the coordinate system parallel to 
the plane of reception, U(x,y,0) the complex amplitude of the signal 
at z = 0,x and U(x,y,z) the received complex amplitude. Because 
the complex exponential In the integral varies rapidly with x-| and 
V, for small values of z> this integral i's difficult to evaluate 
numerically. This function is written in terms of the convolution of 
the initial signal with the Fresnel kernel 

U(x,y,z) = exp[j2«/x] * u(x,y.O)_ (2) 

Use the convolution theorem of Fourier transform theory to obtain 


U(x,y,z) - exp[j27rz/x]T‘'‘[exp[-jXZTT(f^ + fy)]?U(x,y,0)] (3) 


where and S are the inverse and forward Fourier .transforms 

respectively. Expand the exponential inside the inverse transform to 


give 

U(x,y,z) =■ exp[i2Trz/x3T ^ 


f)(-jXZ7f)^(f^ + flf/5i JU(x.y,0) 
S=0 X y 


( 4 ) 
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The differentiation theorem of Fourier Transform theory is 


V^^yU(x,y) = j‘^[(-4TT)^f^ + fj)^ ?U{x,y)] (5) 


2 9 ^ 9 ^ 

where V.. .. = — ^ — 2 • application of this theorem results in 


x,y 


9x 3y 


uu rt 

U(x,y,z) = exp[j2irz/A] X) {jW4ir)^ U(x,y,0)]/s! (6) 

s=0 


2 2 

If U is radially symmetric (i.e., U(z=0) is a function of x + y 
2 

r only) then 


U(r,z) = exp[j2iTZ/X] 2 

s=0 



The received complex amplitude signal is an infinite sum of spatial 
derivatives of the signal at z=0. Only the first few terms of this 
series are significant for small values of z. 

If the initial complex amplitude signal is continuously differ- 
entiable, the multiple derivatives in (7) simplify at r=0 giving 


U(0,z) = exp[j2TTz/A.] X) (jWir)' 

s=0 


(dr^)' 


U(r,0) 


(8) 


A radially symmetric Gaussian phase perturbation is considered, 
U(r,0) = exp[j(f)(r)] where (j)(r) = e exp[-4^/2a^] (9) 


It can be shown that 


‘'2-7 [expCjee-'^ )] 


(dr‘) 


2a ' 


iee 


-r 


2/2a' 


00 ) 
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s ^ 

where Pg(x) = ^ and are the Sterling numbers 

m=l . . 

of the second kind [Abramowitz and Stegun, 1965]. satisfies the' 

following recursive equation 


3 


= m2<"> + 


n+1 


In closed form 


M 



m-k. n 

(- 1 ) 

k! (m-k) ! 


( 11 ) 


The observed complex amplitude at r=0 found from (8) and (10) is • 


u(0,z) = expCj(27rz/X+e)] 2 (J-) 

s=0 WaV ^ W^/ 


( 12 ) 
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